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Final  Report  for  Work  Unit  2300/FF/05 

This  task  addressed  basic  research  to  better  understand  laser  ablation  deposition  of  thin 
films,  nonlinear  optical  (NLO)  effects  in  poled  quartz,  and  the  waveguiding  properties  of 
Ge  doped  silica  planar  waveguides.  These  efforts  expanded  the  understanding  of  the 
interaction  between  light  and  matter  and  thus  advanced  the  envelope  of  the  technology 
base  of  NLO  materials  used  in  optical  device  applications.  Research  directly  supported 
Air  Force  Objectives  for  Optoelectronic  Materials,  as  outlined  in  the  current  research 
Technology  Area  Plan  for  Project  2305,  Electronics. 

The  research  was  organized  around  the  following  areas:  (1)  laser  ablation  deposition  of 
thin  films,  (2)  NLO  effects  in  poled  quartz,  and  (3)  the  waveguiding  properties  of  Ge 
doped  silica  planar  waveguides.  For  the  laser  ablation  effort,  various  techniques  were 
used  to  investigate  the  optical  losses  and  the  properties  of  the  laser  ablation  films.  For 
NLO  effects  in  poled  quartz,  two  possible  sources  of  SHG,  which  are  the  non-bridging 
oxygen  hole  centers  (associated  with  Si  and  Ge)  and  peroxy  bridge  structures  were  found 
to  be  supported  by  experimental  evidence.  Through  ion  implantation  which  creates  an 
internal  poling  field  in  quartz,  second  order  NLO  effects  were  introduced  into  otherwise 
amorphous  films.  For  the  waveguiding  properties  of  Ge  doped  silica  planar  waveguides, 
credibility  of  computation  efforts  were  increased  significantly  with  high  degree  of 
correlation  between  experiment  and  theoretical  predictions.  Results  for  these  topics  are 
summarized  in  the  next  several  sections. 

For  a  more  detailed  account  of  the  research,  see  the  attached  Work  Unit  publications  and 
presentations  listing. 
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Time-integrated  optical  emission  studies  of  plumes  generated  from  laser 
ablated  germania  glass 

Paul  J.  Wolf 

Frank  J.  Seiler  Research  Laboratory,  Materials  Physics  Division ,  2354  Vandenberg  Drive,  Suite  6H79, 

USAF  Academy,  Colorado  80840-6272 

(Received  3  February  1994;  accepted  for  publication  9  April  1994) 

The  optical  emission  from  plumes  induced  by  ArF  laser  irradiation  of  Ge02  was  characterized  as  a 
function  of  laser  fluence,  distance  from  the  target  surface,  and  ambient  02  pressure.  Dispersion  of 
the  light  emitted  by  the  plume  in  a  vacuum  revealed  emission  from  both  neutral  and  singly  ionized 
Ge  atoms  as  well  as  neutral  O  atoms.  The  spatial  variation  showed  that  the  ion  concentration 
decreased  exponentially  from  the  target  surface  while  the  neutral  atom  number  density  reached  peak 
intensities  at  distances  of  ^1. 5-2.5  cm  from  the  target  surface.  Interactions  between  the  plume 
constituents  and  the  ambient  molecular  oxygen  increased  the  excited  Ge  atom  and  Ge  ion 
populations  in  the  plume  and,  most  notably,  significantly  enlarged  the  excited  O  atom  concentration 
over  that  produced  directly  from  the  ablation  process. 


I.  INTRODUCTION  distances  comparable  to  typical  target-substrate  separations 

(^4-8  cm)  should  provide  information  on  the  dynamics  of 
Pulsed  laser  ablation  has  become  a  very  successful  the  laser  ablated  material  immediately  prior  to  collisions 

method  for  depositing  high-quality  thin  films  and  modifying  with  a  substrate  surface.  These  studies  are  particularly  im- 

surfaces.  The  congruent  removal  of  material  from  a  target  portant  when  reactive  gases  are  required  for  stoichiometric 

and  the  relatively  high  kinetic  energies  of  the  plume  constitu-  film  growth  since  the  presence  of  this  ambient  gas  can  affect 

ents  are  important  characteristics  of  this  method  that  can  in-  both  the  nature  and  the  energy  of  the  material  impinging  on 

fluence  thin-film  growth  processes.1  Although  much  is  the  film  growth  surface.  Reactive  gases  are  commonly  used 

known  about  the  properties  of  films  deposited  under  typical  to  maintain  the  oxygen  content  in  films  produced  in  ablation 

laser  ablation  conditions,  knowledge  is  generally  lacking  in  experiments  that  employ  oxidized  targets  such  as  YBaCuO.4 

understanding  the  dynamics  of  laser  generated  plumes  and  The  formation  of  metal-oxide  diatomic  molecules  has  been 

the  kinetics  that  govern  thin  film  growth  by  this  method.  The  directly  correlated  with  improved  film  stoichiometry  and,  in 

collection,  analysis,  and  interpretation  of  spectroscopic  data  general,  better  film  quality,5  and  the  distant  separation  be- 

obtained  from  plumes  created  above  a  surface  irradiated  by  tween  the  target  and  the  substrate  allows  for  potentially  more 

an  intense  laser  pulse  details  the  plume  composition,  the  en-  gas-phase  reactive  encounters  which  can  produce  these  spe- 

ergy  content  of  the  ablation  products,  and  the  dynamics  of  cies.  Several  optical  emission  studies  have  been  performed  to 

the  plume  species  which  can  subsequently  lead  to  an  im-  understand  the  ablation  of  high  Tc  superconducting  materials 

proved  understanding  of  the  laser  ablation  phenomenon.  In  that  ultimately  result  in  the  transfer  of  material  onto  a 

principle,  the  evolution  of  the  plume  can  be  theoretically  substrate.4'8  Although  the  bulk  of  both  the  temporal  and 

predicted  with  an  understanding  of  the  initiation  mechanisms  time-integrated  optical  emission  studies  have  been  devoted 

and  kinetic  processes  (electron  temperature,  electron  density,  to  high  Tc  superconducting  materials,  optical  emission  spec- 

ionization  state,  etc.).  troscopy  has  also  been  performed  on  other  laser  ablated  ma- 

The  spatial  and  temporal  evolution  of  transient  species  terials  such  as  carbon,9  dielectrics, 10,1 1  metals  such  as  Al, 2  512 

produced  during  the  laser-target  interaction,  such  as  excited  and  various  polymers.13 

atoms,  ions,  and  molecules,  can  be  monitored  as  a  function  Recently,  stoichiometric  Ge02  films  have  been  fabri- 

of  the  irradiation  conditions,  ambient  gas  pressure  and  com-  cated  on  ambient  temperature  substrates  by  ablating  Ge02 

position,  and  the  axial/radial  distance  from  the  target  surface  targets  in  150  mTorr  of  02.14  The  high  02  pressures  required 

by  investigating  the  optical  emission  produced  in  the  plumes.  to  form  stoichiometric  films  suggested  that  gas  phase  reac- 

Many  of  these  spectroscopic  studies  have  been  performed  tions,  in  addition  to  the  presence  of  energetic  atomic  species, 

near  the  target  surface  where  interesting  plasma  dynamics  may  be  important  for  incorporating  stoichiometric  quantities 

can  be  explored.  Recently,  Mehlman  et  a/.2’3  reported  results  0f  oxygen  in  the  film.  Thus  the  purpose  of  this  work  was  to 

of  a  study  in  the  vacuum  ultraviolet  to  ultraviolet  regions  at  study  the  dynamics  of  material  ablated  from  a  Ge02  target  in 

distances  less  than  1  cm  from  the  target  surface  aimed  at  the  presence  of  molecular  oxygen  to  correlate  the  ablation 

obtaining  insight  into  the  early  stages  of  the  laser-target  and  process  and  the  plume  composition  with  the  formation  of 

laser-plasma  interaction.  In  contrast,  probing  the  plume  at  stoichiometric  thin  films.  The  light  emitted  by  a  plume  after 

_  ArF  laser  irradiation  of  a  Ge02  target  was  collected  and  ana- 

“Present  address:  Dept,  of  Engineering  Physics,  Air  Force  Institute  of  Tech-  '^d  using  optical  emission  spectroscopy  (OES).  The  time- 

noiogy,  Wright  Patterson  afb,  OH  45433.  integrated  emission  spectra  were  analyzed  to  identify  the 
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FIG.  I.  Schematic  diagram  of  the  experimental  arrangement. 


composition  of  the  plume,  and  the  behavior  of  neutral  and 
ionic  species  was  mapped  when  variations  were  made  in  the 
0?  pressure,  laser  fluence,  and  distance  from  the  target  sur¬ 
face.  In  a  related  study,  Vega  et  al. 15  investigated  the  ablation 
products  of  a  Ge  target  as  a  function  of  02  pressure.  These 
experiments,  however,  were  performed  at  distances  less  than 
5  mm  from  the  target  surface  and,  therefore,  could  be  insen¬ 
sitive  to  possible  reactions  occurring  downstream. 

II.  EXPERIMENT 

Figure  1  shows  the  experimental  arrangement  with  the 
ablation  chamber  and  the  data  collection  system.  The  experi¬ 
ments  were  performed  in  a  vacuum  chamber  that  typically 
achieved  a  base  pressure  of  2X10  6  Torr  with  a  turbomo- 
lecular  pump.  An  ArF  laser  [\=  193  nm;  28  ns  full  width  at 
half  maximum  (FWHM)]  was  focused  at  normal  incidence 
onto  the  surface  of  a  GeCU  target  to  generate  the  plume.  The 
laser  fluence  was  varied  between  0.5-6. 5  Jem  during  the 
experiments  by  adjusting  the  combination  of  the  laser  energy 
and  the  spot  size  on  the  target  surface.  All  data  were  col¬ 
lected  with  laser  repetition  rates  ^1  Hz.  The  experiments 
were  performed  both  in  vacuum  (2X 10  6  Torr)  and  02  pres¬ 
sures  ranging  from  10  4  Torr  to  200  mTorr  which  spanned 
the  pressure  range  used  in  previous  GeOt  thin-film  growth 
studies.14 

The  target  was  mounted  on  an  adjustable  shaft  which 
allowed  translation  perpendicular  to  the  target  surface.  The 
optical  emission  was  sampled  at  several  locations  from  the 
front  surface  of  the  target,  z,  by  moving  the  target  away  from 
the  “0"  position  while  simultaneously  adjusting  the  laser 
focusing  lens  to  keep  the  laser  fluence  constant.  Distances 
from  0.2  to  5  cm  were  covered  in  these  studies.  The  “0” 
position  was  determined  empirically  by  translating  the  target 
until  the  collected  intensity  lost  features  associated  with  dis¬ 
crete  or  continuum  emission. 

The  interaction  of  laser  radiation  with  the  target  pro¬ 
duced  a  visible  plume  extending  normal  to  the  target  surface. 
The  optical  emission  was  collected  normal  to  the  plume’s 
main  expansion  direction  with  a  spherical  lens  positioned  a 
focal  length  (23  cm)  from  the  plume  axis.  A  15  cm  focal 
length  cylindrical  lens  subsequently  focused  the  light  onto  a 


l()()-/xm-wide  entrance  slit  of  a  0.5  m  spectrometer.  The 
spectrometer  was  equipped  with  an  linear  diode  array  (PAR 
OMA  III)  to  detect  the  spectrally  resolved  emission.  A  glass 
microscope  slide  was  positioned  in  front  of  the  spectrometer 
to  eliminate  second  order  uv  lines  from  being  recorded  at 
wavelengths  above  350  nm.  The  OMA  III  signal  was  re¬ 
corded,  stored,  and  analyzed  on  a  computer.  Several  visible 
and  UV  emission  lines  were  recorded  in  the  following  wave¬ 
length  regions:  260  nm^X^340  nm  in  20  nm  increments; 
420  nm^\^520  nm  in  20  nm  increments;  and  580,  600,  and 
780  nm.  These  spectral  regions  were  chosen  based  on  the 
known  spectral  line  positions  of  neutral  Ge  atoms  (Ge  l),16 
singly  ionized  Ge  atoms  (Ge  II),16  molecular  GeO,17  and  neu¬ 
tral  O  atoms  (O  l).lx  Ten  laser  shots  were  required  to  achieve 
good  signal-to-noise  ratios  in  these  experiments.  The  light 
collection  system  was  not  calibrated  for  absolute  intensity 
measurements  since  only  relative  intensities  were  important 
in  this  work.  The  intensity  of  each  spectral  line  was  obtained 
by  simply  measuring  peak  heights  after  correcting  for  base¬ 
line  variations. 

The  ablation  rate  of  Ge02  was  also  measured  as  a  func¬ 
tion  of  the  incident  laser  energy  density  by  irradiating  fixed 
targets  in  vacuum.  The  laser  fluence  on  target  was  deter¬ 
mined  by  measuring  the  total  energy  delivered  by  the  laser  to 
the  target  position  through  all  the  optics  and  dividing  this 
value  by  the  area  of  the  laser  spot  on  the  target.  The  crater 
depths  were  measured  with  a  surface  profiler  at  several  loca¬ 
tions  along  the  length  of  the  crater  to  obtain  averaged  values. 
The  material  removal  rate  was  subsequently  plotted  as  depth/ 
pulse  versus  laser  fluence. 


III.  RESULTS 
A.  Ge02  ablation  rate 

Germanium  dioxide  targets  were  irradiated  at  energy 
densities  between  0.5  and  7.0  Jem  \  Well  defined  craters 
with  10-15  fim  hole  depths  were  produced  on  the  target 
surface  using  50-80  laser  pulses  at  fluences  F  above  2 
Jem'2.  At  fluences  between  1-2  J  cm”2,  75-100  laser 
pulses  were  required  to  achieve  measurable  crater  depths. 
The  irradiated  areas  at  F<  1  Jem  2  formed  very  shallow 
craters  and  they  appeared  to  be  heat  damaged  on  the  surface 
as  opposed  to  other  irradiated  sites  produced  with  higher 
laser  fluences  which  clearly  showed  crater  formation.  The 
difficulty  in  measuring  the  area  of  these  indentations  over 
large  laser  spot  sizes  introduced  large  relative  errors  in  the 
mass  removal  rates  at  F*sl  J  cm'2.  Figure  2  shows  the  re¬ 
sults  plotted  as  crater  depth  in  /xm  per  laser  pulse  versus 
laser  fluence.  The  data  in  Fig.  2  were  fit  to  a  phenomological 
expression  for  the  ablated  depth  as  a  function  of  fluence  de¬ 
scribed  in  Eq.  (1): 

d  =  a~]  ln(F/Fthrcshold),  (1) 

where  d  is  the  crater  depth,  a  is  the  absorption  coefficient  at 
193  nm,  and  F,hrcshojd  is  the  ablation  threshold  fluence.  As 
shown  in  Fig.  2,  this  expression  produced  a  reasonable  fit  to 
the  data  with  an  absorption  coefficient  of  1.3X  105  cm'1  and 
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FIG.  2.  The  crater  depth  per  laser  pulse  as  a  function  of  laser  fluence  for 
GeOi  targets  irradiated  in  vacuum.  The  Ht  to  the  data  produced  a  ablation 
threshold  fluence  of  0.46  J  cm  '-. 


a  threshold  fluence  of  0.46  J  cm  “.  The  absorption  coeffi¬ 
cient  from  the  fit  was  within  35%  of  the  accepted  value  for 
Ge02  at  193  nm.'y 

The  ablation  threshold  determined  using  the  above  pro¬ 
cedure  was  corroborated  by  examining  the  optical  emission 
produced  from  laser  ablated  targets.  By  varying  the  laser 
fluence  on  the  target  and  monitoring  the  intensities  of  neutral 
Ge  transitions  using  the  304,  288,  265,  and  252  nm  emission 
lines,  the  minimum  detectable  atomic  emission  occurred  near 
0.5  J  cm'2. 

B.  Plume  emission  characteristics 

Time-integrated  emission  spectra  of  a  plume  generated 
by  irradiating  a  Ge02  target  with  a  laser  fluence  of  6.5 
J  cm”2  at  a  base  pressure  of  4X 10”6  Torr  are  shown  in  Fig. 
3.  The  transitions  associated  with  each  spectral  feature  in  this 
figure  are  listed  in  Table  1.  This  figure  illustrates  the  dramatic 
change  in  the  plasma  emission  as  the  collection  position  is 
changed  from  the  target  surface  to  5  cm  away.  At  z=0.2  cm, 
the  emission  consisted  of  features  associated  with  both  Ge  I 
and  Ge  II  atoms  superimposed  onto  a  continuum  background. 
The  ion  lines  appeared  very  broad  and  had  intensities  com¬ 
parable  to  the  neutral  line  intensities.  The  continuum  back¬ 
ground  significantly  decreased  and  the  Ge  I  emission  lines 
became  much  narrower  1  cm  away  from  the  target  surface. 
The  Ge  II  features  were  reduced  in  intensity  relative  to  the 
neutral  lines  and  remained  broad.  At  3  cm  from  the  “0” 
position,  the  emission  was  dominated  by  Ge  I  transitions 
with  a  small  contribution  from  the  ionic  species.  Finally,  the 
Ge  ll  features  disappeared  at  z=5  cm  and  the  plume  radiated 
solely  from  electronically  excited  states  of  Ge  I  atoms. 

The  discrete  emission  lines  observed  in  this  study  prima¬ 
rily  consisted  of  transitions  originating  from  electronic  states 
involving  neutral  and  singly  ionized  Ge  atoms.  Most  of  the 
neutral  transitions  originated  from  the  Ap2-Ap5s  configura¬ 
tion.  As  shown  in  Table  I,  however,  transitions  originating 
from  higher  lying  electronic  states  in  Ge  I  were  detected  in 


WAVELENGTH  (nm) 


FIG.  3.  Plasma  emission  spectra  of  material  ablated  from  a  GcO:  target  in 
vacuum  (4X  10  h  Torr)  at  F= 6.5  J  cm  :.  The  spectra  arc  labeled  according 
to  distance  from  the  target  surface:  (a)  z~ 0.2  cm,  (b)  z~  1.0  cm,  (c)  2=3.0 
cm,  and  (d)  2=5.0  cm.  The  transitions  labeled  a-h  arc  summarized  in  Table 
I.  Although  the  intensities  have  arbitrary  units,  the  relative  magnitudes  are 
preserved  in  the  plots  in  order  to  compare  line  intensities  at  different  loca¬ 
tions. 


these  plumes.  The  emissions  from  Ge  II  were  comprised  of 
the  following  transitions:  5s2S~5p2P{\  5p2Pl)-6s2S, 
5 p2P{)-5d2D,  and  Ad2D-Af2F{) .  Higher-energy  transitions 
associated  with  Ge  ill  could  not  be  observed  due  to  the  de¬ 
tection  limits  imposed  by  the  spectrometer  and  the  diode 
array.  Examination  of  the  spectrum  at  long  wavelengths  also 


TABLE  I.  Details  of  the  atomic  transitions'  associated  with  the  emission 
lines  depicted  in  Fig.  3. 


Label 

Species 

Transition 

wavelength 

(nm) 

Transition 

r 

l-  luwcr 

(cm " l) 

upper 

(cm  ') 

a 

Ge  i 

269.13 

4 p-  V.-Sv’P, 

557 

37  702 

b 

Ge  i 

270.96 

4 pl  V.-.Vv1/*,, 

557 

37  452 

c 

Ge  i 

274.04 

4 p2  lS„-4l/'P, 

16  367 

48  962 

d 

Ge  i 

275.46 

4 p2  ’P;-5.s'P| 

1410 

37  702 

c 

Ge  i 

279.39 

4 />'  l5’,|-6.v'P| 

16  367 

52  148 

f 

Ge  i 

282.90 

4  p1  'S(,-4(/'P| 

16  367 

51  705 

g 

Ge  li 

283.18 

4 pl  ;~4/2F5/2 

65015 

100  317 

h 

Ge  M 

284.55 

4/;:  :Ds  ;-4/:F7 

65  184 

KM) 316 

aFrom  Ref.  16. 


1482 


J.  Appl.  Phys.,  Vol.  76,  No.  3,  1  August  1994 


Paul  J.  Wolf 


OXYGEN  PRESSURE  (mTorr) 


100  t- 


1  r~  i  : 

o  j-  ,  , _ 

250  270  290  310 

WAVELENGTH  (nm) 


FIG.  4.  Plume  emission  features  recorded  at  three  laser  fluences:  (a)  0.5 
J  cm  (b)  1.0  J  cm  and  <c)  6.5  J  cm  :.  The  spectra  were  obtained  in 
vacuum  at  z  -  1  cm. 


revealed  emissions  originating  from  electronically  excited, 
neutral  atomic  oxygen  (0*1).  (Note:  the  *  symbol  will  be 
used  to  specifically  designate  electronically  excited  states.) 
Molecular  emission  from  GeO  was  not  detected  under  any  of 
the  conditions  utilized  in  these  experiments. 

The  main  spectral  features  of  the  plume  remained  the 
same  over  the  range  of  fluences  utilized  in  this  work.  That  is, 
the  nature  of  the  emitting  states  from  the  neutral  species  and 
ions  with  the  same  degree  of  ionization  did  not  change  with 
distance  from  the  target  surface  regardless  of  laser  fluence. 
However,  the  intensities  of  the  atomic  lines  were  sensitive  to 
the  laser  fluence  as  illustrated  in  Fig.  4.  No  discrete  ion  emis¬ 
sion  features  were  observed  at  F<  1  J  cm  “  regardless  of  the 
position  from  the  target  surface  or  the  pressure  of  the  ambi¬ 
ent  gas.  However,  neutral  emission  did  appear  at  laser  flu¬ 
ences  as  low  as  0.5  J  cm"2  but  the  emission  was  extremely 
weak  at  z^0.5  cm.  Since  no  ion  state  emission  was  observed 
at  these  laser  fluences,  the  electronically  excited  neutrals 
were  presumably  dissociated  directly  from  the  target  surface 
and  these  states  were  probably  not  formed  via  ion  recombi¬ 
nation  reactions.  Continuum  emission,  which  has  been  fre¬ 
quently  observed  in  laser  generated  plumes  near  the  target 
surface,  was  conspicuously  absent  during  ablation  in  vacuum 
at  F^\  J  cm  As  the  fluence  increased  from  0.5  to  1.0 
J  cm"2,  the  neutral  line  intensities  increased  by  a  factor  of  15 
while  a  much  slower  increase  (factor  of  2)  was  observed 
between  1.0  and  6.5  J  cm-2  for  both  the  Ge  I  and  Ge  li  emis¬ 
sion  lines.  Although  the  intensities  of  the  individual  atomic 
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FIG.  5.  Relative  intensity  behavior  of  the  Ge  i  269.1  nm  emission  line  (  +  1 
and  the  Get:  283.2  nm  emission  line  (Cl  as  a  function  of  the  oxygen 
pressure.  The  Ge  i  intensities  are  obtained  at  e  -3  cm  and  the  Ge  n  intensi¬ 
ties  are  recorded  at  z  =  \  cm.  Both  data  have  been  obtained  with  a  laser 
tluence  of  6.5  J  cm  These  trends  are  common  to  all  other  neutral  and  ion 
Ge  emissions  recorded  in  these  studies. 

lines  varied  with  laser  fluence,  the  relative  intensities  of 
spectral  features  remain  unchanged.  This  implied  that  the 
local  equilibrium  was  maintained  among  the  energy  levels 
regardless  of  the  laser  fluence. 

The  presence  of  an  ambient  gas  dramatically  affected  the 
plasma  characteristics.  Figure  5  summarizes  the  general 
trend  in  the  relative  emission  intensities  versus  02  pressure. 
The  neutral  line  intensities  in  this  figure  were  recorded  at 
z  =  3  cm  while  the  ion  data  were  obtained  at  z=l  cm  and 
both  data  sets  were  produced  with  a  laser  fluence  of  6.5 
J  cm-2.  The  intensities  of  both  Ge  I  and  Ge  II  lines,  which 
were  proportional  to  the  excited  state  number  densities  of  the 
respective  species,  remained  relatively  constant  over  four  or¬ 
ders  of  magnitude  change  in  02  pressure.  The  neutral  emis¬ 
sion  intensities  increased  by  a  factor  of  10  between  10  and 
100  mTorr,  and  subsequently  approached  a  constant  value 
above  100  mTorr  of  added  02.  Clearly,  the  higher  oxygen 
pressures  caused  a  considerable  enhancement  in  the  neutral 
populations.  The  ion  emission  behavior  was  different  at  02 
pressures  above  10  mTorr:  the  ion  state  intensity  monotoni- 
cally  increased  but  showed  no  indication  of  leveling  or  de¬ 
creasing  at  the  highest  pressures  used  in  these  studies.  The 
intensities  recorded  at  120  mTorr  were  a  factor  of  10  greater 
than  those  obtained  at  the  lowest  pressures.  Similar  trends 
were  observed  for  data  recorded  with  a  laser  fluence  of  1 
J  cm”2  for  both  Ge  l  and  Ge  II. 

The  dynamics  of  laser  produced  plumes  are  partially 
governed  by  the  backing  pressure  of  an  ambient  gas.20  In 
vacuum,21  the  plume  is  allowed  to  expand  freely  in  a 
forward-directed  pattern  with  a  (cos  0)n  functional  depen¬ 
dence,  where  n  typically  ranges  from  10  to  40.  However,  the 
propagation  of  the  plume  is  modified  even  at  moderate  back¬ 
ing  pressures  resulting  in  a  narrowing  of  the  plume  and  a 
subsequent  increase  in  a?  to  a  value  as  high  as  200.  In  this 
study,  the  addition  of  100-150  mTorr  of  02  forced  the  plume 
to  appear  constricted  and  more  spatially  confined  than  in 
vacuum.  This  effect  could  cause  the  number  of  emitters  per 
unit  volume  to  increase  within  the  observation  region  result¬ 
ing  in  an  apparent  increase  in  the  emission  intensities  at 
higher  CT>  pressures.  This  mechanism,  which  could  account 
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FIG.  6.  The  relative  line  intensities  of  Ge  I  and  Gc  n  vs  distance  from  the 
target  surface:  (a)  shows  the  Get  (312.5  nm  transition)  behavior  at  two 
pressures  (  +  )  vacuum  and  (O)  200  mTorr  of  0:;  (b)  illustrates  the  Ge  ii 
(602.1  nm  transition)  dependence  in  vacuum  (  +  h  /:>(0:)=I  mTorr  (A), 
P(0;)=84  mTorr  (O).  and  PtCM^nO  mTorr  (□).  A  laser  fluencc  of  6.5 
J  cm"2  has  been  utilized  to  obtain  these  data.  All  other  neutral  and  ion 
emission  features  follow  the  trends  depicted  in  this  figure. 

for  the  trends  displayed  in  Fig.  5,  was  tested  by  simulating  a 
change  in  the  pressure  gradient  in  the  plume  by  using  differ¬ 
ent  laser  spot  sizes  at  a  constant  fluence.  A  tight  laser  focus 
causes  a  high-pressure  gradient  across  the  plume  which 
forces  a  stronger  expansion  than  a  softer  focus  at  the  same 
fluence.  Thus,  focusing  the  laser  spot  more  softly  should  pro¬ 
duce  a  similar  effect  as  that  created  by  increasing  the  back¬ 
ground  gas;  namely,  the  shorter  extent  of  the  expansion 
should  keep  more  emitting  atoms  confined  to  a  smaller  vol¬ 
ume.  Under  these  conditions,  the  experiments  showed  that 
the  ratios  of  spectral  line  intensities  for  the  soft  focus  to  the 
tight  focus  increased  by  10%-30%.  Spatial  confinement 
could  partially  account  for  the  pressure-dependent  intensi¬ 
ties,  but  other  processes  uncovered  after  examining  the  be¬ 
havior  of  the  oxygen  atom  line  intensities  versus  02  pressure 
may  dominate  as  discussed  below. 

Figure  6  shows  relative  axial  intensity  distributions  at 
successive  distances  from  the  target  surface  at  two  pressures. 
In  Fig.  6(a),  the  Ge  I  line  intensities  recorded  in  vacuum 
reached  their  maximum  values  1. 5-2.0  cm  from  the  target 
surface  and  the  position  of  the  maximum  intensity  shifted  to 
2.5-3.0  cm  at  an  02  pressure  of  200  mTorr.  Although  the 
data  were  recorded  only  at  z^5  cm,  the  emission  was  ob¬ 
servable  to  at  least  7  cm  from  the  target  surface  at  all  pres¬ 
sures.  The  ion  emission,  however,  had  a  peak  intensity  near 
the  target  surface  and  exponentially  decayed  over  distance  as 
illustrated  in  Fig.  6(b).  In  vacuum,  the  ion  emission  was 
recorded  out  to  —0.5- 1.0  cm  from  the  target  surface,  where 
the  signal  levels  fell  below  the  detection  limits  of  the  experi¬ 
ment,  whereas  increases  in  the  02  pressure  resulted  in  detect¬ 


FIG.  7.  The  evolution  of  the  844.6  nm  atomic  oxygen  line  vs  distance  from 
the  target  surface  at  F~ 6.5  J  cm  (a)  illustrates  the  behavior  in  0:  while 
(b>  shows  the  effect  of  N:  as  the  ambient  gas.  The  gas  pressures  in  each 
figure  are:  (  + )  vacuum.  (A)  l  mTorr,  (O)  84  mTorr.  and  (□)  170  mTorr.  All 
data  have  been  normalized  to  the  intensity  measured  at  r  =0.2  cm. 


able  signals  as  far  out  as  z  =  3  cm.  Exponential  density  pro¬ 
files  have  been  used  in  gas  dynamic  simulations  in  the 
adiabatic  expansion  region  of  a  plasma,22  and  these  results 
may  support  this  picture  of  a  nonthermal  plume  expansion. 
The  laser  fluences  affect  the  intensity  distribution  only  to  a 
minor  extent;  the  brightest  regions  of  the  plume  appear  in  the 
same  location  and  to  the  same  extent. 

The  trends  observed  for  atomic  oxygen  are  contrary  to 
those  observed  with  Ge.  Figure  7(a)  shows  the  evolution  of 
the  844.6  nm  emission  line  from  atomic  oxygen.  In  the  02 
pressure  range  0  mTorr^PfO^)^!  mTorr,  the  O  l  emission 
rapidly  decreased  to  negligible  values  at  —3  cm  from  the 
target  surface.  The  O  I  emission  observed  under  these  condi¬ 
tions  was  clearly  a  product  of  electronic  excitation  of  atomic 
oxygen  due  to  the  ablation  process  because,  at  these  pres¬ 
sures,  the  relatively  low  collision  frequency  prohibited  any 
significant  interaction  between  the  oxygen  added  to  the 
chamber  and  the  plume  material.  At  higher  pressures,  how¬ 
ever,  the  O  I  emission  intensity  displayed  a  small  decrease 
but  attained  a  constant  value  in  the  region  1.5  cm^z^3  cm. 
The  0*1  populations  subsequently  decayed  at  z>3  cm  and 
finally  became  undetectable  near  z=5  cm.  The  behavior  of 
the  emissions  at  higher  pressures  could  be  caused  by  an  in¬ 
creased  residence  time  in  the  observation  region  as  discussed 
previously.  However,  the  40  ns  radiative  lifetime  of  this  tran¬ 
sition  (\S0-3P)  limited  the  travel  distance  to  =sl  mm  before 
decaying,  assuming  a  typical  velocity  of  10ft  cm  s“l  found  in 
laser  generated  plumes.  The  high-pressure  results  demon¬ 
strated  that  additional  0*1  production  mechanisms  must  exist 
to  account  for  the  observed  behavior. 
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In  order  to  better  understand  the  intensity  behavior  de¬ 
scribed  in  Fig.  7(a),  the  O  l  emission  was  monitored  as  a 
function  of  N->  pressure  to  check  whether  or  not  additional 
excited  oxygen  atoms  were  being  produced  in  reactions  be¬ 
tween  the  plume  material  and  the  ambient  molecular  oxygen. 
Figure  7(b)  shows  that  the  relative  0*1  number  density  re¬ 
mained  constant  at  a  fixed  location  from  the  target  as  the  N: 
pressure  increased  from  0  to  150  mTorr.  However,  the  O  I 
emission  decreased  exponentially  with  distance  with  the 
same  functional  dependence  at  each  pressure.  Thus,  the  pres¬ 
ence  of  02  acted  as  a  source  of  additional  atomic  oxygen  and 
this  trend  also  revealed  that  the  0*1  dependence  on  the  02 
pressure  was  not  a  residence  time  effect. 

IV.  DISCUSSION 

The  ablation  of  Ge02  glass  can  be  partially  explained  by 
considering  the  optical  and  thermal  properties  of  the  irradi¬ 
ated  material.  The  removal  of  material  from  a  target  by  laser 
irradiation  depends  on  the  coupling  efficiency  of  the  laser 
beam  energy  to  the  target.  The  absorption  coefficient,  a  of 
Ge02  at  6.4  eV  (193  nm)  is  2X105  cm"1  which  corresponds 
to  an  absorption  depth  of  -50  nm.19  For  “congruent  evapo¬ 
ration,"  the  heat  diffusion  distance  8  must  be  smaller  than 
the  absorption  depth.  We  can  estimate  the  heat  diffusion 
depth  using  typical  values  for  the  thermal  conductivity3  ~4 
(AT =0.015  W  cm" 1  K~!)  and  the  specific  heat25  (C/;=0.19 
cal  g  1  K"1)  for  glass  at  300  K.  On  a  time  scale  comparable 
to  the  laser  pulse  duration  (r=28  ns),  8[~(2Dt)v~]  is  cal¬ 
culated  to  be  —170  nm,  where  D-KlpCf)  and  p=3.68 
gem"3  for  Ge02  glass.25  Thus,  the  energy  absorbed  during 
the  laser  irradiation  time  is  confined  to  a  surface  layer  corre¬ 
sponding  to  the  thermal  diffusion  length. 

When  energy  densities  above  the  vaporization  threshold 
are  applied  to  the  target,  the  laser  energy  is  not  only  absorbed 
in  a  surface  layer  of  **50  nm  but  intense  local  heating  also 
occurs  within  a  distance  of  170  nm  which  can  lead  to  addi¬ 
tional  material  vaporization.  In  these  experiments,  the  laser 
ablation  threshold  was  near  0.46  J  cm  ~2  based  on  the  results 
of  both  the  optical  emission  studies  and  crater  formation  on 
the  target  surface.  The  surface  temperature  on  the  target  is 
estimated  at  this  fluence  using  a  simplified  approximation  to 
the  heat  diffusion  equation  in  one  dimension  described  by 
Ready26  given  in  Eq.  (2). 

T{z-{),t)  =  (2<t>0/AT)(Dr/  7r)l/2,  (2) 

where  <f>n=(<t>/r)  (1-/?),  <t>  is  the  laser  flux,  and  R  is  the 
reflectivity  (AM).  12—0. 13  for  Ge02  at  6.4  eV).1M  Using  the 
thermodynamic  values  defined  above,  the  upper  limit  to  the 
surface  temperature  at  the  end  of  the  laser  pulse  for  a  fluence 
of  0.46  J  cm-2  is  — 7000  K.  This  value  is  greater  than  the 
liquid  transition  temperature  (rliq^2625  K)  for  Ge02,  so 
removal  of  material  from  the  target  should  occur  within  the 
thermal  diffusion  length.  By  scaling  T(z  =  0,r)  with  <!>,  Eq. 
(2)  suggests  that  ablation  threshold  may  be  lower  than  0.46 
J  cm  2.  In  order  to  achieve  the  liquid  transition  temperature 
of  Ge02,  Eq.  (2)  predicts  that  the  laser  must  deposit  an  en¬ 
ergy  density  of  200  mJ  cm~~. 

The  interaction  of  the  laser  with  the  solid  surface  pro¬ 
duces  electrons  and  ions  in  addition  to  copious  amounts  of 


neutrals.  The  number  density  of  these  species  has  a  complex 
evolution  and  they  depend  on  the  irradiation  conditions,  the 
collisional  kinetics  in  the  plume,  and  the  distance  from  the 
target  surface.  The  most  important  result  of  this  study  is  the 
pronounced  effect  on  the  plume's  luminescence  and,  thus, 
the  number  densities  of  individual  species,  caused  by  the 
presence  of  molecular  oxygen.  At  02  pressures  below  10 
mTorr.  both  Ge*ll  and  0*1  decay  exponentially  from  the  tar¬ 
get  surface.  In  addition,  the  Ge  I  emission  peaks  when  the 
Ge*ll  and  0*1  concentrations  reach  their  minimum  values 
near  z  =  2  cm.  These  results  imply  the  existence  of  a  neutral- 
ion  conversion  mechanism  to  explain  the  delayed  peaking  of 
the  Ge  l  emission  intensity.  At  these  low  pressures,  the  equi¬ 
librium  mean-free  path  (MFP)  between  02  and  Ge  is  >1  cm 
so  collisional  interactions  between  these  species  or  other 
plume  constituents  with  the  ambient  could  be  eliminated. 
Therefore,  electron-ion  or  electron-atom  interactions  within 
the  plume  may  explain  the  increase  in  the  Ge*I  number  den¬ 
sity  at  the  expense  of  the  ions  and  atomic  oxygen. 

Both  the  neutral  and  ion  state  concentrations  begin  to 
increase  near  20—30  mTorr  of  added  02.  This  threshold  also 
coincides  with  an  increase  in  the  0  l  emission.  At  these  pres¬ 
sures,  the  MFP  is  typically  1-2  mm;  thus,  the  increase  in  the 
neutral  atom  concentrations  appear  to  be  linked  with  higher 
collision  frequencies  which  makes  the  reaction  Ge(l  or 
n)+0*l— >Ge*(l  or  n)-fO  more  probable.  The  ion  emission, 
however,  remains  confined  to  within  about  1-1.5  cm  of  the 
target  surface  where  electron  densities  are  higher  which  pos¬ 
sibly  promote  c”-atom  ionization  processes  to  account  for 
the  increase  in  the  ion  state  populations. 

The  ablation  of  Ge02  at  high  ambient  pressures 
[P(O2)>10  mTorr]  not  only  produces  0*1  directly  as  a  result 
of  the  ablation  process,  but  also  increases  the  excited  O  atom 
concentration  through  plume-02  interactions.  The  magnitude 
of  the  cross  sections  for  e"-02  collisions  resulting  in  the 
dissociation  of  02  into  various  excited  atomic  products  range 
from  2.2X10" 17  to  6.6X10”17  cm2  for  electron  energies 
ranging  between  13.5  and  33.5  eV.^7  These  values  translate 
into  relatively  efficient  collisional  dissociation  if  the  electron 
kinetic  energies  are  in  the  appropriate  range.  The  ablation  of 
targets  containing  Ge  is  known  to  produce  atomic  species 
with  average  kinetic  energies  of  40-100  eV.  i5  ”x  The  accom¬ 
panying  electrons  are  expected  to  have  higher  kinetic  ener¬ 
gies  than  the  atomic  species  which  makes  02  dissociation  by 
electron  impact  highly  probable. 

The  maximum  emission  intensity  is  determined  by  the 
collisional  excitation  mechanisms  and  the  radiative  lifetimes 
at  the  higher  ambient  gas  pressures.  Collisional  interactions 
between  the  plume  and  the  ambient  02  clearly  increase  the 
concentrations  of  excited  Ge  I  and  Ge  II  over  those  observed 
at  the  low  ambient  pressures.  In  addition  to  e~-atom  colli¬ 
sions  which  can  produce  more  ions  in  the  plume,  excited  O 
atoms  can  interact  with  Ge  atoms  to  produce  electronically 
excited  Ge.  The  presence  of  an  0*1  creation  process,  there¬ 
fore,  makes  this  event  more  probable.  The  excited  species 
are  also  observed  farther  from  the  target  surface  at  high  pres¬ 
sures  which  implies  that  the  processes  responsible  for  pro¬ 
duction  of  these  excited  states  have  a  spatial  dependence. 

The  results  of  this  study  impact  the  understanding  of 
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stoichiometric  Ge02  thin-film  formation.  Different  deposi¬ 
tion  conditions  typically  lead  to  variations  in  plasma  lumi¬ 
nescence  and  the  details  of  the  emission  spectra  can  identify 
conditions  for  producing  high-quality  thin  films.  In  an  earlier 
work,14  we  found  that  stoichiometric  Ge02  thin  films  could 
be  produced  at  high  ambient  02  pressures  (150  mTorr).  Fur¬ 
thermore,  these  studies  showed  that  stoichiometric  films 
were  produced  only  when  the  substrate  was  physically  lo¬ 
cated  within  the  plume  and  films  grown  outside  the  plume 
have  been  found  to  be  oxygen  deficient.  The  results  pre¬ 
sented  here  demonstrate  that  high  pressures  of  02  increase 
the  number  density  of  excited  neutrals  atoms  and  excited  ion 
states.  This  implies  that  the  higher-energy  content  of  the 
plume  is  necessary  for  producing  quality  films.  By  placing 
the  substrate  within  the  plume  [typically  4  cm  at  P(02)-  150 
mTorr],  the  substrate  is  bombarded  by  O  I,  Ge  1,  and  Ge  n 
atoms  in  various  electronic  states  as  well  as  by  molecular 
oxygen  in  the  ambient.  Reactions  of  these  species  in  the  gas 
phase  and/or  on  the  substrate  surface  promote  the  proper 
oxygen  incorporation  in  these  films. 

V.  CONCLUSIONS 

The  plume  generated  from  193  nm  irradiation  of  Ge02 
targets  was  examined  by  recording  dispersed  emission  in  the 
UV  to  visible  regions.  The  emission  spectra  were  recorded  as 
a  function  of  laser  fluence,  distance  from  the  front  surface  of 
the  target,  and  molecular  oxygen  pressure.  The  ablation  rate 
of  Ge02  was  also  determined  in  this  work.  Neutral  and  sin¬ 
gly  ionized  atomic  species  were  observed  in  the  plume  and 
the  emission  intensities  of  these  species  depended  on  both 
the  02  pressure  and  distance  from  the  front  surface  of  the 
target.  The  behavior  of  the  emission  intensities  of  these  spe¬ 
cies  were  directly  connected  to  the  excited  atomic  oxygen 
concentrations  which  was  produced  as  a  by-product  of  the 
ablation  process  in  vacuum  and  as  a  product  of  plume-02 
interactions  at  higher  ambient  gas  pressures.  The  change  in 
the  plume’s  concentration  of  excited  species  by  ambient  02 
affected  the  film  growth  by  producing  additional  energetic 
species  that  impacted  the  substrate  which  subsequently  in¬ 
corporated  stoichiometric  proportions  of  oxygen  in  the  film 
structure.  Because  the  excited  species  responsible  for  the 
plume  luminescence  represents  only  a  small  fraction  of  the 
material  removed  from  the  target,  we  are  currently  utilizing 
laser  induced  fluorescence  to  obtain  information  on  the  dy¬ 
namics  and  kinetics  of  the  nonemitting  species  in  the  plasma. 
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Emission  from  several  electronically  excited  states  of  NS  is  observed  when  the  energetic 
molecule  S4N4  is  photolyzed  with  radiation  from  an  excimer  laser.  Photolysis  at  248  nm 

generates  fluorescence  from  the  5  2n1/2,3/2> /f2n  i/2,  G  ,  and  /  2  states  o  . 

NS(fl2n1/23/2)  and  NS(C22  +  )  fluorescence  is  observed  when  the  photolysis  wavelength  is 
changed  to  222  nm.  The  NS (H)  and  NS(C)  spectra  are  postulated  to  arise  from  a 
resonant,  interaction  between  the  KrF  and  KrCl  excimer  photons,  respectively,  and 
vibrationally  hot  ground  state  NS.  LIF  excitation  scans  on  the  NS  X  n1/2.3/2  -  B  “1/2.3/2 
system  confirm  the  production  of  rotationally  and  vibrationally  excited  NS(-T)  up  to 
y'^4.  A  mechanism,  based  on  the  experimental  data  (i.e.,  spectral  composition,  laser  fluence 
studies,  excited  state  time  histories),  calculated  heats  of  formation,  and  Gaussian 
molecular  orbital  calculations,  is  proposed  to  account  for  the  observed  emissions.  For 
photolysis  at  248  nm  it  is  hypothesized  that  a  two  photon  absorption  promotes  the  ground 
singlet  state  of  S4N4  to  an  upper  repulsive  singlet  state,  which  rapidly  dissociates 
(r<30  ns),  producing  an  acyclic  S3N3  fragment  and  vibrationally  excited  monomeric  NS(A). 
The  photofragments  can  interact  further  with  the  excimer  radiation  to  produce  NS (B) 
and  NS(tf),  respectively.  A  similar  mechanism  is  proposed  to  account  for  the  presence  of  the 
NS(5)  and  NS(C)  excited  states  for  the  222  nm  photolysis. 


I.  INTRODUCTION 

The  nitrogen  sulfide  radical  is  an  unstable  radical 
which  cannot  be  stored  as  a  monomeric  liquid  or  solid,  and 
must  be  produced  in  situ  in  the  gas  phase.  NS  is  isoelec- 
tronic  with  NO  and  PO,  and  as  expected,  it  shares  many 
analogous  states  with  these  diatomics.  For  example,  the  NS 
B  2n  and  C22+  electronic  states  have  the  well-known  NO 
B  2I1  and  A  22+  electronic  states  as  their  analogs,  respec¬ 
tively.  In  both  molecules,  these  states  give  rise  to  the  P  and 
y  band  emission,  respectively.  However,  the  NO  f3  bands 
lie  significantly  higher  in  energy  (approximately  40  000 
cm-1)  than  the  NS(52II)  which  lie  near  30  000  cm  . 
The  y  bands  of  both  molecules  are  nearly  resonant  in  en¬ 
ergy,  lying  near  43,000  cm-1. 

Practical  investigations  involving  NS  are  not  nearly  as 
numerous  as  the  isovalent  nitric  oxide,  which  is  known  to 
play  key  roles  in  both  combustion  and  atmospheric  chem¬ 
istry.  However,  the  radical  has  been  detected  in  atmo¬ 
spheric  pressure  flames  by  LIF  methods,1  and  it  has  been 
detected  in  interstellar  clouds  via  radio  frequency  emis¬ 
sion.2,3  The  suggestion  has  also  been  made  that  NS  may  act 
as  the  active  gain  media  in  an  electronic  transition  laser. 
Specifically,  the  NS(52II)  state  of  the  molecule  has  been 
characterized  as  an  attractive  upper  laser  level  primarily 
due  to  its  significant  Franck-Condon  shift  with  respect  to 
the  NS  ground  state  potential  well. 

Fluorescent  emissions  from  NS  are  typically  generated 
by  reacting  active  nitrogen  with  a  sulfur  compound  (e.g., 
H2S,  S2CI2,  SC12)  in  a  discharge  flow  reactor.  Such  emis- 
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sions  have  produced  a  wealth  of  spectroscopic  and  kinetic 
data.  For  example,  the  spectroscopy  of  the  B  2II  state  has 
been  well  characterized  with  detailed  measurements  of 
Franck-Condon  factors,6  r  centroids,7  vibrational  transi¬ 
tion  probabilities,8  and  perturbations.9  Collisional  quench¬ 
ing  and  vibrational  energy  transfer  rates  have  also  been 
measured  for  several  collider  gases.10'"  Less  information  is 
available  on  the  higher  lying  doublet  states  (e.g.,  -  A  2A, 
H  2n,  C  21+ ),  although  studies  utilizing  pulsed  laser  tech¬ 
niques  have  begun  to  detail  the  spectroscopy  of  these  more 
energetic  states.8  In  this  regard,  we  have  recently  reported 
on  the  KrF  laser  photolysis  of  the  inorganic  cage  com¬ 
pound  tetrasulfur  tetranitride.12  S4N4  may  be  regarded  as  a 
labile  polymer  composed  of  four  thiazyl  (NS)  units  (see 
inset  to  Fig.  I).13  The  photolysis  produced  fluorescent 
emission  from  several  of  the  more  energetic  NS  states  in¬ 
cluding  the  tf2nl/2,  G22-,  and  I21  +  .  Strong 
NS(fi  2II)  emission  was  also  observed,  and  measurements 
indicated  the  overall  efficiency  of  conversion  of  absorbed 
photons  by  S4N4  into  NS(52H)  was  2.6% ±0.7%.  The 
facile  photodecomposition  process  observed  at  248  nm  may 
be  attributed  to  the  energetic  nature  of  S4N4 
[Aff/ solid)  =  110  kcal  mol-1],14  the  molecules  large  ab¬ 
sorption  cross  section  at  250  nm  (cr  =  4.8  X  10 
cm2),  and  the  excellent  energy  overlap  of  the  KrF  excimer 
laser  with  the  250  nm  absorption  band. 

This  study  extends  and  supplements  our  initial  inves¬ 
tigation,  and  we  report  on  further  measurements  at  the  248 
nm  photolysis  wavelength,  our  LIF  studies  on  the  produc¬ 
tion  of  ground  state  NS(.V),  and  on  the  fluorescent  emis¬ 
sion  and  products  observed  in  the  photolysis  at  222  nm.  A 
mechanism,  based  on  the  experimental  observations  (spec¬ 
tral  composition,  excited  state  time  histories,  fluence  stud- 
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FIG.  1.  UV-visible  absorption  spectrum  of  2.92  X  10  5  M  S4N4  in 
acetonitrile.  The  inset  figure  shows  the  cage  like  structure  of  S4N4. 


ies,  etc.),  thermochemistry,  and  ab  initio  calculations,  is 
proposed  to  account  for  the  observed  excited  state  and 
ground  state  products. 


II.  EXPERIMENT 

The  experimental  apparatus  consists  of  three  general 
sections  including  the  photolysis  chamber,  the  detection 
ensemble,  and  the  photolysis  and  probe  lasers. 

The  photolysis  chamber  is  constructed  from  a  cylindri¬ 
cal  piece  of  aluminum  tubing  25  cm  in  length  with  an 
inside  diameter  of  3.5  cm.  MgF2  windows  glued  to  the 
tubes  end  allow  for  the  transmission  of  the  laser  radiation. 
LIF  emissions  originating  within  the  cell  were  viewed  off 
axis,  at  the  cells  center  through  MgF2  windows.  The 
S4N4  (typically  <  1  g)  was  contained  in  a  small  10  ml  pyrex 
flask,  which  was  connected  to  the  chamber  through  a  12 
mm  side  port.  Nichrome  wire  wrapped  around  the  neck  of 
the  flask,  and  a  small  heating  mantle  placed  underneath  the 
flask  allowed  for  gentle  heating/increased  sublimation  of 
the  S4N4.  Recrystalization  of  the  tetramer  on  the  cell  sur¬ 
faces  was  prevented  by  wrapping  the  cell  with  heating  tape 
and  heating  the  surface  to  near  70  °C.  Flow  through  the 
cell  was  generated  by  a  mechanical  pump  whose  free  air 
displacement  was  500  1  min-1.  To  insure  a  new  fill  of 
S4N4  on  every  laser  shot  a  flow  of  uhp  helium  buffer  gas 
sufficient  to  flush  the  cell  every  0.1  s  was  introduced 
slightly  upstream  of  the  S4N4  reservoir.  The  temperature  of 
the  He  bath  gas,  measured  at  the  center  of  the  cell,  was 
near  65  °C.  The  total  pressure  in  the  photolysis  zone  was 
measured  with  a  capacitance  manometer  and  was  main¬ 
tained  in  the  range  0.8  to  1.2  Torr.  The  pressures  of  sub¬ 
limed  S4N4  were  typically  <10  mTorr.  Helium  flow  rates 
were  measured  with  a  calibrated  mass  flow  meter.  ' 

The  NS  radicals  produced  in  the  photolysis  were  mon¬ 
itored  by  recording  the  LIF  excitation  spectrum  of  the 
X2U  B2 n  transition  using  a  Nd:YAG  pumped  pulsed 
dye  laser  (PDL)  with  attendant  frequency  doubling.  The 


bandwidth  of  the  frequency  doubled  radiation  was  near 
0.01  nm.  Pulse  energies  <100  /xJ/pulse  were  employed  in 
the  LIF  scans. 

Photolysis  of  S4N4  was  carried  out  with  a  excimer  laser 
operated  on  KrF  (248  nm)  or  KrCl  (222  nm).  The  full- 
width  at  half  maximum  (FWHM)  pulse  width  of  this  laser 
was  1  nm.  The  fluence  of  the  excimer  laser  was  varied  from 
1-85  mJ  cm-2.  The  PDL  could  be  delayed  from  0.5  to  255 
(i s  from  the  photolysis  pulse.  The  actual  minimum  delay 
was  limited  to  —2  /.is  due  to  the  presence  of  prompt  NS(i?) 
fluorescence.  The  photolysis  and  probe  lasers  were  ar¬ 
ranged  to  counter  propagate  through  the  gas  medium,  the 
cross  section  of  the  photolysis  beam  being  much  larger 
than  that  of  the  probe  beam.  The  LIF  was  filtered  and 
detected  by  either  a  1P28  or  a  R 1463-01  photomultiplier 
tube  (PMT),  with  the  latter  tube  specific  for  photon  count¬ 
ing.  The  signals  from  the  PMT  were  sent,  without  ampli¬ 
fication,  straight  to  a  gated  photon  counter  operated  in  the 
boxcar  mode.  The  photon  counter  was  interfaced  to  a 
80286  personal  computer  for  data  collection  and  manipu¬ 
lation.  The  spectral  response  of  the  dye  laser/filter/PMT 
combination  was  measured,  and  used  to  correct  the  LIF 
excitation  spectra.  Spectra  of  the  prompt  emissions  pro¬ 
duced  by  the  photolysis  were  recorded  by  replacing  the 
bandpass  filter  with  a  0.5  m  monochromator.  Alterna¬ 
tively,  prompt  emissions  were  dispersed  with  a  0.25  m 
polychromator  and  detected  with  a  1024  element  silicon 
diode  array.  Calibration  of  the  systems  spectral  sensitivity 
was  performed  with  an  Optronics  Laboratory  deuterium 
lamp  and  standard  lamp. 

The  time  histories  of  the  emissions  were  collected  with 
a  100  MHz  oscilloscope  which  was  interfaced  via  a  GPIB 
to  the  80286  PC.  Storage  and  manipulation  of  the  temporal 
traces  were  performed  with  Asystant/GPIB  software. 

S4N4  was  prepared  by  reacting  sulfur (II) chloride  with 
NH3  in  CC14  solution,  according  to  the  method  given  by 
Villena-Blanco  and  Jolly.15  The  S4N4,  formed  in  the  reac¬ 
tion 


6SC12  +  1 6NH3  -  S4N4  +  2S  +  1 2NH4CI, 


was  isolated  by  a  soxhlet  extraction  with  dry  dioxane.  On 
cooling  large  orange  S4N4  crystals  appeared,  which  were 
filtered  and  dried  in  air.  The  purity  of  the  material  was 
checked  by  FTIR  and  uv-vis.  and  indicated  no  appreciable 
impurities  aside  from  traces  of  sulfur.  A  typical  uv-vis. 
spectrum  of  S4N4  in  acetonitrile  is  given  in  Fig.  1,  and 
shows  one  primary  unresolved  vibronic  band  centered  near 
250  nm  and  a  secondary  much  less  intense  feature  f  —  310 
nm)  to  the  red  of  the  primary.  The  main  uv  band  at  250 
nm  has  been  attributed  to  either  an  n-^rr*  transition  in 
NS16  or  to  a  a->a*  transition  in  the  S-S  interaction.17  The 
250  nm  band  has  a  large  absorption  cross  section  of  4.8 
X  10“ 17  cm2.  Deeper  in  the  uv,  an  additional  strong  ab¬ 
sorption  is  observed  close  to  185  nm.  The  gas12  and  solu¬ 
tion  phase  spectra  are  nearly  identical. 
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FIG.  2.  Emission  spectrum  produced  at  the  248  nm  photolysis  wave¬ 
length.  The  broad  band  near  248  nm  is  due  to  the  KrF  laser  spike. 


III.  RESULTS 
A.  Photolysis  at  248  nm 

On  photolysis  of  ~  10  mTorr  of  S4N4  a  bright  blue 
emission  was  clearly  visible  to  the  eye.  Figure  2  shows  the 
spectrum  of  this  prompt  emission  collected  with  the  0.25  m 
polychromator/diode  array  assembly.  The  emission  ex¬ 
tends  from  the  uv  ( —  200  nm)  into  the  green  ( —550  nm) 
and  comprises  two  distinct  band  systems.  The  broader  and 
more  extensive  double  headed  band  system  was  readily 
identified  as  NS (B  lX\u2Ml  -  X  2n l/2i3/2),  and  shows  the 
characteristic  1/2-3/2  spin  orbit  components  (A0  =  90 
cm-1).  A  second  much  sharper  spectral  progression  is  ev¬ 
ident  in  the  200-330  nm  region.  By  comparison  of  the 
vibrational  term  values  calculated  for  these  bands  and  by 
comparison  of  the  band  frequencies  with  the  known  elec¬ 
tronic  states  of  NS,  a  single  progression  originating  out  of 
»'=0  for  NS(//2n1/2  -  X  2nl/2)  was  identified.  We  have 
recently  completed  a  detailed  analysis  of  this  system,  re¬ 
porting  on  the  Franck-Condon  factors,  r  centroids,  radia¬ 
tive  lifetime,  and  collisional  quenching  constants  of  various 
species  with  the  NS(7/2n)  state.12  As  is  evident  in  the 
spectrum,  the  248  nm  laser  radiation  strongly  overlaps  the 
//1/2(0, 3)  transition  near  248.41  nm.  The  large  Franck- 
Condon  factor  (FCF)  calculated  for  this  band12  (0.1342) 
and  the  sharp  appearance  of  the  progression  suggests  that 
this  state  is  produced  via  a  resonant  interaction  with  the 
laser  radiation. 

A  Stem-Volmer  quenching  study  of  the  NS (77)  state 
indicated  a  collision  free  radiative  lifetime  of  87  ±  14  ns  and 
facile  quenching  by  He,  CF4,  and  N2,  with  the  latter  two 
gases  quenching  near  gas  kinetic  rates.12  Several  bands  to 
the  blue  of  the  KrF  laser  spike  are  not  attributable  to 
NS  (77)  and  have  been  assigned  to  progressions  out  of 
NS(7  22+ )  and  NS (G  21+ ),  although  unequivocal  assign¬ 
ment  of  these  bands  is  made  difficult  by  the  large  number 
of  nested  energy  states  in  the  40  000-45  000  cm  1  range. 

Time  profiles  of  the  NS  ( B )  and  NS  ( H )  emission  were 


FIG.  3.  Temporal  profile  of  the  NS(5— X)  emission. 


recorded  for  various  He  buffer  gas  pressures.  Figure  3 
shows  the  time  profile  of  the  NS (B-X)  intensity  recorded 
for  the  0,4  band  near  394.4  nm  with  no  added  quencher. 
The  profile  reveals  a  rapid  rise  of  -30  ns  followed  by  a 
much  slower  decay  over  several  microseconds.  The  lifetime 
of  the  decay  of  1.1  ±0.2  /zs  compares  favorably  with  the 
known  radiative  lifetime  of  NS(J5,  t/  =  0)  of  1.04 ±0.1  /xs, 
and  we  associate  the  decay  with  radiative  and  collisional 
removal  of  the  excited  state.  The  rise  was  found  to  track 
the  excimer  laser  pulse  (FWHM  pulse  time  =  17  ns),  and 
was  invariant  with  added  quencher.  The  time  history  of  the 
NS (H^X)  fluorescence  was  similar,  and  also  indicated  a 
pressure  invariant,  laser  limited  rise  but  with  a  consider¬ 
ably  faster  decay  rate  as  noted  above,  i.e.,  r l/e  =  87  ns. 
Both  the  NS(77)  and  NS(5)  intensities  were  found  to  vary 
as  the  third  power  of  the  laser  flux. 

From  Fig.  2  it  is  readily  apparent  that  progressions  out 
of  higher  v '  levels  of  the  B  state  are  obscured  due  to  the 
strong  overlap  by  the  H  fluorescence.  However,  since  the  B 
radiative  lifetime  is  considerably  longer  than  the  77,  the  H 
fluorescence  was  easily  discriminated  against  by  setting  a 
photon  counter  gate  delay  equivalent  to  several  77  lifetimes, 
e.g.,  625  ns.  This  gate  delay  also  discriminates  against  any 
overlapping  transitions  from  other  NS  states  in  this  region, 
namely  the  NS(G)  and  NS(7)  whose  lifetimes  are  esti¬ 
mated  to  be  <20  ns.  The  time  delayed  B-X  spectrum 
recorded  with  the  0.5  m  mono/lP28/gated  photon  counter 
is  shown  in  Fig.  4.  The  spectrum,  collected  for  a  gate  width 
of  2  /xs  and  for  an  excimer  flux  of  70  mJ  cm-2,  reveals  a 
highly  excited  nascent  vibrational  distribution  with  popu¬ 
lation  in  levels  up  to  at  least  t’  =10.  It  is  known  that 
NS(2?)  ceases  to  fluoresce  from  vibrational  levels  above 
i/  =  12  because  of  predissociation  to  N(4S)  +  S(  7*). 
However,  population  in  higher  vibrational  levels  could  not 
be  confirmed. 

As  an  aid  in  measuring  the  excited  vibrational  popu¬ 
lations  and  rotational  temperatures  a  spectral  simulation 
code19  was  employed  to  fit  synthetic  spectra  with  the  ex¬ 
perimentally  obtained  spectra,  and  a  relevant  fit  to  the 
NS(2?)  spectrum  is  shown  in  Fig.  4.  The  vibrational  dis¬ 
tribution  yielding  this  fit  is  plotted  in  Fig.  5  which  plots  the 
relative  vibrational  population  against  the  vibrational  term 
energy.  We  note  that  the  intensity  discrepancies  between 
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NS  (B-X)  SPECTRUM  (DELAYED) 
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FIG.  4.  Time  delayed  NS(Z?— X)  spectrum  with  a  superimposed  syn¬ 
thetic  spectrum.  The  peaks  marked  with  an  asterisk  do  not  belong  to  the 
NS(Z?)  system  and  are  probably  due  to  emission  from  higher  lying  states 
of  NS,  eg.,  —NS (//),  NS(G),  and  NS(/). 


the  synthetic  and  experimental  spectrum  in  the  390—450 
nm  region,  i.e.,  notably  the  (0,4),  ( 1,5),  (2,6),  (0,5),  and 
(0,6)  may  be  due  to  photon  counter  pulse  pile  up.  Conse¬ 
quently,  larger  uncertainties  in  the  relative  populations  of 
the  i/=0,  1,  and  2  levels  are  indicated  in  Fig.  5.  From  Fig. 
5  it  is  evident  that  the  photolysis  produces  a  non- 
Boltzmann  vibrational  distribution  with  large  relative  pop¬ 
ulations  in  high  NS(2?)  v  levels.  In  contrast  to  the  vibra¬ 
tional  excitation,  a  fairly  cool  rotational  temperature  of 
600  ±200  K  was  calculated.  Hence,  under  these  photolysis 
conditions  the  nascent  NS(£)  is  generated  vibrationally 
hot  yet  rotationally  cool.  Further,  an  analysis  of  the  rota¬ 
tional  temperature  for  the  H  state12  revealed  a  similarly 
cool  rotational  temperature  of  330  ±  100  K. 


Vibrational  Level 


FIG.  5.  Nascent  NS(f?2ll)  vibrational  population  distribution  produced 
from  the  248  nm  photolysis  of  S4N4.  The  total  pressure  in  the  cell  was  0.8 
Torr. 


Other  photolysis  products  such  as  sulfur  or  metastable 
states  of  nitrogen  were  not  detected.  For  example,  no 
prompt  emission  from  N*  or  S*  was  detected.  Production 
of  the  first  excited  state  of  nitrogen,  N 2(A  3X  + ),  was  tested 
for  by  adding  a  small  amount  of  NO  to  the  photolysis 
chamber.  The  presence  of  N2(T)  would  typically  be  sig¬ 
naled  by  strong  NO  y  band  emission  due  to  the  rapid 
energy  transfer  reaction:  N2M)  +  NO(X)  —  N2(X) 
+  NO(/f).20  No  y  band  emission  was  detected.  From 
NO(/4)  detection  limits  determined  in  previous  flow  tube 
experiments21  we  estimate  the  N2(,4)  density  to  be  <  1 
X  106  mol  cm“3.  The  production  of  S2(X)  was  investi¬ 
gated,  as  described  below,  by  LIF  excitation  scans  on  the 
S2(X  3X~  -  B  3X“ )  system. 

B.  Photolysis  at  222  nm 

Photolysis  of  S4N4  at  222  nm  probes  vibronic  struc¬ 
ture,  lying  on  the  wing,  to  the  blue  of  the  primary  absorp¬ 
tion  centered  near  250  nm.  The  absorption  cross  section  at 
222  nm  is  equal  to  8.14  X  10~18  cm2,  and  while  consider¬ 
ably  reduced  from  the  248  nm  value  is,  nonetheless,  sub¬ 
stantial  and  provides  a  further  probe  of  the  photodissocia¬ 
tion  dynamics  of  the  molecule. 

Photolysis  of  10-20  mTorr  of  the  tetramer  at  the  KrCl 
wavelength  again  produced  a  visible  blue  fluorescence,  al¬ 
beit  much  weaker  in  intensity  than  that  observed  in  the  248 
nm  case.  Wavelength  dispersion  of  the  emission  produced 
a  spectrum,  similar  to  that  shown  in  Fig.  3,  with  two  dis¬ 
tinct  band  systems.  NS(2?— X)  emission  was  again  ob¬ 
served  in  the  270-550  nm  region  with  fluorescence  origi¬ 
nating  from  v  levels  as  high  as  eight.  However,  the  NS(£) 
emission  intensity  generated  via  222  nm  photolysis  was 
<0.1  of  the  NS(i?)  intensity  observed  during  the  248  nm 
photolysis,  for  similar  conditions  of  fluence.  More  strik¬ 
ingly,  NS(ff)  emission  was  no  longer  observed  in  the  short 
wavelength  region  being  replaced  with  the  double  banded 
spectrum  shown  in  Fig.  6(a). 

This  spectrum  is  attributable  to  the  NS(C2X  + 
X  2n1/2(3/2)  transition  and  shows  only  a  single  progres¬ 
sion  originating  out  of  v'  —  2.  Each  band  in  the  spectrum 
shows  two  components  separated  by  221.5  cm"1  equiva¬ 
lent  to  the  spin  orbit  energy  splitting  of  the  NS  ground 
state.22  The  short  and  long  wavelength  components  are  due 
to  the  1/2  and  3/2  transitions,  respectively. 

Temporal  histories  recorded  for  the  (2,4)  and  (2,5) 
transitions  with  no  added  diluent  gas  showed  very  spiked 
Gaussian  intensity  profiles  with  a  FWHM  of  17  ns.  The 
profiles  were  unaffected  by  addition  of  Ar  up  to  9  Torr. 
From  these  time  histories  the  NS(C2X  +  ,  v'  =  2)  fluores¬ 
cent  lifetime  is  estimated  to  be  <10  ns.  This  is  consistent 
with  a  6.5  ns  radiative  lifetime  determined  by  Silvers  and 
Chiu  for  the  NS(C)  t/  =  0  level  via  a  Hanle  effect  measure¬ 
ment.23  In  a  more  detailed  study  of  the  NS(C)  state  Chiu 
and  Silvers24  observed  absorption,  emission,  and  fluores¬ 
cence  from  the  t/  =  0,  1,  and  2  levels.  They  noted  a  weak¬ 
ening  of  the  bands  originating  from  u'  =  1  and  2  relative  to 
that  from  the  </=0  level.  They  ascribed  this  diminution  in 
intensity  to  a  weak  predissociation  in  the  v'  =  1  and  2  levels, 
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FIG.  6.  (a)  Emission  spectrum  collected  in  the  200-300  nm  region  for 
photolysis  at  222  nm.  The  broad  band  near  222  nm  is  the  KrCl  laser 
spike,  (b)  Synthetic  spectrum  of  the  NS(C22V  =  2)  progression  show¬ 
ing  the  relative  magnitude  of  the  (2,1)  band. 


an  idea  supported  by  configuration  interaction  calculations 
on  the  C22+  adiabatic  potential  well.25  Hence,  the  short 
radiative  lifetime  we  measure  for  the  y'=2  level  may  be 
due  to  predissociation. 

It  is  of  interest  to  compare  the  photolytically  induced 
spectrum  in  Fig.  7(a)  with  Chiu  and  Silver’s  fluorescence 
spectrum  (their  Fig.  3)  obtained  by  exciting  the  NS(X) 
produced  in  an  active  nitrogen  discharge  with  a  deuterium 
lamp.  They  observe  strong  emission  out  of  u'=0  and  much 
weaker  v'=\  and  2  progressions.  Moreover,  their  (2,u") 
band  intensities  are  barely  discemable  from  the  baseline 
noise,  and  only  the  (2,1),  (2,3),  and  (2,4)  are  observable. 
In  marked  contrast,  the  photolytically  induced  spectrum 
shows  no  emission  from  the  i/=0  and  1  levels  and  intense 
emission  only  from  o' =  2,  with  this  progression  extending 
to  at  least  o'  =  8.  This  behavior  suggests  interaction  of  the 
222  nm  excimer  radiation  with  an  energetically  resonant 
state  of  NS.  From  the  C  spectrum  it  is  clear  that  the 
(2,1)  1/2  and  (2,1 ) J/2  bands  at  222.46  and  223.56  nm  are 
resonant  with  the  222  nm  laser  pulse  (FWHM=1  nm), 
and  that  excimer  laser  pumping  of  these  transitions  may 
give  rise  to  the  observed  spectrum,  provided  the  photode¬ 
composition  creates  a  nascent  ground  state  population 
which  is  vibrationally  excited,  i.e.,  N S ( T ,  o"  =  l)  must  be 
populated.  As  noted  above,  in  the  248  nm  case  an  analo- 
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FIG.  7.  The  apparent  dependence  of  the  electronic  transition  moment  on 
the  r-centroid  for  NS(C22\  v'  =  2)  as  given  by  a  plot  of 
(/V  P-) 1/2  vs  r ■  The  solid  iine  rePresents  a  fourth-order  polynomial  fit 

to  the  data.  D — estimated  value  for  the  (2,2)  band  based  on  pii 
~  tooo — 2  +  — Jeffries  et  al.  data  for  NS(C22  +  ,  u’  =  0);  this 

data  was  not  included  in  the  polynomial  fit, 

gous  LIF  process  is  hypothesized  to  account  for  NS (H) 
fluorescence,  i.e.,  the  248  nm  radiation  pumps  the 
NS(JT,  um  =  3)  -  NS (//j/2,  V  =  0)  transition  at  248.4 
nm,  thereby  requiring  nascent  population  in  NS(X  v" 
=  3).  LIF  excitation  scans  on  the  NS  X-*B  transition, 
described  below,  confirm  the  production  of  vibrationally 
excited  ground  state  populated  up  to  at  least  v"  =4  in  both 
the  222  and  248  nm  photolysis. 

A  number  of  experiments  were  performed  in  which  the 
relative  intensity  of  the  NS  C-X  emission  was  recorded  as 
the  incident  fluence  of  the  222  nm  pulse  was  varied.  The 
data  indicate  a  quadratic  dependence  of  the  NS  (C-X)  in¬ 
tensity  on  excimer  fluence. 

C.  Determination  of  Franck-Condon  factors,  R 
centroids,  and  vibrational  transition  probabilities  for 
the  NS (C  2S+)  state 

Due  to  the  enhancement  of  the  NS(C,  v”  =  2)  emission 
spectrum  relative  to  others  reported  to  date,24  a  more  ex¬ 
tensive  spectral  analysis  of  the  emission  was  undertaken, 
with  calculations  of  FCFs,  r  centroids,  vibrational  transi¬ 
tion  probabilities  and  transition  moments.  An  RKR  com¬ 
puter  code  was  employed  in  calculating  the  potential  func¬ 
tion  U(r)  and  classical  turning  points  for  both  NS(C)  and 
NS  (JO  using  existing  Dunham  coefficients  for  each  state.2” 
FCF’s  and  r  centroids  were  thus  obtained  for  the  first 
three,  and  only  observed,  vibrational  levels  of  the  C  state, 
and  are  presented  in  Table  I.  These  values  compare  favor¬ 
ably  with  past  values  calculated  using  Morse  wave  func¬ 
tions.8,24 

Vibrational  transition  probabilities  Ov.,,0  were  calcu¬ 
lated  from  the  relation 

where  A  is  the  Einstein  coefficient  for  spontaneous  emis¬ 
sion  and  is  proportional  to  the  intensity  of  fluorescence  in 
a  given  vibrational  band  at  wavelength  Relative 
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TABLE  I.  RKR  Franck-Condon  factors  and  r  centroids  for 


NS(C22+). 


Transition 

(q'n+qM)/  r 

(r  -)/2 

0,0 

612 

1.474 

0,1 

280 

1.425 

0,2 

80 

1.379 

0,3 

18 

1.335 

0,4 

3.8 

1.296 

0,5 

0.64 

1.257 

0,6 

0.11 

1.215 

1,0 

311 

1.530 

1,1 

161 

1.485 

1,2 

301 

1.432 

1,3 

158 

1.388 

1,4 

50 

1.342 

1,5 

14 

1.303 

1,6 

3.2 

1.270 

1,7 

0.65 

1.220 

2,0 

68 

1.590 

2,1 

382 

1.539 

2,2 

18 

1.505 

2,3 

202 

1.435 

2,4 

200 

1.400 

2,5 

88 

1.350 

2,6 

31 

1.304 

2,7 

9 

1.287 

2,8 

2 

1.242 

“Values  multiplied  by  1000. 


Av>v -  were  obtained  from  area  integration  of  the  discreet 
v\v"  spectral  bands  taken  from  the  0.25  m  polychromator/ 
diode  array  response  corrected  spectrum  shown  in  Fig. 
6(a)  [excluding  the  (2,1)  band].  Using  the  known  A,.-/'2 
Eq.  (1)  was  used  to  calculate  relative  p^y  for  the  i/  =  2 
progression.  The  results,  presented  in  Table  II,  are  com¬ 
pared  with  the  RKR  computed  Franck-Condon  factors 
(qvy)f  and  indicate  a  gradual  increase  in  pv>y  over  qv-y 
with  increasing  v*.  This  implies  a  variation  of  the  elec¬ 
tronic  transition  moment  [He(JV-0]  with  internuclear  sep¬ 
aration.  The  r  centroid  relation,  6 

Po\i>*  —  )  1  Qv\v*> 


TABLE  II.  Comparison  of  NS(C2£+.  u'  =  2)  vibrational  transition 
probabilities  with  the  RKR  Computed  Franck-Condon  factors. 

v'  —  2,  v*  = 

iq'^  +  q3'2)/  2‘ 

0 

68 

109 

1 

382 

(254)b 

2 

18 

16 

3 

202 

202c 

4 

200 

194 

5 

88 

107 

6 

31 

50 

7 

9 

36 

8 

2 

30 

*p  and  q  values  multiplied  by  103. 

hThe  pi ,  value  was  estimated  from:  p2l  cr  1000 — X 

cThe  pi,)  value  is  normalized  to  the  }  value. 

was  used  to  quantify  this  variation.  Figure  7  plots  the  elec¬ 
tronic  moment  as  (pv'y/qv>y) 1/2  vs  r  centroid.  According 
to  the  RKR  calculations,  the  r  centroid  varies  from  1.24  A 
for  the  (2,8)  band  to  1.59  A  for  the  (2,0).  Over  this  in¬ 
terval  a  smooth  decrease  in  Re{rv>y)  is  readily  apparent.  A 
similar  variation  was  noted  by  Jeffries  et  a/.8  in  their  LIF 
spectroscopic  study  of  the  NS(C,  u'  =  0)  state,  and  we  have 
included  their  data  for  comparison  in  Fig.  7.  A  fourth 
order  polynomial  fit  to  the  averaged  1/2  plus  3/2  data 
yielded  the  following  expression: 

Re(rv.y )  =  const  [  1 —2.76  r+2. 86  r 2— 1.31  r3 

+0.226  r*  (3) 

and  gives  the  variation  of  the  transition  moment  for  the 
region  sampled:  1.2  A  <  rv<y  <  1.6  A. 

Introducing  the  functional  form  for  Re(rv>y )  given  by 
Eq.  (3),  and  using  the  RKR  calculated  qv>y  (see  Table  I) 
in  the  spectral  simulation  code  a  best  fit  NS(C)  spectrum 
was  computed  and  is  graphically  overlaid  on  the  experi¬ 
mental  spectrum.  The  simulation  correctly  predicts  the 
(2,2)  vibrational  node  and  fits  the  higher  v*  bands  cor¬ 
rectly  in  accord  with  the  increasing  Re(rv<y)  in  this  region. 
Based  on  the  width  of  the  vibronic  bands  a  rotational  tem¬ 
perature  of  2500  ±  500  K  was  determined.  Since  the  (2,1) 
band  should  have  considerable  intensity  [i.e.,  q{2, n 
=  0.382),  and  since  this  band  is  spectrally  obscured  by  the 
222  nm  laser  light  a  synthetic  spectrum  of  the  entire  i/  =  2 
progression  is  presented  in  Fig.  6(b).  This  spectrum  indi¬ 
cates  the  relative  magnitude  of  the  (2,1)  band. 

D.  NS(X  2n  -  B  2n)  LIF  excitation  scans  for  248  nm 
photolysis 

The  operation  of  photolytic  channels  which  produce 
rotationally  and  vibrationally  excited  NS(Z)  was  probed 
by  searching  for  LIF  from  NS  fragments.  Excitation  spec¬ 
tra  were  recorded  by  monitoring  NS(2?— X)  emission 
through  bandpass  filters.  Two  separate  wavelength  regions 
were  scanned:  300-310  nm  and  360-372  nm  using  Rh  640 
and  LDS  750  dyes,  respectively.  These  experiments  clearly 
showed  that  rotationally  and  vibrationally  excited  NS(A0 
is  produced  by  the  248  nm  photolysis  as  indicated  by  the 
LIF  spectra  shown  in  Figs.  8  and  13.  All  of  the  spectra 
were  corrected  for  the  wavelength  dependence  of  the 
PDL/filter/PMT  combination.  As  mentioned  previously, 
the  minimum  excimer  pulse  to  PDL  pulse  delay  time  was 
limited  to  2  ps  due  to  the  prompt  NS (B)  fluorescence.  No 
NS  LIF  was  observable  with  the  PDL  probe  pulse  alone. 

As  a  potential  photolytic  product,  the  presence  of  di¬ 
atomic  sulfur  was  tested  by  scanning  for  LIF  on  the 
-  £32“)  transition.  Two  wavelength  regions 
were  examined:  300-310  nm  and  360-370  nm.  Since  the 
lifetime  of  S 2(B)  (37  ns)  is  considerably  shorter  than 
NS (B)  the  NS (B)  fluorescence  was  discriminated  against 
by  setting  a  small  gate  delay  with  a  gate  width  equivalent 
to  several  S 2(B)  lifetimes,  i.e.,  100  ns.  Further,  several  of 
these  experiments  were  carried  out  with  high  PDL  fluences 
(laser  induced  pumping  rates  >  1  X  1010  s"1)  in  order  to 
saturate  the  S2  transitions.  Under  these  experimental  con- 
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NS  (X->B)  Excitation  Scan 


3010  3020  3030  3040 

WaVftUllflth  (A) 

FIG.  8.  Experimental  (Expt)  and  simulated  (Sim)  LIF  excitation  spectra 
of  the  (4,0)  transition  of  the  B  2Hi/2.3/2  system  of  NS  pro- 

duced  via  the  photolysis  of  S4N4  at  248  ran.  1  and  2  designate  the  1/2  and 
3/2  spin-orbit  components,  respectively.  The  total  pressure  in  the  pho¬ 
tolysis  cell  was  0.8  Tort,  and  the  probe  dye  laser  delay  time  was  70  /is. 


ditions  no  S 2(B  -  X)  fluorescence  was  detected.  Based  on 
a  two  level  model  for  near  saturated  fluorescence  in  di¬ 
atomic  molecules27  we  calculate  a  S2(Ar32g  )  detection 
limit,  for  our  experimental  conditions,  of  1  X  10 
mol  cm-3. 

E.  Internal  energy  distribution  of  the  NS(X) 
fragments 

Since  the  linewidth  of  the  probe  laser  (0.01  nm)  was 
not  sufficient  to  fully  resolve  the  individual  rotational  lines 
of  the  NS  LIF  spectra,  the  spectral  simulation  code  was 
modified  and  was  used  as  an  aid  in  determining  vibrational 
and  rotational  population  distributions.  The  simulation 
employed  a  Gaussian  line  shape  function,  a  laser  line  width 
parameter  and  the  following  equation  for  rotational  line 
intensities: 

n-PvQ~ W -BmN"  (AT  + 1  )/kT},  ^ 

where  k  is  a  proportionality  constant,  q  is  a  Franck- 
Condon  factor,  vuV  is  the  transition  frequency,  $>.*•  is 
the  Honl-London  factor  appropriate  to  a  AA=0  transi¬ 
tion,  Pv-  is  the  population  of  the  u"  vibrational  level,  Q  is 

J.  Chem.  Phys.,  Vol. 
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Rotational  Energy  {cm'1) 


FIG.  9.  Quantum  state  distribution  of  NS(t»w— 0)  as  measured  from 
(4,0)  excitation  spectra  collected  at  photolysis  to  probe  pulse  delay  in¬ 
ternals  of  2  and  70  /zs,  respectively.  The  latter  spectrum  is  shown  in  Fig. 
8. 


the  partition  function,  BH  is  the  rotational  constant,  and  T 
is  the  rotational  temperature.  Equation  (4)  is  applicable  to 
a  Boltzmann  rotational  distribution.  An  excitation  scan  for 
the  (4,0)  band  of  the  B-X  system  is  presented  in  Fig.  8. 
The  appearance  is  typical  for  an  unperturbed  band  of  this 
parallel  transition,  and  should  show  an  R  branch  head,  a 
red  degraded  P  branch,  and  no  Q  branch.  However,  due  to 
the  limited  laser  resolution  the  P  and  R  branches  remain 
unresolved.  The  spectrum  collected  after  a  70  /is  delay 
appears  thermalized  (A^ax  —  ^LS),  and  is,  indeed,  well  fit 
by  a  spectral  simulation  utilizing  a  340  K  rotational  tem¬ 
perature.  A  excitation  spectrum  recorded  for  a  shorter 
photolysis  to  probe  pulse  delay  interval  (2  /xs)  showed  a 
much  flatter  and  broader  rotational  distribution  with  a 
much  less  distinct  R  branch  head  and  population  in  higher 
jV"  (tv"'  =58.5),  behavior  consistent  with  a  more  highly 
excitedTotational  population  in  the  ground  state.  A  syn¬ 
thetic  fit  to  this  spectrum  indicated  a  rotational  tempera¬ 
ture  of  800  ±200  K.  These  “best  fit”  iys  are  in  close 
agreement  with  the  Tr  values  calculated  directly  from  the 
observed  (response  corrected)  spectra.  This  is  indicated  in 
Fig.  9,  which  plots  log (//5y\,y)  vs  N  (N  +1)  for 
the  (4,0)  excitation  scans  described  above.  The  plots  are 
quite  linear,  indicative  of  a  Boltzmann  distribution  in  this 
region.  From  the  slopes  of  the  lines  rotational  temperatures 
of  800  and  350  K  were  determined  for  the  2  and  70  (is 
delay  times,  respectively. 

Rotational  relaxation  was  evident  as  the  pump  to 
probe  delay  time  was  increased.  This  is  shown  in  Fig.  10, 
which  plots  the  rotational  temperature  of  the  (4,0)  x/2  band 
vs  number  of  collisions  with  the  helium  bath  gas  for  low 
fluence  (F<10  mJcm^2)  photolysis.  Evidently,  NS(Ar,  v 
=0)  is  formed  with  a  near  nascent  rotational  temperature 
of  800  ±200  K,  and  a  few  hundred  collisions  are  required 
to  thermalize  the  distribution. 

The  variation  in  the  (4,0)  1/2  LIF  intensity  as  a  func¬ 
tion  of  laser  fluence  was  recorded  for  conditions  of  both 
high  and  low  laser  fluence.  For  these  measurements  the 
photolysis  to  probe  pulse  delay  was  set  to  20  ^.s,  sufficient 
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FIG.  10.  Collisional  cooling  of  NSCX  v”  =0,  jV")  for  low  fluence  pho¬ 
tolysis  conditions.  H — 248  nm  photolysis,  A — 222  nm  photolysis. 


FIG.  12.  Variation  of  the  NSCX,  u*=0)  rotational  temperature  with  KrF 
laser  fluence.  The  solid  line  represents  an  exponential  fit,  of  the  form 
y  —  j0exp(fl[X),  to  the  data. 


time  to  allow  for  the  relaxation  of  the  rotational  manifold. 
In  the  low  fluence  regime  (0.8  —  25  mJ/cm2)  the  rotational 
distributions  were  indeed  thermalized.  Figure  1 1  shows  a 
log-log  plot  of  the  integrated  fluorescent  intensity  vs  exci- 
mer  laser  fluence  for  the  low  fluence  photolysis.  A  least 
squares  fit  to  the  data  yielded  a  slope  of  2.0 ±0.2.  This 
quadratic  dependence  indicates  that  2  KrF  photons  are 
required  to  generate  an  NS(X)  photofragment.  For  the 
high  fluence  measurements  the  rotational  distributions  did 
not  appear  relaxed,  and  a  substantial  variation  in  the  rota¬ 
tional  quantum  state  distribution  was  evident  as  the  laser 
flux  was  varied.  This  effect  is  quantified  in  Fig.  12,  which 
graphs  the  rotational  temperature  of  the  (4,0)  l/2  band  vs 
laser  fluence.  It  is  apparent  that  under  high  fluence  condi¬ 
tions  the  rotational  distribution  is  much  more  excited  than 
would  be  predicted  for  a  quenched  distribution.  A  nearly 
exponential  decrease  in  rotational  temperature  with  fluence 
is  seen,  and  as  the  fluence  is  reduced  the  temperature  ap¬ 
proaches  the  thermal  limit  (near  340  K).  This  behavior  is 
suggestive  of  a  fluence  dependent  rotational  excitation  pro¬ 
cess  in  the  gas  medium. 


In  Fluence  (mJ/cm2) 

FIG.  II.  A  plot  of  the  log  of  the  integrated  emission  intensity  of  the 
NS(4,0)1/2  band  vs  log  of  the  KrF  laser  fluence.  The  slope  of  2.0±0.2 
indicates  two  photons  are  required  to  generate  the  NS  ground  state. 


The  variation  in  the  LIF  intensity  of  the  (0,2)  1/2  band 
was  also  recorded  as  a  function  of  laser  fluence.  In  this 
experiment  the  delay  separating  photolysis  to  probe  pulse 
was  set  to  70  /zs,  and  the  KrF  fluence  was  varied  in  the 
range  3  to  17  mJ  cm“2.  The  results,  for  this  higher  vibra¬ 
tional  level,  also  indicated  a  quadratic  dependence  of  in¬ 
tensity  with  fluence. 

Population  in  vibrational  levels  above  v"  =0  was 
probed  by  scanning  the  360-372  nm  region.  Figure  13  dis¬ 
plays  a  typical  LIF  excitation  scan  of  this  region,  collected 
for  a  delay  interval  of  70  jzs  and  at  a  total  pressure  of  1.2 
Torn  Under  these  experimental  conditions  the  rotational 
distribution  is  equilibrated  (i.e.,  rrot  =  340  K)  while  the 
vibrational  distribution  remains  unperturbed.  Hence,  the 
spectral  simulation  designed  to  determine  relative  Pv -  was 
simplified  as  only  a  single  rotational  temperature  in  Eq. 
(4)  needed  to  be  specified. 

The  rate  constant  for  quenching  vibrationally  excited 
NS  (JO  by  He  has  not  been  reported.  However,  we  note 
that  excitation  spectra  collected  under  low  collision  condi¬ 
tions  ( ~30  collisions)  and  under  high  collision  conditions 
( —  1300  collisions)  showed  no  significant  differences  in  the 
vibrational  distribution.  This  indicates  the  NS(Ar)-He 
VET  is  not  a  facile  process. 

Population  in  the  u"  =  2,  3,  and  4  levels  is  evident  in 
the  spectrum  of  Fig.  13(a).  Population  in  u"=4  is  indi¬ 
cated  by  the  presence  of  the  weak  (3,4)  bandhead  near  361 
nm.  The  spectral  code  was  used  in  determining  the  relative 
vibrational  populations  in  the  2,  3,  and  4  levels  and  an 
optimized  fit  is  shown  in  Fig.  13(b).  The  fractional  vibra¬ 
tional  populations  obtained  in  this  manner  are  0.89,  0.07, 
and  0.03  for  v”=  2,  u"  =  3,  and  v"=4f  respectively.  This 
distribution  was  found  to  be  independent  of  laser  fluence. 
The  fit  was  such  that  changing  Pv*= 4  from  0.03  to  0  re¬ 
moved  the  (3,4)  bandhead,  in  clear  disagreement  with  the 
experimental  data.  Population  in  higher  v"  levels  was  not 
evident.  The  region  of  the  (5,6)1/2  band  in  the  vicinity  of 
371.0  nm  was  scanned.  The  FCF  for  this  band  is  equivalent 
to  the  (0,2)  transition,  with  strong  LIF  signals  from  the 
latter  being  observed.  However,  no  LIF  signals  were  de- 
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FIG.  13.  Experimental  (a)  and  synthetic  (b)  NS  {X- B)  LIF  excitation 
spectra  showing  population  in  v"  =  2,  v"  —  3 ,  and  n#=4  vibration  levels. 

tected  in  this  region  even  with  substantially  higher  pres¬ 
sures  of  S4N4  in  the  photolysis  cell.  From  determinations  of 
signal  to  noise  ratios,  i.e.,  -2S/N  in  the  LIF  scans  we 
estimate/V  =b  <  0.01 . 

IV.  DISCUSSION 

The  photolytic  process  at  248  nm  leading  to  the  for¬ 
mation  ofNS(fl)  and  NS  (H)  is  complicated  by  the  variety 
of  (SN)X  intermediates  which  may  be  formed,  and  by  the 
energetics  of  these  species.  However,  in  view  of  the  data, 
presented  above,  it  appears  that  the  NS(2?)  and  NS (H)  are 
formed  from  different  precursor  species.  For  example,  the 
emission  spectrum  (Fig.  4)  indicates  the  nascent  NS (B)  is 
produced  vibrationally  hot,  in  a  non-Boltzmann  distribu¬ 
tion  with  v'  levels  as  high  as  ten  populated.  The  NS {B)  is 
also  formed  with  some  degree  of  rotational  excitation  as 
rrot  =  600  ±  200  K.  In  contrast,  the  NS  (tf)  appears  rota- 
tionally  cold,  and  the  collected  emission  is  reminiscent  of  a 
fluorescence  scan  with  only  a  single  vibrational  level  ( v 
=0)  emitting.  Additional  considerations,  as  previously 
noted,  such  as  the  energy  overlap  of  the  H{ 0,3)  transition 
with  the  KrF  frequency,  the  population  detected  in  ground 
state  vibrational  levels,  and  the  disappearance  of  the  H 
fluorescence  when  the  excimer  wavelength  is  changed 
clearly  indicates  the  H  state  is  formed  by  laser  excitation  of 
vibrationally  excited  ground  state  NS  fragments.  Since  the 


KrF  wavelength  is  above  the  dissociation  limit  for  NS (B), 
the  B  fluorescence  cannot  arise  from  an  analogous  KrF 
excitation  of  NSW  molecules.  Rather,  some  higher  order 
(SN)X  homolog  must  be  the  precursor  to  NS (B).  We  fur¬ 
ther  note  that  the  temporal  history  of  the  B  emission  indi¬ 
cates  that  the  NS (B)  is  not  generated  via  a  collisional 
process.  Rather,  its  production  appears  photolytically  in¬ 
duced.  The  fluence  data  indicates  that  NS (B),  as  well  as 
the  more  energetic  NS(/f),  is  formed  in  a  three  photon 
process.  Aside  from  the  prompt  NS(5)  and  NS(H)  emis¬ 
sions,  the  excitation  scans  reveal  that  rotationally  and  vi¬ 
brationally  excited  NSW  is  also  a  product  of  the  photol¬ 
ysis,  and  is  produced  in  a  two  photon  event.  In  light  of 
these  considerations,  the  following  general  mechanism  is 
proposed  to  qualitatively  account  for  the  NS(5),  NS( //), 
and  NS(A»  states  produced  in  the  KrF  photolysis  of 


S4N4: 

C-S4N4  +  2h  v  -  a-S3N3  +  NS  ( X,v ) , 

Atf=-132  kcal  mol-1, 

(Rl) 

a-S3N3  +  Av->NS(5,y)  +c-S2N2, 

A//=  —21  kcal  mol-1, 

(R2) 

NSWu"  =3)  +Av— NS(//,u'=0), 

(R3) 

where  reaction  (R3)  shows  the  resonant  excitation  which 
produces  the  NS(//)  state,  and  where  c  and  a  represent  the 
cyclic  and  acyclic  isomers,  respectively.  The  heats  of  reac¬ 
tion  have  been  calculated  from  A Hf(SxNx)  values  previ¬ 
ously  determined  via  a  MOPAC  calculation. ‘8 

Reaction  (Rl)  depicts  a  multiphoton  excitation  of  cy¬ 
clic  S4N4  to  an  excited  singlet  state  which  subsequently 
dissociates  yielding  vibrationally  excited  NSW  fragments 
and  an  acyclic  S3N3  fragment.  The  large  value  of  the  exo- 
thermicity  for  this  process  is  not  surprising,  considering 
the  2  uv  photon  excitation  contributes  230  kcal  mol  at 
the  248  nm  wavelength.  The  heat  of  reaction  is  sufficient  to 
provide  for  the  vibrational  excitation  of  the  NS(A')  frag¬ 
ment  (u"=4  requires  13.5  kcal),  and  for  a  large  remnant 
energy  balance  of  over  100  kcal  mol-1.  For  comparison,  a 
one  photon  excitation  would  be  much  less  exothermic 
(-17  kcal  mol-1)  and  would  be  nearly  thermal  neutral 
when  the  vibrational  energy  in  the  NS  fragment  is  ac¬ 
counted  for.  That  two  photons  are  required  to  reach  the 
dissociative  state  is  also  supported  by  the  results  of  MO 
calculations  on  c-S4N4,  as  described  below. 

If  the  partitioning  of  excess  energy  between  the  NSW 
and  S3N3  photofragments  in  reaction  one  is  statistical  in 
nature,  the  more  massive  S3N3  fragment,  with  12  vibra¬ 
tional  modes,  would  consume  the  majority  of  the  vibronic 
energy.  Further,  provided  the  excited  S4N4  electronic  state, 
from  which  dissociation  occurs,  has  approximately  the 
same  S-N  geometry  as  the  S4N4  ground  state,  the  energy 
available  to  the  NS  fragment  would  not  be  expected  to 
exceed  that  available  from  simple  geometric  relaxation  ot 
the  NS  bond  length  in  c-S4N4.  The  energy  available  in 
contracting  the  N-S  bond  length  in  C-S4N4  (1.63  A)  to 
its  equilibrium  value  in  monomeric  NS  (1.495  A)‘“  is 
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readily  obtained  from  the  RKR  potential  well  of  the  NS 
ground  state.  We  calculate  that  8.6  kcal  mol -i  of  vibra¬ 
tional  energy  can  be  imparted  to  the  released  NS  (AT  frag¬ 
ment.  This  represents  enough  energy  to  populate  NS  (AT 
up  to  v”=3.  This  result  agrees  reasonably  well  with  the 
experimental  measurements  which  showed  small  relative 
populations  of  0.07  and  0.03  in  the  v"  =  3  and  v"  =  4  levels, 
respectively,  and  no  measurable  population  in  higher  lev¬ 
els.  Thus,  it  appears  that  the  excited  state  is  sufficiently 
long  lived  (on  the  order  of  a  few  vibrational  periods)  to 
allow  for,  at  least,  partial  vibrational  energy  randomization 
in  the  transitional  S4N4  structure. 

Evidently,  from  the  rotational  excitation  observed  in 
the  NS  (A",  y=0)  fragment,  a  relatively  small  fraction  of 
the  excess  energy  available  in  the  primary  photodissocia- 
tion  ends  up  in  NS  rotation.  Rotational  energies  much 
higher  than  1000  cm"1  (2.9  kcal)  were  not  detected.  In¬ 
tuitively,  the  presence  of  some  rotational  excitation  in  the 
NS  photofragments  indicates  that  the  dissociation  is  not 
linear,  and  that  there  is  some  excitation  of  the  bending 
modes  in  the  transition. 

The  photofragments  produced  in  reaction  one  may  un¬ 
dergo  secondary  photolytic  processes,  as  indicated  by  re¬ 
actions  (R2)  and  (R3).  Reaction  (R2)  shows  the  NS (B) 
being  formed  from  the  single  photon  excitation/ 
dissociation  of  an  a-S3N3  intermediate.  The  observation  of 
NS(2?,  umax  =  10)  places  an  energy  constraint  on  reaction 
(R2),  and  we  are  required  to  choose  the  cyclic  form  of 
S2N2  over  the  acyclic  form,  since  the  latter  would  lead  to 
j-f 23  kcalmol-1  endothermic  process.  The  reaction  exo- 
thermicity  thus  calculated  (—21  kcal  mol-1)  is,  consider¬ 
ing  the  error  of  the  MOPAC28  calculation,  just  sufficient  to 
populate  NS  up  to  its  observed  vibrational  maximum. 
Thus,  it  would  appear  that  the  reaction  exothermicity  is 
completely  consumed  in  NS  (2?)  vibration.  Since  it  is  rea¬ 
sonable  to  assume  that  some  energy  is  partitioned  to 
NS (2?)  translation  and  S2N2  vibration-translation,  the  A H 
for  reaction  (R2)  is  likely  to  be  <—21  kcalmol.  This 
would  require  a  larger  energy  content,  than  determined  by 
the  MOPAC  calculation,  for  the  a- S3N3  precursor.  As  ex¬ 
plained  above,  the  <z-S3N3  molecules  are  likely  to  be 
formed  with  considerable  internal  excitation. 

Other  S4N4  photodecomposition  mechanisms  may  cer¬ 
tainly  be  written.  However,  the  majority  of  these  can  be 
ruled  out  based  on  considerations  of  energetics  and  the 
experimental  data.  For  example,  a  mechanism  involving 
the  nonlinear  absorption  of  three  photons  by  the  tetramer 
to  produce  NS(B)  and  NS (7/), 

c-S4N4  +  3  A v248  -  NS  ( B)  +  NS  ( H)  +  S2N2,  ( R4 ) 

is  endothermic  by  37  kcal  mol-1,  and  does  not  account  for 
the  specific  NS (//,  */  =  0)  fluorescence  or  the  two  photons 
required  to  produce  the  ground  state. 

The  data  for  the  KrCl  photolysis  is  qualitatively  simi¬ 
lar  to  that  obtained  in  the  KrF  photolysis.  For  example, 
the  222  nm  emission  spectrum  indicates  formation  of  vi- 
brationally  hot  NS  (2?),  the  LIF  scans  show  that  rotation- 
ally  and  vibrationally  excited  ground  state  NS  photofrag¬ 
ments  are  produced,  and  the  NS(C)  emission  indicates  a 


resonant  excitation  by  the  laser  on  an  NS (X-C)  transition. 
Further,  the  time  histories  of  the  emissions  suggest  the 
NS (B)  and  NS(C)  formation  rates  are  photolytically 
driven  and  are  independent  of  collisional  processes.  Given 
these  similarities  an  analogous  mechanism  is  proposed  for 
the  photodissociation  at  222  nm: 

c-S4N4  +  hr 222  -  NS  ( Ar, V )  +  g-S3N3, 

AT/—  —  32  kcalmol-1,  (R5) 

g-S3N3  +  hvl22  —  NS  (B,v)  -f  c-S2N2, 

AT/=  —35  kcalmol-1,  (R6) 

NS(Z,i;w=l)+/?v222  -NS(C,y,  =  2),  (R7) 

where  the  linear  dissociation  of  the  tetramer  [reaction 

(R5)]  is  supported  by  the  quadratic  dependence  of  the 
NS(C)  fluorescence  intensity  on  KrCl  excimer  fluence. 
The  energy  release  for  reactions  (R5)  and  (R6)  is  high 
enough  to  account  for  the  observed  NS(AT  and  NS(2?) 
vibrational  excitation,  i.e.,  NS(A",  v  =  4)  requires  13.5 
kcalmol-1  and  NS  (2?,  u  =  8)  requires  17.5  kcalmol-1. 
The  observation  of  NS  (2?,  umax  =  8)  again  restricts  the 
choice  of  the  S2N2  photofragment  to  the  cyclic  form. 

The  strong  variation  in  Trot  with  KrF  fluence  that  is 
depicted  in  Fig.  12  suggests  the  operation  of  a  fluence  de¬ 
pendent,  rotational  excitation  process  in  the  photolysis 
zone.  A  second  order  process  involving  an  efficient  V R 
energy  transfer  could  account  for  the  observed  behavior. 
Due  to  the  large  absorption  cross  section  of  S4N4  at  248  nm 
large  densities  of  NS  (A",  vu)  and  S3N3(u)  may  be  produced 
from  the  primary  photolysis  step  [reaction  ( R 1 )].  These 
species  could  react  further  in  energy  transfer  processes 

NS(y)  +  NS(i;,N)-NS(z;-l)  +  NS(»,N  +  /i/)  +  A£, 

(R8) 

S3N,(u)  +  NS(u,N)  -S3N,(u-  1 )  +  NS(u,N  +  w/)  +  A£, 

(R9) 

where  nJ  indicates  n  additional  rotational  quanta.  The  den¬ 
sities  of  vibrationally  excited  photofragments  will  depend 
on  laser  fluence.  For  high  fluence  photolysis  (T>75 
mJ  cm"2)  the  calculated  dissociation  is  near  unity.29  If  we 
assume  the  V-R  transfer  has  a  —20  ^s  lifetime,  i.e.,  the 
rotational  excitation  process  was  observed  at  a  20  /is  delay, 
and  that  the  high  fluence  photolysis  of  10  mTorr  S4N4 
produces  —  1  X  1015  mol  cm-3  of  NS(u)  -f  NxSx(y) 
photofragments  then  the  calculated  second  order  rate  con¬ 
stant  ~  10-11cm-3s-1.  This  V-R  rate  coefficient  is 
not  untypical  of  those  reported  by  Moore  for  a  large  class 
of  polyatomics,30 

Molecular  orbital  calculations,  both  semiempirical31 
and  ab  initio>n  indicate  the  S4N4  ground  state  to  be  a  sin¬ 
glet,  and  specify  the  intense  250  nm  absorption  band  as  the 
molecules  highest  occupied  molecular  orbital-lowest  unoc¬ 
cupied  molecular  orbital  (HOMO-LUMO)  transition; 
giving  it  a  2?2  —  2T  designation.  The  LUMO  is  a  degenerate 
IT  orbital  and  is  antibonding,  with  respect  to  the  S-S  link¬ 
age.  The  ab  initio  calculations  indicate  that  the  B2-E  en¬ 
ergy  separation  is  224  kcal  mol-1,  and  that  there  are  four, 
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nearly  resonant,  energy  states  (designated  by  *)  situated 
230  kcal  above  the  E  level  (LUMO).  Since  these  calcula¬ 
tions  have  not  incorporated  electron  correlation,  the  en¬ 
ergy  separations  are  clearly  too  large  (~2x);  this  may  be 
seen  by  comparing  the  measured  B2-E  A E  of  120 
kcal  mol-1  (obtained  from  the  250  nm  spectral  feature 
shown  in  Fig.  1 )  with  the  calculated  B2-E  A E  value  of  224 
kcal  mol-1.  Applying  a  correction  factor  of  X0.5  to  the 
calculated  values  yields  a  B2  to  *  state  A E  of  227 
kcal  mol"1,  which  is  in  good  agreement  with  the  fluence 
and  thermochemical  data  indicating  that  two  248  nm  pho¬ 
tons  (230  kcal)  are  required  to  reach  the  dissociative  state. 
The  calculations  suggest  that  a  two  photon  absorption  by 
S4N4  is  a  resonant  process  (sequential  excitation)  in  which 
the  first  photon  promotes  an  electron  to  the  E  antibonding 
n  orbital  (perhaps  with  rupture  of  the  weak  S  linkage), 
and  the  absorption  of  a  second  photon  promotes  the  elec¬ 
tron  from  the  E  state  to  one  of  the  four  upper  *  states.  It  is 
likely  that  the  *  states  are  repulsive  or  very  weakly  bound, 
and  rapidly  dissociate  to  yield  a  nonlinear  <z-S4N4  interme¬ 
diate.  The  temporal  histories  of  the  emissions  indicate  that 
the  lifetimes  of  these  states  is  r< 30  ns. 

V.  SUMMARY  AND  CONCLUSIONS 

Emission  from  several  electronically  excited  states  of 
NS  is  observed  when  the  energetic  molecule  S4N4  is  pho- 
tolyzed  with  radiation  from  an  excimer  laser.  Photolysis  at 
248  nm  generates  fluorescence  from  the  £(2nI/2t3/2), 
H2 n1/2,  G21~ ,  and  /2X+  states  of  NS.  NS(B  2nl/2i3/2) 
and  NS(C22+ )  fluorescence  is  observed  when  the  photol¬ 
ysis  wavelength  is  changed  to  222  nm.  The  NS (H)  and 
NS(C)  spectra  are  postulated  to  arise  from  a  resonant  in¬ 
teraction  between  the  KrF  and  KrCl  excimer  photons,  re¬ 
spectively,  and  vibrationally  hot  ground  state  NS.  LIF  ex¬ 
citation  scans  confirm  the  production  of  rotationally  and 
vibrationally  excited  NS(X)  up  to  v=4.  An  analogous 
mechanism  is  thought  to  generate  some  of  the  i/  =  2  bands 
of  NS(/,<7).12 

The  experimental  data,  thermochemistry,  and  ab  initio 
calculations  on  the  tetramer  support  a  nonlinear  (two  pho¬ 
ton)  dissociation  step  for  parent  S4N4  at  the  248  nm  wave¬ 
length.  It  is  postulated  that  the  dissociation  produces  vi¬ 
brationally  excited  NS  and  an  acyclic  S3N3  photofragment. 
The  observation  of  the  maximum  vibration  in  NS(^T),  i.e., 
y"=4  is  consistent  with  statistical  energy  partitioning  in 
the  photofragments  and  suggests  the  a-S3N3  segment  is 
formed  with  considerable  internal  excitation.  The  data  in¬ 
dicates  the  a~S3N3  fragment  may  undergo  further  photo¬ 
dissociation  leading  to  the  production  of  vibrationally  hot 
NS {B)  and  an  S2N2  molecule,  and  that  these  dissociations 
occur  on  a  time  scale  much  shorter  than  the  laser  pulse. 
Surprisingly,  no  evidence  for  the  formation  of  sulfur  or 
nitrogen  was  obtained.  This  may  be  contrasted  with  the 
explosive  or  thermal  decomposition  mechanisms  which  are 
known  to  generate  both  nitrogen  and  sulfur.33,34 

The  short  radiative  lifetime  measured  for  the  NS(C, 
u'  =  2)  state  (r<  10  ns)  may  be  indicative  of  some  predis¬ 
sociation  in  this  vibration.  The  NS(C)  intensity  distribu¬ 
tion  and  the  calculated  FCF’s  are  used  to  evaluate  the  C 


states  electronic  transition  moment  variation  with  intemu- 
clear  separation;  a  small  dependence  was  found  with  Re 
increasing  with  decreasing  r. 

The  photolysis  of  S4N4  at  excimer  wavelengths  may  be 
used  as  an  alternate  method  for  producing  the  NS  radical 
for  kinetic  and  spectroscopic  study.  In  comparison  to  the 
standard  method  for  NS  production  (reaction  of  active 
nitrogen  plus  sulfur  compound)  the  photolysis  is  largely 
free  from  interferences,  i.e.,  S2,  N2,  and  creates  a  relatively 
high  distribution  of  vibrationally  excited  ground  state.  The 
latter  effectively  red  shifts  the  absorption  spectrum  of  NS 
such  that  NS  electronic  states  with  energies  in  excess  of 
43  000  cm-1  are  more  readily  accessed  via  optical  excita¬ 
tion. 
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Second-harmonic  generation  was  produced  in  germanium-doped  silica  planar  waveguides  prepared  by  simulta¬ 
neous  illumination  with  1064-  and  532-nm  laser  light.  During  preparation  using  prism  coupling  to  specific 
waveguiding  modes,  the  film-generated  second-harmonic  intensity  grew  as  a  function  of  preparation  time  until 
it  saturated.  The  growth  rate  and  saturation  level  for  p-polarized  second-harmonic  intensity  was  an  order  of 
magnitude  greater  than  that  observed  for  the  s  polarization.  The  efficiency  for  a  2-cm  waveguide  length  was 
at  least  0.5%.  The  comparison  of  experimental  results  indicates  a  mechanism  for  this  planar  geometry  that  is 
similar  to  that  producing  harmonic  effects  in  optical  fibers. 

The  efficient  production  of  frequency-doubled  light 
in  optical  fibers  has  been  observed  through  optical 
processing  with  the  fundamental  frequency  co  (Ref.  1) 
or  additionally  seeded  with  frequency-doubled  light 
2(o. 2  The  process  is  attributed  to  the  formation  of  a 
spatially  varying  dc  field  within  the  fiber,  resulting 
in  quasi-phase-matching  between  co  and  the  fiber¬ 
generated  2<o. 3  Whereas  the  process  for  localized 
modification  of  optical  properties  is  not  fully  un¬ 
derstood,  it  has  been  linked  with  the  atomic  and 
molecular  structure  associated  with  dopants  and  de¬ 
fects  within  the  fiber.4"6  These  effects  have  been 
characterized  by  nonlinear  time-dependent  growth1,7 
followed  by  saturation  at  efficiencies  greater  than 
5%.8  Because  of  their  one-dimensional  nature,  no 
input  polarization  dependence  has  been  reported  for 
second-harmonic  generation  (SHG)  from  fibers8  or 
ridge  waveguides.9 

This  Letter  describes  the  observation  of  second- 
harmonic  light  generated  in  two-dimensional  wave- 
guiding  films  that  have  a  composition  similar  to  the 
optical  fibers.  We  present  measurements  of  second- 
harmonic  growth,  saturation,  and  erasure  in  specific 
waveguiding  modes  of  optically  processed  films.  In 
addition,  the  dependence  of  these  phenomena  on  the 
polarization  of  the  input  coupled  light  is  detailed. 

The  waveguides  were  produced  by  ion-beam 
sputtering  of  high-purity  silica  and  germanium 
onto  commercial  optical-grade  fused-silica  substrates. 

The  sputter  deposition  was  carried  out  by  an  argon- 
ion  beam  at  a  rate  of  0.11  nm/s  in  a  partial  pressure 
of  1  X  10"4  Torr  of  oxygen.  The  two  films  used  in 
this  study  were  approximately  2  and  4  /xm  thick 
and  had  an  average  surface  roughness  of  1.5  nm. 

Auger  analysis  indicated  that  the  germanium  content 
was  approximately  6  at.  %  and  was  substantially 
oxidized. 

The  experimental  geometry  is  shown  in  Fig.  1. 

The  film  was  illuminated  with  1064-nm  (<o)  and/or 


532-nm  (2 a>)  radiation  which  was  Q  switched  at 
1.22  kHz  and  mode  locked  with  an  average  power  of 
as  much  as  500  mW.  This  maximum  average  power 
produced  laser  pulses  with  a  peak  power  of  150  kW. 
A  half-wave  plate  controlled  the  polarization  of  both 
the  (o  and  2(o  light  with  a  polarization  ratio  of  greater 
than  98%  for  both  frequencies  and  polarizations. 
The  laser  light  was  focused  with  a  23-cm  focal-length 
lens  and  prism  coupled  into  the  film.  Scattering 
from  waveguided  light  at  the  film-air  interface  was 
collected  by  a  fiber-optic  waveguide  and  detected  by  a 
Hamamatsu  1P28  photomultiplier  tube  (PMT).  The 
waveguided  light  emitted  from  the  film  edge  was 
focused  by  a  10  X  microscope  objective  and  filtered 
for  detection  of  only  532-nm  light  by  the  PMT. 
Typically  300-1000  laser  pulses  were  detected  by 
the  PMT,  averaged  by  a  boxcar,  and  collected  by 
a  computer.  Larger  SHG  signals  were  detected 
by  a  photodiode  that  was  sensitive  to  oj  and 
2(o .  During  measurement  of  film-generated  second- 


to 

computer 


Fig.  1.  Experimental  setup  for  thin-film  waveguiding  ex¬ 
periments.  Laser  light  is  coupled  into  different  wave¬ 
guiding  modes  by  rotation  of  the  film  and  prism  assembly 
P.  To  optically  prepare  the  film,  filter  F  is  removed  to 
allow  illumination  of  film  by  (o  and  2<o  light.  Tb  measure 
film-generated  second-harmonic  light,  F  blocks  2 (o  light. 
Filter  FI  passes  only  second-harmonic  light  to  the  PMT. 
BS,  beam  splitter;  L,  lens;  FO,  fiber-optic  waveguide. 
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prism.  The  input  beam  was  s  polarized  with  a  532-nm 
wavelength. 

harmonic  intensities,  an  IR-transmitting  filter  was 
placed  before  the  sample  assembly  to  block  the  532- 
nm  seed  light. 

Waveguiding  was  initiated  by  coupling  a>  and/or 
2 o)  light  into  the  film  at  the  right-angle  corner  of 
the  prism-film  interface.  The  angles  for  maximum 
coupling  efficiency  were  determined  by  detecting  the 
intensity  of  scattered  light  from  waveguiding  modes 
as  a  function  of  the  external  coupling  angle  (Fig.  2). 
The  angular  coupling  ranges  for  the  lowest-order  co 
and  2(o  modes  were  0.64°  and  0.35°,  respectively. 
The  index  of  refraction  for  p -polarized  co  and  2 co 
illumination  was  calculated10  to  be  1.5234  and  1.5401, 
respectively,  by  using  the  coupling  angle  for  peak 
scattering  intensities.  Small  anistropies  (0.002) 
were  calculated  for  s-polarized  indices  that  were 
slightly  greater  than  the  experimental  error  (0.001). 
These  index  values  agree  with  the  dispersion 
measurements  obtained  by  using  visible  and  near- 
IR  transmission  measurements.11  The  waveguide 
scattering  intensity  as  a  function  of  distance  along 
the  waveguiding  path  indicated  that  losses  in  the 
waveguide  were  approximately  5  dB/cm. 

Initial  experiments  attempted  to  produce  harmonic 
generation  from  films  prepared  with  only  1064-nm 
illumination.  After  preparation  times  exceeding  8  h, 
no  second-harmonic  signal  was  observed.  In  our  film 
experiments,  incident  peak  laser  powers  of  as  much 
as  100  kW  were  incident  upon  the  coupling  prism,  but 
because  of  the  absence  of  an  air  gap  in  this  geometry 
the  maximum  light  coupled  into  the  waveguide  was  of 
the  order  of  1  kW.  The  majority  of  the  incident  beam 
was  reflected  at  the  prism— film  interface.  It  would 
appear  that  the  required  threshold  of  approximately 
5-kW  peak  power  to  initiate  efficient  second-harmonic 
intensities  with  only  the  fundamental  wavelength8 
was  never  reached. 

When  1064-  and  532-nm  light  was  simultaneously 
coupled  into  waveguiding  modes  of  a  germania-doped 
silica  film,  film-generated  second-harmonic  (FGSH) 
light  grew  as  a  function  of  time  after  initiating  illumi¬ 
nation,  as  shown  in  Fig.  3.  The  second-harmonic  sig¬ 
nal  was  detected  by  briefly  blocking  the  532-nm  seed 
beam  during  the  growth  process  and  measuring  the 


FGSH  signal  with  only  IR  illumination.  The  data  in¬ 
dicate  an  initial  quadratic  dependence  of  the  second- 
harmonic  intensity  for  approximately  20  min  before 
becoming  nearly  linear  for  as  much  as  4  h  when 
the  second-harmonic  intensity  saturates.  Multiple 
paths  within  the  film  were  prepared,  and,  in  general, 
saturation  occurred  between  2  and  10  h.  The  time  to 
reach  saturation  and  the  maximum  FGSH  intensity 
level  is  highly  dependent  on  the  intensity  I(a>)  and  the 
ratio  of  intensities  I{2(o)/I(co)  coupled  into  the  film. 
The  long  preparation  times  before  saturation  in  our 
thin  films  relative  to  fiber  preparation  times  were 
attributed  to  the  low  preparation  power  levels  within 
the  film,  which  are  due  to  low  coupling  efficiency  of 
the  laser  light. 

After  saturation,  the  FGSH  effect  was  stable  at 
room  temperature.  Prepared  films  retained  approx¬ 
imately  the  same  FGSH  efficiency  over  a  period 
of  10  days.  One  growth  study  was  stopped  and 
restarted  15  h  later,  with  the  same  efficiency  and 
growth  rate  as  before  the  interruption. 

To  ensure  that  the  second  harmonic  was  generated 
in  the  film,  an  output  prism  coupler  was  added  to 
intercept  both  the  a>  and  2  co  light  in  an  unprepared 
waveguiding  path.  After  the  output-coupled  2 co 
light  reached  saturation,  the  prism  was  removed,  and 
the  second-harmonic  light  was  then  observed  from 
the  edge  of  the  film,  as  shown  in  Fig.  1.  The  second- 
harmonic  light  again  grew  as  a  function  of  time 
and  saturated.  Because  the  only  change  was  the 
addition  of  a  longer,  unprepared  waveguiding  path, 
we  conclude  that  the  additional  second-harmonic 
light  was  generated  within  the  newly  illuminated 
waveguide  path. 

At  the  saturation  point,  the  efficiency  of  SHG  was 
determined  by  comparing  the  collected  output  power 
of  FGSH  with  the  fundamental  light  intensity.  After 
compensating  for  detector  sensitivity,  the  output  at 
the  fundamental  frequency  was  190  /iW  and  the  out¬ 
put  at  the  second  harmonic  was  approximately  1  yuW, 
which  gave  an  efficiency  of  approximately  0.5%.  The 


pled  into  the  film  for  preparation  was  approximately  1  kW 
for  the  1064-nm  light  and  5  W  for  the  532-nm  light. 
The  FGSH  intensity  was  measured  after  coupling  only 
1064-nm  radiation  into  the  film  with  a  peak  power  of 
0.8  kW. 
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Fig.  4.  Intensity  of  FGSH  versus  preparation  time  for  s- 
and  p-polarized  input  preparation  beams.  The  output- 
coupled  intensities  of  the  1064-  and  532-nm  preparation 
beams  were  the  same  for  the  two  polarizations. 


FGSH  light  could  be  observed  visually  near  the  sat¬ 
uration  intensities. 

Both  the  growth  rate  and  the  saturated  intensity 
levels  of  the  FGSH  light  produced  in  these  doped 
thin  films  displayed  a  dependence  on  the  polarization 
of  the  input  beams,  as  illustrated  in  Fig.  4.  The 
coupling  angle  for  these  data  was  chosen  so  that 
the  ratio  I(2co)/I(co)  was  the  same  for  both  s  and 
p  polarizations,  where  I(co)  ~  200  p,W  and  1(2  co)  « 
4  pW.  The  film  was  initially  illuminated  with 
s-polarized  co  and  2<o  light,  and  measurements 
were  made  as  described  above  for  time-dependent 
growth  measurements.  After  the  film  was  irradi¬ 
ated  with  the  s-polarized  beams  for  3  h,  a  half¬ 
wave  plate  was  adjusted  to  produce  p-polarized 
co  and  2(o  illumination.  The  change  in  the  po¬ 
larization  to  p-polarized  input  light  produced  a 
p-polarized  FGSH  signal  that  grew  an  order  of 
magnitude  faster  than  with  s-polarized  preparation 
light  and  a  corresponding  s-polarized  output.  The 
saturated  intensities  (not  shown)  of  the  FGSH  light 
were  also  an  order  of  magnitude  greater  for  the 
p- polarization  compared  with  the  s-polarized  light, 
but  the  time  for  the  FGSH  intensity  to  reach 
saturation  was  similar  for  both  polarizations.  This 
anisotropy  in  the  polarization  dependence  may  arise 
when  preferentially  ejected  electrons  predicted  by 
some  models12  are  effectively  confined  normal  to  the 
plane  of  the  waveguide,  while  those  ejected  in  the 
plane  would  diffuse,  further  producing  a  smaller 
effective  x 2*  This  model  would  require  either  long- 
lived  or  multiple  photoexcitations  to  move  electrons 
to  the  waveguide  boundary. 

For  the  experiments  described  above,  FGSH 
signals  were  observed  only  when  the  input  beams 
propagated  along  their  respective  lowest-order  wave- 
guiding  modes.  The  width  of  the  waveguiding  modes 
for  co  and  2<o  allowed  a  substantial  fraction  of 
the  power  at  the  optimal  coupling  angles  to  be 
propagated  simultaneously  in  the  film.  Typically, 
the  FGSH  light  was  observed  after  only  a  few  minutes 
of  initiating  propagation.  In  contrast,  after  30  min 
of  illumination,  no  FGSH  light  was  observed  for 
coupling  into  the  second-order  mode  of  the  1064-nm 
light.  The  lack  of  observed  FGSH  light  for  this  mode 


could  be  due  to  a  mismatch  of  the  optimal  coupling 
angles  for  the  co  and  2  co  light,  thereby  reducing  the 
intensity  of  copropagating  beams.  Alternatively,  the 
antisymmetric  nature  of  the  second-order  1064-nm 
mode  and  the  third-order  532-nm  mode  could  produce 
poor  overlap  and  thus  no  observable  harmonic 
intensity. 

In  fiber  experiments,  grating  erasure  was  demon¬ 
strated  by  a  number  of  techniques,  including  prop¬ 
agating  only  the  2 co  seed  beam  through  the  fiber.7 
We  attempted  to  erase  the  light-induced  grating  in 
the  film  by  initially  preparing  the  doped  film  with 
the  p-polarized  input  beam  in  the  lowest-order  IR 
waveguiding  mode  and  subsequently  illuminating  the 
film  with  only  2  co  light  after  removing  the  1064-nm 
input  beam.  The  green  seed  beam  illuminated  the 
film  at  various  2 co  coupling  angles.  After  illumina¬ 
tion  at  532  nm  for  periods  as  long  as  6  h,  the  initial 
preparation  angles  and  I(co)  conditions  were  restored, 
and  no  decrease  in  the  FGSH  was  observed. 

In  conclusion,  we  have  made  what  is  to  our  knowl¬ 
edge  the  first  observation  of  harmonic  generation 
in  two-dimensional  amorphous  doped  glass  films. 
These  results  in  thin  films  are  in  agreement  with 
observations  of  growth  and  saturation  of  harmonic 
generation  in  fibers.  The  anisotropic  response  to  the 
input  polarization  dependence  produced  by  the  thin- 
film  geometry  offers  new  insight  into  the  mechanism 
of  the  second-harmonic  effect.  For  future  studies 
these  films  offer  the  capability  of  more  controlled 
testing  of  the  optical  processing  because  of  the  ability 
to  excite  and  detect  specific  waveguiding  modes 
through  prism  coupling.  In  addition,  the  simplified 
geometry  allows  in  situ  characterization  that  can 
identify  the  microscopic  cause  of  these  effects. 

We  thank  the  Developmental  Optical  Facility  at 
Phillips  Laboratory,  Kirtland  Air  Force  Base,  New 
Mexico,  for  film  production  and  Scott  Urban  and  Skip 
Pribyl  for  their  assistance  in  data  accumulation. 
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Second  harmonic  generation  (SHG)  in  poled  fused  silica,  which  was  first 
reported  by  Myers,  Mukherjee  and  Brueck™  in  1991,  has  increased  in  importance 
because  of  a  multitude  of  potential  device  applications.  Such  applications  include 
electro-optic  switching,  frequency  doubling  crystals,  parametric  amplification,  parametric 
oscillation,  and  linear  electro-optic  modulation/frequency  conversion  which  can  be 
monolithicaily  integrated  into  optical  fibers  and  planar  integrated  circuit  geometries. 
Because  of  these  potential  applications,  much  work  has  been  done  to  understand  the 
phenomenon  of  SHG  in  poled  fused  silica. 

In  the  past,  a  second-order  nonlinearity  has  been  induced  in  fused  silica  by 
either  parallel  plate(1)  or  corona  poling.(2)  Through  parallel  plate  poling,  Myers  et  al.(1) 
were  able  to  induce  an  optically  active  region  of  thickness  5-10  pm  on  the  positively 
biased  face  of  a  sample.  The  second-order  nonlinearity,  x(2).  associated  with  this  region 
in  natural  fused  quartz  was  found  to  be  on  the  order  of  1  pnW.  But,  for  the  higher  purity 
UV  grade  synthetic  fused  silica  having  only  10  percent  as  many  impurities  as  natural 
fused  quartz,  the  magnitude  of  the  nonlinearity  was  found  to  be  an  order  of  magnitude 
lower.™  Based  on  this  data,  Myers  proposed  that  the  SHG  originated  from  a  space 
charge  region  created  by  migrating  impurity  ions  and  electrons.  This  was  further 
supported  by  the  ratio  of  electro-optical  coefficients,  r33/r3i,  which  was  found  by  Long  et 
al.™  to  have  a  value  of  3.0.  This  is  consistent  with  °omm  symmetry  and  a  space  charge 
region. 


Besides  parallel  plate  and  corona  poling,  nonlinearity  has  also  been  induced  in 
various  glasses  by  electron  implantation.  Kazansky  et  al.(4)  reported  second-order 
nonlinearities  as  large  as  0.7  pm/V  upon  implantation  of  lead  silicate  glass  with 
electrons.'  Russell  et  al.(5)  also  reported  the  generation  of  nonlinearity  in  lead  glass,  Ge 
doped  silica  glass,  Nb  doped  silica  glass,  Ti+Zr  doped  silica  glass,  lead  germanate  glass 
and  phosphate  glass  by  electron  implantation.  Although  nonlinearity  was  induced  in  the 
above  mentioned  materials,  attempts  to  induce  nonlinearity  in  fused  silica  by  electron 
implantation  have  failed.  In  fact,  Kazansky  et  al.(6)  reported  the  erasure  of  nonlinearity 
induced  by  thermal  poling  upon  electron  implantation. 

Because  of  the  apparent  role  of  the  space  charge  region  in  the  generation  of  a 
second  order  nonlinearity,  fused  silica  with  a  low-water  content  was  implanted  with 
protons  in  an  attempt  to  break  the  centrosymmetry  and  produce  a  space  charge  region. 
Protons  were  implanted  using  a  pelletron  at  energies  ranging  from  200  to  725  keV  and 
dosages  ranging  from  0.25  to  3.0  mC.  The  depth  of  the  implanted  region  increased 
from  approximately  1.5  urn  to  8.0  pm  at  proton  energies  of  200  and  725  keV, 
respectively.  The  level  of  SHG  was  studied  in  two  ways:  (1)  as  a  function  of  proton 
implantation  energy  at  a  constant  dosage  of  2  mC  (approximately  5x1 016  protons/cm2), 


and  (2)  as  a  function  of  dosage  (0.25  to  3.0  mC)  at  a  fixed  proton  energy  of  526  keV. 
Upon  implantation  of  the  desired  dosage  of  protons,  the  SHG  profile  was  obtained  by 
scanning  a  polarized  1.06  pm  Q-switched  Nd/YAG  laser  beam  at  an  incident  angle  of 
60°  across  a  diameter  of  the  irradiated  region  of  the  sample.  Prior  to  striking  the 
sample,  the  beam  was  reduced  to  between  10  and  60  mW  of  power  and  was  filtered  to 
remove  532  nm  radiation.  After  traversing  the  sample  and  generating  second  harmonic, 
the  beam  was  filtered  to  remove  all  radiation  except  wavelengths  of  530  ±10  nm.  After 
measurement  of  the  strength  of  the  frequency  doubled  radiation  with  a  PMT,  the 
absolute  level  of  d^  (d obtained  by  comparison  with  a  y-cut  quartz  standard. 

To  evaluate  the  distribution  of  optical  activity  in  the  irradiated  region  of  the 
sample,  linear  scans  were  performed  across  a  diameter  of  the  irradiated  area,  and  each 
point  recorded  was  the  time  averaged  response  of  a  photomultiplier  tube  to  the  SHG 
signal.  For  a  proton  implantation  energy  of  253  keV,  Figure  la,  we  found  the 
distribution  of  nonlinearity  induced  within  the  irradiated  area  to  be  a  maximum  at  the 
center  of  the  region,  falling  off  on  either  side.  But,  as  the  proton  implantation  energy 
increased,  the  strong  central  region  became  an  optically  active  annulus  of  increasing 
diameter  (Figure  1b).  The  central  region  which  had  no  visible  damage,  showed  very 
little  nonlinearity.  The  annular  region,  on  the  other  hand,  was  quite  optically  active.  The 
peaks  in  Figure  1b  correspond  to  the  linear  scan  passing  through  opposite  sides  of  the 
annulus.  The  optically  active  annulus  and  less  optically  active  central  region  are 
believed  to  originate  from  an  accelerator  beam  profile  which  became  annular  with 
increasing  diameter  as  the  implantation  energy  increased.  As  a  result,  the  portion  of  the 
sample  corresponding  to  the  beam  center  probably  received  very  little  irradiation  and, 
therefore,  had  much  less  optical  activity.  Kazansky  et  al. 5  observed  a  similar 
distribution  of  nonlinearity  in  lead  silicate  glass  that  had  undergone  electron  beam 
irradiation.  A  second  harmonic  signal  was  observed  at  the  edges  of  the  irradiated 
region  which  was  twice  that  in  the  center. 
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Figure  1  -  Relative  intensity  of  the  SHG  [a.u.]  vs  distance  [mm]  along  the  sample 
surface  for  a  cumulative  proton  dosage  of  2  mC.  Proton  implantation 
energies  are:  a.  250  keV,  b.  463  keV. 

The  strength  of  the  nonlinear  optical  signal  was  first  evaluated  as  a  function  of 
the  proton  implantation  energy  which  ranged  from  200  keV  to  725  keV  at  a  constant 
dosage  of  2  mC.  Upon  evaluation  of  the  nonlinear  optical  coefficient,  d-^,  for  the  various 
proton  energies  (Figure  2)  d^  was  found  to  be  dependent  on  the  proton  energy.  Each 


point  in  Figure  2  represents  a  different  sample.  Each  d33  value  was  calculated  from  the 
maximum  SHG  signal  found  along  a  single  linear  scan  through  the  implanted  region  of 
the  sample.  Because  a  linear  scan  through  a  different  line  of  the  implanted  arsa  might 
produce  a  slightly  different  result,  these  values  are  believed  to  be  good  within  ±0.1 
pm/V.  Saturation  of  d$3  in  the  range  0.5  to  0.6  pm/V  appeared  to  occur  for  proton 
energies  in  excess  of  410  keV.  One  exception  was  found  to  occur  for  a  sample  of  fused 
silica  implanted  with  568  keV  protons  where  a  dM  value  of  0.68  was  generated.  Finally, 
from  Figure  2,  a  minimum  proton  implantation  energy  of  approximately  350  keV  was 
found  to  be  necessary  to  generate  substantial  nonlinearity  in  fused  silica. 

One  conclusion  that  can  be  drawn  based  on  the  above  data  is  that  proton 
implantation  of  fused  silica  is  capable  of  creating  a  nonlinear  optical  region  exhibiting  a 
value  of  d33  at  least  as  large  as  obtainable  from  parallel  plate  poling  (i.e.,  approximately 
0.5  pm/V).  In  addition,  the  potential  exists  for  even  larger  nonlinearities  as  evidenced  by 
the  d33  value  of  0.68  pm/V  generated  by  a  proton  implantation  energy  of  568  keV. 
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Figure  2  -  Maximum  d^  [pm/V]  vs  proton 
implantation  energy  [keV]. 
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Figure  3  -  Max.  d33  [pm/V]  vs  dosage  [mC]  for 
a  proton  implantation  energy  of  526  keV. 


The  effect  of  proton  dosage  at  a  fixed  implantation  energy  of  526  keV  was 
studied  next.  This  particular  proton  energy  was  selected  because  of  its  relatively  large 
nonlinear  .response  at  a  dose  of  2  mC.  The  dosages  studied  ranged  from  0.25  mC  to 
3.0  mC.  As  is  apparent  from  Figure  3,  values  of  the  nonlinear  optical  coefficient,  d33, 
ranging  from  0.48  to  0.68  pm/V  were  obtained  at  dosage  levels  between  0.5  and  2.0 
mC.  When  compared  to  the  nonlinear  optical  coefficient,  d^,  produced  by  parallel  plate 
poling,  the  values  of  d33  produced  by  proton  implantation  using  these  parameters  are  up 
to  36  percent  greater.  Another  observation  apparent  from  Figure  3  was  that  as  the 
dosage  level  increased,  the  level  of  d33  achieved  decreased  after  a  certain  point.  This 
was  most  likely  due  to  extensive  lattice  damage  brought  about  from  prolonged  proton 
bombardment. 


In  conclusion,  the  distribution  of  SHG  on  the  surface  of  a  sample  of  irradiated 
fused  silica  was  found  to  become  an  annulus  of  increasing  diameter  as  the  energy  of 
the  proton  beam  increased.  This  damage  pattern  was  most  likely  the  result  of  an 
annular  beam  profile  brought  about  by  repulsion  of  the  positively  charged  protons.  By 
smoothing  out  the  beam  profile  or  moving  the  sample  during  irradiation,  it  should  be 


possible  to  induce  a  uniform  area  of  nonlinearity  stronger  than  that  achievable  through 
traditional  parallel  plate  poling. 

In  addition,  for  a  constant  dosage  of  2  mC,  a  threshold  proton  implantation 
energy  of  400  keV  was  found  to  be  necessary  to  generate  a  nonlinear  region  in  a 
sample  of  fused  silica.  As  the  proton  implantation  energy  increased  above  358  keV,  x 
(where  x(2>  =  sdsj)  appeared  to  generally  saturate  between  1 .0  and  1 .2  pm/V  for  proton 
implantation  energies  greater  than  400  keV.  One  exception  was  found  for  a  proton 
energy  of  568  keV  which  produced  a  value  of  x(2)  of  1 .4  pm/V.  For  a  proton  implantation^ 
energy  of  526  keV,  the  optimum  dosage  level  was  found  to  be  between  0.25  and  2.0" 
mC.  At  these  levels,  values  of  between  1.0  and  1.4  pm/V  were  found  to  occur. 
Finally,  values  of  x(2)  larger  than  1.4  may  be  achievable  by  finding  the  optimal 
combination  of  proton  implantation  energy  and  dosage. 


References 

1.  R.  A.  Myers,  N.  Mukherjee  and  S.  R.  J.  Brueck,  Opt.  Lett.  16  1732  (1991). 

2.  A.  Okada,  K.  Ishii,  K.  Mito  and  K.Sasaki,  Appl.  Phys.  Lett.  60  2853  (1992). 

3.  X.  -C.  Long,  R.  A.  Myers  and  S.  R.  J.  Brueck,  Opt.  Lett.  19(22)  1819-1821  (1994). 

4.  P.  G.  Kazansky,  A.  Kamal  and  P.  St.  J.  Russell,  Opt.  Lett.  18(9)  693-695  (1993). 

5.  P.  St.  J.  Russell,  P.  G.  Kazansky  and  A.  Kamal,  SPIE  (Photosensitivity  and  self¬ 

organization  in  optical  fibers  and  waveguides)  2044,  August  1 993,  Quebec, 
Canada. 

6.  P.  G.  Kazansky,  A.  Kamal  and  P.  St.  J.  Russell,  Opt.  Lett.  18(14)  1141-1142  (1993). 


Optical  Society  of  America  Annual  Meeting 
%j‘y  Portland  OR,  Sep  95.  Photosensitivity  and 

Quadratic  Nonlinearity  in  Glass  Waveguide's: 
Fundamentals  and  Applications 

First-Principles  Characterization  of  Structure  and  Properties  of  E’  Centers 

in  Silica  Glass 

Shashi  P.  Kama 

Frank  J.  Seiler  Research  Laboratory 
2354  Vandenberg  Dr.,  Suite  2A35 
United  States  Air  Force  Academy  CO  80840-6272 
Phone:(719)472-2655 
Fax:  (719)472-3649 
e_mail:  kama@gauss.usafa.af.mil 

> 

John  J.  Kester 

Frank  J.  Seiler  Research  Laboratory 
2354  Vandenberg  Dr.,  Suite  6H79 
United  States  Air  Force  Academy  CO  80840-6272 
Phone:  (719)  472-3122 
Fax:  (719)472-3649 

e_mail:  KesterJJ%FJSRL%USAFA@dfmail.usafa.af.mil 

Electronic  structure  and  properties  of  E *  center  and  its  precursor  in  a-Si02  have 
been  calculated  for  the  first  time  by  ab  initio  Hartree-Fock  method. 


v3 0 


First-Principles  Characterization  of  Structure  and  Properties  of  E'  Centers 

in  Silica  Glass 

Shashi  P.  Kama  and  John  J.  Kester 
Frank  J.  Seiler  Research  Laboratory 
2354  Vandenberg  Dr.,  Ste  2A35 
U.S.  Air  Force  Academy  CO  80840-6272 
Phone:  (719)  472-2655 
FAX:  (719)472-3649 

The  recent  observations  of  second  harmonic  generation  (SHG)  of  the  infrared  (IR)  light 
in  doped  silica  glass  fibers,1'2  planar  waveguides,3,4  and  poled  silica  glass5  have  attracted  a 
great  deal  of  attention  in  understanding  the  mechanism  and  origin  of  nonlinear  optical  (NLO) 
processes  in  a-Si02.  A  number  of  possible  mechanisms  have  been  proposed6,7  to  explain  the 
experimental  observations.  However,  in  so  far  as  these  phenomenological  models  describe  the 
observed  SHG  in  fibers  and  waveguides,  they  do  not  address  the  microscopic  mechanism  of 
optical  nonlinearity  in  silica  glass.  An  understanding  of  the  atomic-scale  mechanism  of  NLO 
properties  in  a-Si02  is  important,  not  only  for  increasing  our  fundamental  knowledge,  but  also 
for  modeling  new  materials  with  yet  enhanced  optical  susceptibility  and  faster  response  time  for 
device  applications.  Here  we  present  the  results  of  first-ever  ab  initio  studies  of  microscopic 
structure-NLO  property  relationships  in  systems  containing  free  spin  electrons,  such  as  a-S02. 

Considering  the  structure  of  a-Si02  as  a  continuous  random  network  of  Si-O,  its 
second  order  NLO  coefficient,  x(2)»  should  vanish  by  virtue  of  symmetry.  The  possibility  of 
observing  SHG  in  silica  glass  based  materials,  therefore,  requires  a  non-centrosymmetric 
structure  created  either  by  intrinsic  defects  or  by  external  dopants.  In  this  work,  we  concentrate 
on  characterizing  the  intrinsic  defect  and  its  possible  contribution  to  the  observed  NLO 
properties  of  a-Si02. 

Among  the  various  intrinsic  defects  identified  to  date,  the  E’  centers  constitute  the  single 
most  important  paramagnetic  point  defect  in  amorphous  silica.9  From  the  similarity  of  the 
electron  spin  resonance  (esr)  spectrum  of  the  y-ray  irradiated  silica  glass  with  that  observed  in 
radiation  induced  a-quartz,10  Griscom  concluded11  that  the  atomic-scale  structure  of  the 
paramagnetic  center  in  a-Si02  was  similar  to  that  of  E/.  The  current  accepted  structure  of  the 
Ei'  centers  in  a-quartz  is  based  on  the  Feigl,  Fowler  and  Yip  model.12  According  to  this  model, 
the  E| '  center  in  irradiated  quartz  is  produced  by  asymmetric  relaxation  of  =Si-Si=+  created  by 
an  O'  vacancy.  The  unpaired  electron  spin  is  located  on  a  tetrahedral  Si  atom  in  the  Si-0  (sp3 
hybridized)  bond  and  the  ‘hole’  containing  Si  (Si+  ion)  relaxes  to  sp2  hybridization,  containing  its 
three  SiO  bonds  in  the  same  plane.  The  two  types  of  Si  centers  in  E’  in  this  model  are 
characterized  by  the  two  distinct  (423  G  and  8G)  hyperfine  splitting  observed  in  the  esr 
spectrum  and  both  attributed  to  29Si  nuclei  (1=1/2). 

The  extensive  work  of  Griscom,9,11  aided  by  the  quantum  mechanical  studies  of  Fowler 
and  coworkers12,13  in  the  past  several  years  have  recognized  the  £/  the  Ei '  equivalent  in  silica 
glass.  It  is  now  also  believed  the  £/  centers  constitute  the  intrinsic  Frenkel  defects  (vacancy- 
interstitial  pairs)  in  amorphous  silica  created  by  melt  quenching  and/or  the  action  of  ionizing 
particle  irradiations.9 

In  order  to  establish  the  structure-NLO  property  relationship  of  £/,  we  have  calculated 
the  geometry,  one-electron  properties,  and  linear  and  nonlinear  optical  properties  of  four  model 
systems  by  the  ab  initio  Hartree-Fock  method.  The  model  systems  selected  (Fig.  1)  for  our 
study,  though  smaller  in  dimension  than  one  would  consider  adequate  for  modeling  purposes, 
contain  all  the  essential  structural  features  of  the  proposed  E/  center9  and  its  precursor.  The 


*3/ 


technical  advantage  of  selecting  the  four  model  systems  of  the  present  study  is  that  it  allowed 
us  to  perform  accurate  calculations  by  the  first-principles  quantum  mechanical  method. 

The  oxygen  vacancies  (ll-IV)  will  be  designated  as  follows:  Structure  II  will  be  denoted 
as  V0°,  structure  III  will  be  denoted  as  V0+1,  and  structure  IV  will  be  denoted  as  V0‘ .  The 
notation  is  such  that  V0  denotes  an  oxygen  vacancy  and  the  superscript  indicates  the  charge 
on  the  vacancy.  Accordingly,  a  positive  charge  indicates  a  ‘hole’  and  a  negative  charge 
indicates  an  ‘electron’.  The  precursor,  (structure  I)  will  be  denoted  as  P.  The  geometry  of  P 
and  the  three  vacancies,  V0°,  V0+1,  and  Vo'1  were  optimized  by  the  ab  initio  Hartree-Fock 
method  using  a  double  zeta  plus  polarization  (DZP)  Gaussian  type  basis  set.  For  the  property 
calculations,  the  DZP  basis  was  further  augmented  by  adding  a  set  of  diffuse  p  (0.025)  £nd 
d(0.025)  functions  on  Si,  diffuse  p(0.05)  and  d(0.05)  functions  on  O,  and  a  diffuse  p(0.1) 
function  on  H.  The  numbers  in  parentheses  are  the  orbital  exponents. 
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The  geometry  optimization  and  the  calculation  of  one-electron  properties  were  performed  by  the 
HONDO-8  suite  of  programs.14  The  linear  and  the  nonlinear  optical  properties  were  calculated 
using  the  recently  developed  XNLOP  program.15  The  XNLOP  program  is  a  general  purpose 
computer  code  for  calculating  the  static  (dc  field)  and  dynamic  (frequency-dependent)  linear 
and  NLO  properties  of  molecules  within  the  framework  of  the  time-dependent  coupled 
perturbed  Hartree-Fock  (TDCPHF)  theory.  The  details  of  the  TDCPHF  method  for  the  closed 
shell  systems-  are  described  in  the  works  of  Kama  and  Dupuis,16  Kama,17  and  those  for  the 
open-shell  systems  have  been  given  in  a  recent  work  by  Kama. 

The  calculated  bond  distance  and  bond  angles  of  the  systems  were  in  reasonable 
agreement  with  the  available  experimental18,19  and  the  most  recent  theoretical  values  For 
example,  the  Rsi^=1.633  A,  RSi-H=1.474  A,  <SiOSi=150  degrees,  and  <HSiO=109.5  degrees 
obtained  for  P  (structure  I)  in  the  present  study,  agrees  well  with  the  corresponding  values 
obtained  recently  in  a  similar  calculation20  and  also  with  the  experimental  data.  We  believe 
the  calculated  geometries  of  other  structures  (ll-IV)  also  have  similar  accuracy.  During  the 
geometry  optimization,  no  attempt  was  made  to  force  the  symmetry  on  the  structures.  This  had 
the  effect  that  in  the  systems  where  symmetry  would  cause  the  dipole  moment  to  vanish,  e.g.  in 
Vo°,  the  relaxed  geometry  leads  to  a  nonzero  value  of  the  dipole  moment. 

Table  1 .  Calculated  values  of  RSi-si  (A)  and  dipole  moment  (D) 


Structure 


Rsi-Si 

H* 

Hz 

iHl 


1 

II 

3.164 

2.367 

-0.00258 

0.01344 

-0.30653 

0.11267 

0.05827 

0.00923 

0.31204 

0.11384 

III  IV 

2.786  3.370 

0.00262  0.14997 

0.02290  0.14714 

-0.00109  0.02770 

0.02307  0.14997 


In  order  to  examine  the  geometry  relaxation  and  its  effect  on  calculated  properties,  the 
calculated  Si-Si  distance  and  the  dipole  moment  of  structure  I  -  IV  are  listed  in  Table  1.  It  is 
noted  that  a  neutral  O  vacancy  brings  the  two  Si  atoms  much  closer  by  about  0.08  A  with 
respect  to  the  precursor  (P).  A  negative  O  vacancy  (O’)  creating  the  Vo+1  structure  also  brings 
the  two  Si  atoms  closer  with  respect  to  P,  but  the  decrease  in  the  Si-Si  bond  length  is  now 
smaller  than  that  in  Vo0.  On  the  other  hand,  a  positive  O  vacancy  (O*)  that  leaves  behind  an 
electron,  thus  creating  Vo’1,  considerably  elongates  the  Si-Si  distance  with  respect  to  P. 

The  calculated  changes  in  the  Si-Si  distance  is  consistent  with  the  composition  of  the 
highest  occupied  molecular  orbital  (HOMO).  In  the  case  of  P,  the  HOMO  is  a  nonbondin tj  pz 
orbital  situated  on  the  oxygen  atom.  In  the  case  of  Vo°  the  HOMO  is  a  a-type  (sp  )  bonding  MO 
and  the  lowest  unoccupied  molecular  orbital  (LUMO)  is  its  antibonding  complement,  a  a-type 
(sp3)  orbital.  In  the  case  of  V0+1,  an  electron  is  removed  from  the  a  bonding  MO  of  the  neutral 
vacancy  which  leads  to  the  weakening  of  the  binding  between  the  Si  atoms  and  elongation  of 
the  Si-Si  bond.  In  the  case  of  Vo'1,  an  electron  is  added  in  the  antibonding  (o‘)  MO  of  the 
neutral  vacancy  which  destabilizes  the  bonding  between  the  Si  atoms,  resulting  in  a  large  Rsi-si- 
The  components  of  the  linear  polarizability  tensor,  a(-®;co),  and  first-hyperpolarizability 
tensor,  P(-2co;co,oo),  for  the  two  neutral  species,  P  (structure  I)  and  Vo°  (structure  II),  calculated 
at  a  fundamental  wavelength,  X=1.064  |im  are  listed  in  Table  2.  It  is  noted  that  whereas  the 
linear  polarizabilities  of  the  precursor  and  the  neutral  oxygen  vacancy  are  very  similar,  the  NLO 
polarizability  tensor  is  vastly  different  for  the  two  structures.  In  fact,  the  observable  vector 
quantity,  p,  (=  Zi  (Pijj+Pjij+Pjji)/5)  for  the  V0°  structure  is  vanishingly  small  because  of  its  nearly 
centrosym metric  structure. 

Table  2.  Components  of  a(-co;co)  (10  24  esu)  and  p(-2co;o),(o)  (10’32  esu)  for  P  (structure  I) 
and  V0°  (structure  II)  calculated  at  X= 1 .064  Jim 
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OCaa 

Obb 

OCcc 

Paab 

Pbbb 

Pbbc 

Pbcc 

Pccc 
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II 

8.27 

10.65 

7.35 

7.68 

7.40 

7.67 

5.85 

-0.24 

6.71 

-7.13 

-0.02 

10.27 

11.86 

7.75 

0.03 

-10.27 

Even  if  the  V0°  had  a  completely  non-centrosymmetric  structure,  as  would  be  expected  in  silica 
glass,  its  NLO  coefficient  is  expected  to  be  quite  small  due  to  the  fact  that  the  outermost 
electrons  in  its  valence  space  occupy  a  strongly  bonded  a  MO. 

Table  3.  The  hyperfine  coupling  constants,  ais0  (G),  and  Tg  (G)  for  Vo+1  and  Vo 1 


Struct.  Center  ais0  Txx  Tw  Tzz 

III  Si,  -17.77  15.84  15.68  -31.52 

Si2  -17.75  15.83  15.67  -31.50 

IV  Sii  -323.21  13.64  12.55  -26.19 

Si2  -323.74  13.75  12.55  -26.30 


The  29Si  hyperfine  coupling  constants  of  the  two  paramagnetic  species,  V0+1  (structure 
III)  and  Vo’1  (structure  IV),  calculated  by  spin-unrestricted  Hartree-Fock  method  are  listed  in 
fable  3.  It  is  immediately  clear  that  the  420  G  spectrum  observed  by  Griscom  could  not  be 


due  to  a  V0+1  center.  Such  a  large  hyperfine  splitting,  according  to  the  present  calculation 
would  require  excess  electron  spin  on  the  Si  nuclei,  as  may  be  possible  by  an  electron 
trapping,  resulting  in  a  Vo'1  (structure  IV)  type  structure.  Of  course,  the  smaller  hyperfine 
splitting  of  8G  observed  by  Griscom11  could  very  well  be  due  to  a  V0+  center  as  the  present 
calculation  does  yield  a  small  cJjS0  value  for  this  species.  In  fact,  the  value  of  ajS0  for  Vo 
calculated  after  removing  the  lowest  contaminating  component  from  the  spin  density  matrix, 
reduces  to  —1 1  G,  in  much  closer  agreement  with  the  experiment.  An  important  implication  of 
these  results  is  that  it  does  not  fully  support  the  hitherto  accepted  structure  of  the  E 
center.9,12,13  Rather,  the  present  study  indicates  an  ‘electron’-’hole’  pair  having  two  Si  nuclei 
with  distinct  electronic  environments  as  the  source  of  the  observed  esr  spectra.  A  more 
detailed  analysis  of  the  electronic  and  NLO  properties  of  the  model  E’  systems,  presently 
underway  in  our  laboratory,  will  be  able  to  provide  additional  information  about  the  structure  of 
this  important  defect  in  silica  glass. 
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The  generation  of  second  harmonic  light  in  optically  modified  fibers1  and  thin  films2 
has  been  shown  to  arise  from  a  periodic  modulation  of  the  second-order  susceptibility, 
%(2),  along  the  direction  of  propagation.3  This  longitudinal  modulation  allows  the  phase 
matching  of  the  generated  second  harmonic  light  with  the  fundamental  wave.  The 
efficiency  of  the  second  harmonic  for  different  waveguiding  modes  and  different 
polarizations  is  governed  by  an  additional  modulation  in  the  transverse  direction.  The 
overlap  of  with  the  transverse  profiles  of  the  optical  electric  fields  across  the 

waveguide,  as  described  by  the  overlap  integral,  governs  the  modal  response  of  the 
waveguide  generated  light/  Unlike  the  traditional  second  harmonic  generation  (SHG) 
processes,  the  overlap  integrals  describing  higher  order  modal  combinations  may  be  on  the 
same  order  as  the  overlap  integral  describing  the  combination  of  the  lowest  order  modes. 
In  planar  waveguides  the  difference  in  the  transverse  charge  density  governs  the 
polarization  dependence. 

Theoretical  predictions  of  second  harmonic  light  from  various  fiber  mode 
combinations  have  indicated  that  relative  efficiencies  within  an  order  of  magnitude  are 
possible  for  lower  order  modes.5  Furthermore,  they  found  the  fiber-generated  SHG  was 
produced  in  the  same  mode  as  that  used  during  the  seeding  process.  Recent 
experimental  work  with  modal  combinations  in  fibers  showed  that  while  the  seeding  was 
typically  carried  out  with  one  second  harmonic  mode,  the  fiber  generated  SHG  was 
reported  in  a  number  of  modes.6  In  one  experiment  the  SHG  was  not  observed  in  the 
second  harmonic  seeding  mode. 

We  report  here  on  the  modal  and  polarization  properties  of  SHG  in  germanium- 
doped  planar  waveguides.  While  there  has  been  experimental  observation  of  polarization 
effects  in  planar  waveguides2  no  theoretical  models  have  been  proposed.  The  film¬ 
generated  second  harmonic  (FGSH)  was  measured  for  different  modal  combinations  and 
polarizations  and  was  compared  to  the  theoretical  predictions  based  on  the  asymmetric 
photoionization  model.  The  cause  of  the  variation  in  the  efficiency  for  different  mode 
combinations  is  attributed  to  the  confinement  of  the  beams  within  the  waveguide  and  the 
overlap  integrals.  These  modal  combinations  can  be  used  for  possible  optical  switch 

applications.  . 

The  waveguides  were  produced  by  ion  beam  sputtering  and  consisted  of  6M 

percent  of  germania  in  silica.  The  experimental  setup  is  shown  in  Fig.  1.  A  fundamental 
beam  at  1.06  pm  was  produced  by  a  Q-switched  mode-locked  Nd:YAG  laser.  A  second 
harmonic  beam,  produced  by  KDP,  was  separated  and  recombined  with  a  slight  offset  to 
allow  prism  coupling  into  specific  waveguiding  modes  of  the  film  with  TM  or  TE  polarization. 
After  simultaneous  illumination  by  the  fundamental  and  second  harmonic  beam,  FGSH  light 
was  observed  by  reading  the  grating  with  a  fundamental  beam  having  a  power  equal  to  the 
power  of  the  writing  beam. 


splitter  A 


polarizer  rotation  stage 

Figure  1 .  Experimental  setup 


Films  prepared  with  both  beams  in  either  TM  or  TE  polarization  show  almost  three 
orders  of  magnitude  difference  in  FGSH  light  (Fig.  2).  This  difference  is  not  due  to 
differences  in  the  overlap  integral  between  the  two  polarizations,  but  rather  to  the  space 
charges  that  arise  due  to  the  asymmetric  photoionization  current7,  Jpi  ^  Em  E2m ,  and  the 
backfield  field  current,  J¥  =  a (1)  EDC,  where  a (1)  is  the  linear  conductivity.  At  saturation 

the  continuity  equation,  |^  +  V  /  =  0,  yields  V-7=0.  However,  the  most  general 

solution  is  not  Jpi  =  -  Jv,  but  rather,  V-  (Jpi  +J¥  )=0.  Only  by  adopting  the  general 

solution  can  one  explain  the  difference  in  FGSH  power  between  the  two  polarizations.  The 
photoionization  current  will  be  the  same  for  either  polarization,  however,  the  divergence  of 
the  current  is  much  larger  for  the  TM  configuration.  The  larger  divergence  in  the  current  for 
the  TM  polarization  will  yield  a  larger  charge  density  and  higher  effective  second  order 
susceptibility.  This  has  been  verified  through  numerical  simulations. 


Time  (minutes) 


Figure  2.  Growth  of  film  generated  second  harmonic  for  TM  and  TE  polarizations 

Using  only  TM  polarization  the  planar  waveguides  were  conditioned  using  the 
fundamental  and  second  harmonic  in  various  lower  order  mode  combinations.  For  each 
mode  combination  the  FGSH  was  allowed  to  grow  until  saturation.  The  output  FGSH 


powers  were  then  normalized  to  the  lowest  order  mode  combination  and  by  the  factor  1/ 
P©4  p2cd  f°r  the  seeding  beams  powers  (Table  1  .). 

When  the  film  was  conditioned  by  seeding  in  one  particular  mode,  the  subsequent 
FGSH  light  was  observed  substantially  in  the  same  order  mode.  Some  light  was  observed 
in  adjacent  modes.  However,  the  percentage  of  FGSH  light  observed  in  other  modes  was 
in  exactly  the  same  percentage  (±  0.2%)  as  that  observed  for  scattering  of  light  from  one 
mode  to  another  within  the  waveguide.  Thus,  there  is  no  appreciable  SHG  light  generated 
in  other  modes  that  is  in  excess  of  the  FGSH  scattered  light  from  the  mode  in  which  it  was 

conditioned.  .  - 

The  relative  power  generated  is  calculated  for  various  low  order  mode  combinations 
by  normal  mode  analysis.8  For  TM  modes,  the  dominant  component  of  the  electric  field  is 
in  the  transverse  direction,  i.e.,  normal  to  the  plane  of  the  waveguide.  In  the  case  where 
the  writing  and  reading  IR  beams  are  in  the  same  mode,  the  overlap  integral  describing  the 
FGSH  produced  in  the  same  waveguiding  mode  as  the  seed  green  is: 

G(z)=  j  jxSKr^+l/^lS  +  (l/3)r1l'Pul4  )<!>,/ +((r2 +2/3) 

— Q  — oo 

xl4>1,/«‘u  'i'i:,+(l/3)(2r1+r1)4>,,il'i'ul2'i'i;i  )®,Xi  <’> 

C2  /  3)^  WC  IJt  I2 1  <P, I2 1  'J'j.,  I2  }bdy 

where  the  IR  beams  are  assumed  to  be  in  mode  k,  the  seed  and  generated  green  are 
assumed  to  be  in  the  mode  /,  and  where  O  and  'F  represents  the  transverse  profiles  of  the 
transverse  and  longitudinal  components  of  the  optical  fields.  The  coefficients  r i  and  are 
obtained  from  experimental  results.9 

The  predicted  FGSH  powers  are  listed  in  Table  1  with  two  different  normalizations. 
The  column  marked  Total  Predicted  Power  calculates  the  power  in  various  modes  when 
normalized  to  the  total  power  propagating  in  the  waveguide  for  the  conditions  used  in  our 
experimental  setup.  This  column  can  then  be  compared  with  the  experimental  results 
observed.  A  very  rough  correlation  can  be  observed  between  the  experimental  and 
calculated  SHG  powers.  The  difference  is  primarily  attributed  to  variability  among  the 
mode  combinations  in  the  temporal  and  spatial  overlap  of  the  fundamental  and  second 
harmonic  used  for  conditioning  the  waveguide.  In  addition,  the  input  and  output  prism 
coupling  efficiency  varies  from  one  mode  combination  to  another.  However,  the  key  point 
is  that  SHG  can  be  produced  in  various  mode  combinations  and  that  they  have  been 
predicted  and  measured  to  be  within  an  order  of  magnitude  for  lower  order  modes. 

Table  1 .  Experimental  and  predicted  SHG  powers  for  various  mode  combinations 


Fundamental 

Mode 

Second  Harmonic 
Mode 

Measured 

Power 

Total  Predicted 
Power 

Waveguide  Predicted 
Power 

0 

0 

1.0 

1.0 

1.0 

0 

1 

0.39 

0.30 

0.31 

0 

2 

0.51 

0.36 

0.38 

1 

0 

0.10 

0.13 

0.28 

1 

1 

0.34 

0.40 

0.89 

1 

2 

0.58 

0.21 

0.48 

Contributions  from  modal  confinement  and  the  overlap  integral  were  examined. 
The  column  marked  Waveguide  Predicted  Power  normalizes  the  output  to  the  power 
•  confined  within  the  waveguide  itself.  This  column  excludes  power  that  is  outside  the 


waveguide  in  the  evanescent  fields.  When  comparing  the  calculated  total  power  and 
waveguide  power  it  is  observed  that  ratios  are  similar  for  a  specific  fundamental  waveguide 
mode.  When  the  fundamental  mode  increases  from  zeroth  to  first  order,  a  large 
discrepancy  is  observed  due  to  the  increased  electric  fields  outside  the  waveguide.  With 
better  confinement  conditions  the  SHG  powers  in  the  TMw,i  -  TM2co,i  could  be  up  to  90%  of 
the  TMoj.o  -  TM2(o,o  mode  combination. 

A  waveguide  was  then  modified  to  produce  multiple  gratings  within  the  same  spatial 
region  of  the  waveguide.  The  waveguide  was  modified  by  directing  a  fundamental  and 
second  harmonic  beam  into  the  film  as  described  above.  The  FGSH  was  allowed  to  grow 
to  a  specific  level.  The  same  region  of  the  waveguide  was  then  illuminated  by  directing  a 
fundamental  beam  into  a  different  waveguide  mode  along  with  the  second  harmonic  beam 
in  the  same  mode  used  previously.  There  was  no  FGSH  initially  observed  with  this  mode 
combination.  A  FGSH  signal  was  allowed  to  grow  from  this  new  mode  combination  to  a 
level  similar  to  the  first  mode  combination.  While  some  erasure  of  the  original  grating 
occurred  during  the  growth  of  the  second  grating,  second  harmonic  light  was  observed  in 
the  same  mode  (the  same  mode  originally  used  to  condition  the  film)  when  illuminated  by 
only  the  fundamental  beams  from  either  of  the  two  waveguiding  modes.  With  proper 
temporal  and  spatial  overlap  the  second  harmonic  beams  can  be  made  to  interfere,  which 
is  the  basis  of  an  optical  switch. 

In  conclusion  we  have  found  that  combining  different  waveguiding  modes  produces 
the  same  order  of  magnitude  FGSH  as  the  lowest  order  mode  combination.  Unlike  fibers, 
the  FGSH  was  produced  in  the  same  waveguiding  mode  as  the  seed  second  harmonic 
wave.  Furthermore,  we  found  that  TM  modes  produced  many  orders  of  magnitude  greater 
SHG  power  than  the  TE  modes.  This  is  due  primarily  to  the  greater  charge  density 
induced  with  TM  polarization,  and  appears  to  be  consistent  with  the  asymmetric 
photoionization  model. 
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Abstract 

Second-harmonic  generation  in  germanium-doped  silica  planar  waveguides  seeded  by  a  fun¬ 
damental  wave  and  its  second  harmonic  is  investigated  using  normal-mode  analysis.  The  effects 
of  self-  and  cross-phase  modulation  are  included  in  this  analysis.  A  set  of  two  coupled  equa¬ 
tions  that  describe  the  evolution  of  the  amplitude  of  the  writing  second-harmonic  beam  along 
the  direction  of  propagation  is  obtained.  The  solution  to  these  coupled  equations  is  used  to 
write  the  effective  x^  grating.  This  grating  is  then  used  to  determine  the  second-harmonic 
power  generated  when  the  beam  is  read  with  only  a  fundamental  wave  propagating  in  a  partic¬ 
ular  waveguiding  mode.  Due  to  diffraction,  the  conversion  efficiency  saturates  as  a  function  of 
waveguide  length.  It  is  also  found  that  if  the  reading  fundamental  wave  is  in  the  same  mode  as 
the  writing  fundamental  wave,  relatively  significant  SHG  is  obtained  with  the  generated  second 
harmonic  wave  propagating  primarily  in  the  same  mode  as  the  seed  second-harmonic  wave.  The 
power  generated  in  any  of  the  higher  order  modes  is  on  the  same  order  as  that  for  the  lowest 
order  mode. 

1  Introduction 

Efficient  doubling  of  1064nm  radiation  was  originally  observed  in  germanium-doped  silica  fibers 
[1].  The  fibers  used  in  the  early  experiments  were  illuminated  with  intense  infra-red  light  for 
several  hours  before  green  light  (532nm)  could  be  observed.  It  was  later  demonstrated  that  the 


i 
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preparation  time  can  be  reduced  considerably  if  the  fibers  were  additionally  treated  with  a  small 
amount  of  green  light  [2].  This  process  is  known  as  external  seeding.  Typically,  in  such  a  process, 
the  second-harmonic  power  grows  for  about  30  minutes  before  finally  saturating.  Efficiencies  of  up 
to  5  percent  have  been  achieved.  More  recently  researchers  have  been  able  to  observe  this  process 
in  germanium-doped  optical  planar  waveguides  [3]. 

Since  silica  is  centro-symmetric,  its  second-order  susceptibility  x^  vanishes.  Thus,  it  appears 
that  efficient  second-harmonic  generation  (SHG)  in  glass  should  not  be  possible.  Several  theories 
have  since  been  proposed  to  explain  this  surprising  phenomenon.  Although  this  phenomenon  is 
not  yet  fully  understood,  it  is  now  generally  agreed  upon  that  a  semi-permanent  dc  electric  field 
is  produced  via  a  third-order  or  higher  optical  process  in  the  fiber  [2,  4,  5,  6].  This  dc  field  breaks 
down  the  inversion  symmetry  and  produces  a  non- vanishing  x^-  A  non- vanishing  x^2\  however, 
is  not  sufficient  for  efficient  SHG.  Due  to  linear  dispersion  effects,  phase-matching  generally  cannot 
be  maintained  as  the  wave  and  its  second-harmonic  travel  along  the  waveguide,  and  thus,  even 
with  a  non-vanishining  x^  efficient  SHG  might  not  be  achieved.  Therefore,  it  appears  that  the 
induced  x^  is  also  periodic  in  the  direction  of  propagation  with  a  periodicity  matching  the  phase 
mismatch  between  the  fundamental  and  its  second  harmonic. 

The  infra-red  peak  power  used  in  the  fibers  and  planar  waveguide  experiments  is  usually  on  the 
order  of  10  KW,  with  the  green  seed  power  typically  being  two  orders  to  three  orders  of  magnitude 
less.  The  effective  transverse  area  is  typically  on  the  order  of  10"’11m2.  Therefore,  the  induced 
optical  dc  field  obtained  through  a  four-wave  mixing  process  is  on  the  order  of  a  few  hundred 
volts  per  meter.  It  can  be  shown  that  the  required  induced  semi-permanent  dc  field  is  a  few 
orders  of  magnitude  higher  than  the  optical  dc  field  [7].  Researchers  were  able  to  generate  efficient 
SHG  by  applying  an  external  periodic  dc  electric  field  to  both  germanium-doped  silica  fibers  and 
planar  waveguides  [8,  9].  The  required  external  field  was  on  the  order  of  a  few  MV/m  [10].  It 
was  therefore  hypothesized  by  some  researchers  that  the  nonlinear  interaction  of  the  optical  fields 
causes  some  form  of  charge  separation  that  leads  to  a  build-up  in  an  internal  dc  electric  field  which 
eventually  saturates.  One  promising  model  that  explains  this  charge  separation  is  the  asymmetric 
photoionization  model  [11,  12,  13].  In  this  model,  one  and  two  photon  absorption  processes  interfere 
resulting  in  electrons  being  ejected  in  a  preferential  direction.  In  its  simplest  form,  the  asymmetric 
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photoionization  current  is  given  by: 

j  oc  E*2E2u  +  c.c.  (1) 

The  resulting  dc  electric  field,  Edc,  then  interacts  with  the  optical  fields  to  produce  SHG  through  a 
four-wave  mixing  process  characterized  by  X(3)(-2u>;  0,o>,o;),  where  X(3)  is  the  third-order  suscep¬ 
tibility.  Therefore,  the  induced  % ^  is  given  by: 

x(2)  =  3x(3)^  (2) 


We  examine  here  this  process  of  second-harmonic  generation  in  thin-film  planar  waveguides  for 
TM  launched  modes.  This  is  motivated  by  experimental  results  that  showed  much  higher  SHG 
efficiencies  for  TM  launched  modes  than  those  obtained  for  TE  launched  modes  [3].  Thus,  the 
induced  dc  electric  field  builds-up  more  effectively  normal  to  the  waveguide  (along  the  optic  axis) 
than  in  the  plane  of  the  waveguide.  Using  normal-mode  analysis,  and  including  the  effects  of 
self-phase  and  cross-phase  modulation,  we  compute  the  efficiency  of  generating  different  waveguid- 
ing  modes  for  various  combinations  of  infra-red  and  green  waveguiding  modes  during  the  writing 
process.  Infra-red  light  in  different  modes  is  also  used  during  the  reading  process. 


2  Derivation  of  amplitude  equations 

We  will  now  derive  the  equation  that  describes  the  evolution  of  the  complex  amplitude  of  the  second- 
harmonic  light  using  normal-mode  analysis.  It  should  be  pointed  out  that  we  are  considering  a 
time-independent  model,  and  thus  we  do  not  address  the  saturation  of  this  process  in  time.  We 
note,  however,  that  a  normal-mode  analysis  will  offer  a  considerable  reduction  in  computational 
time  over  more  numerically-based  techniques.  This  is  crucial  for  the  time- dependent  case  in  which 
coupled  equations  for  the  induced  dc  electric  field  and  the  optical  fields  have  to  be  solved  at  each 
time  step,  as  the  computations  could  prove  to  be  exhaustive. 

The  method  used  here  is  as  follows:  Using  a  simple  model  to  describe  E we  will  initially  solve 
for  the  amplitude  of  the  writing  fundamental  and  second-harmonic  waves  launched  in  particular 
modes,  and  then  use  the  result  to  determine  the  x^  grating.  Once  the  grating  is  written,  the 
power  of  the  green  light  generated  in  the  different  waveguiding  modes  when  it  is  read  with  only  an 
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infra-red  beam,  also  launched  in  a  particular  mode,  is  computed.  It  is  assumed  that  the  generated 
second-harmonic  wave  does  not  affect  the  x^  grating. 

We  consider  a  germanium-doped  silica  planar  waveguide  sandwiched  between  two  semi-infinite 
dielectric  media,  both  possessing  zero  second-order  susceptibility.  The  direction  of  propagation  is 
taken  to  be  in  the  z-direction.  The  optic  axis  is  perpendicular  to  the  direction  of  propagation  and 
is  taken  to  be  in  the  y-direction.  The  wave  diffracts  in  the  unconfined  x-direction  and  we  assume 
that  such  diffraction  has  a  fundamental  Gaussian  profile. 


Using  normal  mode  analysis,  one  can  express  the  electric  field  for  TM  modes,  e^-,  oscillating 
at  a  given  frequency  a;*,  i.e.  Ui  =  u  and  u2  =  2a;,  and  propagating  in  a  particular  waveguiding 
mode  j  as: 


e,j(x,i)  =  (l/2)Ei,i(x)e  ,Uit  +  C.c.  (3) 

The  spatial  parts  of  the  field  components,  and  are  given  by: 

$8  M  =  <«> 

bW(x)  =  AiAWUvy U<*.  »)«*•'*  (5) 


where  accounts  for  the  modulation  of  the  amplitude  along  the  direction  of  propagation  due 
to  nonlinear  interactions  as  well  as  linear  losses,  and  Tpij  describe  the  distribution  of  the  y 
and  z  components  of  the  field  along  the  optic  axis,  respectively,  Oij  describes  the  profile  of  these 
components  along  the  unconfined  x-direction  at  a  given  point  z  and  accounts  for  the  diffraction  of 
the  wave,  and  fiij  is  the  propagation  constant.  Typically,  in  optical  planar  waveguides  E M  E^v\ 

The  field  distribution  in  the  y-direction  for  e\^  is  given  by: 


Mv)  =  5exP 


-\lPh-nhkl(y-°) 

c  cos  [\Jnlokl  -  fiijy  -  s^j 
=  \/ Pi,j  -  nhkl  {y  +  a) 


y  >  a 

—  a  <  y  <  a 

y  <  -a 


(6) 


where  the  waveguide  is  taken  to  lie  between  —a  and  a,  ki  =  Wi/c  is  the  free  space  wave  number, 
and  where  n^o,  n^i,  and  are  the  indices  of  refraction  at  frequency  of  the  waveguide,  the 
dielectric  above  the  waveguide,  and  the  dielectric  below  the  waveguide,  respectively.  The  coefficients 
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B,  C,  D  and  the  phase  shift  Si  are  determined  from  the  boundary  conditions.  The  profile  in  the 
y-direction  of  the  z  component  of  the  field,  ipij,  is  related  to  faj,  by:  iphJ  =  dy4>ij.  The  profile 
in  the  x-direction,  is  taken  to  be  a  fundamental  Gaussian  slowly  varying  in  the  z-direction, 


<  2 iPi,jdz6itj  <  PljOij,  and  is  given  by: 


Oi,j(x>  z)  — 


1  +  (2  iz/Pijwlj) 


1/2  exp  [wl.  +  (2iz/Pij)\ 


—  X * 


(7) 


where  represents  the  spot  size  of  the  beam  in  the  x-direction  at  z  =  0. 


The  fundamental  wave  is  treated  in  the  undepleted  pump  limit.  Since,  the  seed  second-harmonic 
power  is  much  smaller  than  that  of  the  seed  fundamental  wave,  self-phase  modulation  (SPM)  dom¬ 
inates  over  cross-phase  modulation  (XPM).  For  a  seed  fundamental  wave  launched  in  a  particu¬ 
lar  mode  k ,  and  in  the  limit  of  a  slowly- varying  amplitude,  ( d2Aitk/dz 2)  <  2ij3iik(dAlik/ dz)  < 
coupled  mode  theory  gives  [14]: 


2if3ltke^z 


dM,k 

dz 


oo  oo 

+  YAlik  J  J  |$1.fc(x>y>2:)|  dxdy  = 


+ 


—jr-2  J  J  p*(x»w)' 

— oo  —  oo 

,k(x>y>z)z]dxdy  (s) 


where  for  brevity  we  use  $1^  =  (/>itk8itk  and  ai  is  the  power  attenuation  constant 

for  the  fundamental  wave  due  to  linear  losses,  and  where  P5(x,u;)  is  the  Fourier  transform  of  the 
nonlinear  polarization  due  to  SPM.  Coupling  to  other  waveguiding  modes  through  SPM  is  for  most 
cases  small  and  is  ignored  here.  The  double  integral  on  the  left  hand  side  of  the  equation  above  is 
independent  of  z.  Double  integrals  of  this  form  are  henceforth  represented  by  Kij: 


00  00 


dxdy 


(9) 


Due  to  the  amorphous  nature  of  glass,  the  non- vanishing  components  of  the  third-order  tensor 
Xijh  are  °f  the  form  xj;2  and  x\ijj  including  all  permutations  of  the  indices  i  and  j.  Furthermore, 
these  non-vanishing  components  are  related  by: 

Xiijj  =  =  Xijji  =  (!/3)xS!  (10) 


Y(?)  =Y(3.). 
Ann  Ajjjj 


(il) 
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Therefore,  the  polarization  responsible  for  SPM  is  given  by: 


if  (*,») = + |K2iVi + !(s,(2)>8)’) 


The  solution  of  Eq.  (8)  is  then  given  by: 


M,k{z)  =  i4lifc(0)  exp  (~YZ  +  ir]l  J  Mz')\A^,k(z  )\2dz' 


I  A1>k(z)\2  =  |Alifc(0)|2e—  (14) 

where  771  =  (3uj2) / (8c2 Ki^Pi^) ,  and  fi(z)  is  the  overlap  integral  given  by: 

oo  oo 

fl(z)  =  //  xiS(l*i,*|4  +  1*1,  fcl4  +  (4/3)|$1>Jfe|2|$1>fc|2  +  (2/3)^lk^k)dxdy  (15) 


For  the  second-harmonic  wave  seed  ,  SPM  is  much  smaller  than  XPM.  For  a  wave  launched  in 


a  particular  mode  Z,  coupled  mode  theory  then  gives: 


+ fA»\  =  -5?  /  /  Nx'2"> 


‘^(x, 2o>)j  •  ]$*2l(x,y,z)y 


*2  i(x,y,z)z\dydx 


where  $2 ,1  —  <t>2ii®2 tl  and  =  ^2,^2,*  >  «2  is  the  power  attenuation  constant  for  the  second 
harmonic  wave  due  to  linear  losses,  Pc(x,  2u;)  is  the  Fourier  transform  of  the  nonlinear  polarization 
responsible  for  XPM,  P3hg(x>  2a; )  is  the  Fourier  transform  of  the  nonlinear  polarization  responsible 
for  SHG.  Coupling  to  other  waveguiding  modes  is  small  and  is  neglected  here.  Since  the  effect  of 
dispersion  on  is  small,  the  polarization  responsible  for  XPM  is  given  by: 

/*>(*,  2a.)  =  |e„x§  (K!|!4:.)+  fK’iN?)  (1T) 

For  Pjh5(x,  2o/),  The  effective  second-order  susceptibility  tensor  xjjl  is  given  by: 

xgl  =  3xg,££  (18) 

where  E^J  is  the  induced  static  electric  field.  We  assume  here  that  the  static  electric  field  is 
effectively  built  up  along  the  optic  axis,  i.e.  the  y-axis.  Furthermore,  we  model  the  component  of 
Edc  along  y  in  terms  of  the  optical  fields  as  [15]: 
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where  6  is  a  constant  [12];  possibly  a  complex  constant,  and  where  T\  and  r2  are  also  constants. 

(2)  (2)  (2) 

The  non- vanishing  components  of  the  induced  second-order  susceptibility  are  then  Xyvv>  Xy zz>  Xzyz> 
and  xi*y  •  These  components  are  related  by: 

xS  =  xg  =  xg  =  (l/3)xS  (2°) 


Strictly  speaking,  the  asymmetric  photoionization  process  involves  a  third-order  conductivity  tensor 
aijkl  ^at  has  the  same  symmetry  properties  as  Xijh  except  that  the  Kleineman  symmetry  does  not 
hold  here  since  photoionization  is  believed  to  be  a  resonant  process.  Experimental  results  indicate 
that  ct^/43,!  ~  0.64  [16].  Since  from  symmetry  =  a one  concludes  that  T\  «  0.64  and 
r2  ~  0.36.  It  should  be  pointed  out  here  that  a  more  exact  treatment  entails  solving  a  time- 
dependent  problem,  with  Edc  being  obtained  from  the  static  charge  distribution  which  is  governed 
by  the  charge  conservation  equation. 


The  evolution  of  the  complex  amplitude  of  the  second-harmonic  wave  along  the  direction  of 
propagation  is  then  given  by: 


dA2i 

-IT  = 1712 


+  8f2(z)\Alik(z)\>  +  bf3(z)\Alik(z)\'1  J  A2,i  +  b*  fA{z)A\M{z)A'2>le 


*  —2iA/3z 


(21) 


where  r)2  =  (3a;2)/(8c2K2,f/32ij),  A/3  is  the  phase  mismatch  given  by  A/3  =  /32,j  —  2/3i f2  is  the 
overlap  integral  for  XPM: 

oo  oo 

Mz)  =  II  {(\*i,k\2  +  (l/3)l*i,fc|2)  1*2, l|2  +  (l^i.fcl2  +  (1/3)1$!, fe|2)  1^2, i|2]  dxdy 

—  oo  — OO 

and  where  /3 (z)  and  /4(z)  are  the  overlap  integrals  for  for  this  particular  process  of  SHG: 

a  oo 

/a(^)  =  //  x![(l*i*r  +  (n  +  l/3)*;,^:5i  +  (l/3)nl*u|4)l#ai,|2+  ((r2  +  2/3) 

—a  —  OO 

X  |$i,fc|2$iiJt¥i)fe  +  (l/3)(2n  +  r2)$lifc|$1>fc|2$t>fc)$2),$;iJ 
+  (2/3)r2|$1,*|2|$1,fc|2|$2,i|2]dxd3/ 


a  oo 

u (z)  =  J  j  x£![($tfc  + 

—a  —oo 

+ 


(n  +  l/3)$2>fc*2ifc  +  (l/3)r!^ifc)^  +  ((r2  +  2/3)$?|Jfe¥1)ife 
(l/3)(2n  +  r2)$  +  W)r2^lk^dxdy 
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Note  that  f2(z)  and  /3(z)  are  real,  while  /4(z)  is  complex.  Also  note  that  the  limits  of  integration 
for  the  overlap  integrals  describing  SHG,  f3{z)  and  /4(z),  are  -a  and  a  as  the  effective  x(2)  grating 
is  written  inside  the  waveguide  with  the  two  dielectrics  on  the  outside  both  possessing  zero  second- 
order  susceptibility. 

The  equation  describing  the  complex  conjugate  of  the  amplitude  is  then: 

(-^+s/2(^i'4i.‘(*)|2+iV3(»)M.,4(^i4)^+V4,w-4;jb(^2,^‘a'’'  (22) 

If  one  then  uses  the  transformation: 

A2ii(z)  =  A'2  l(z)exp^-Yz  +  i'n2 /[8/2(,')|Alifc(z')|  +  bf3(z)\Ahk(z  )|  dz  j  (23) 

Eqs.  (21)  and  (22)  can  be  expressed  as: 

dA' 

=  imb'f^AitWA^e-W  (24) 

^  =  -imbfl{*)A{%{z)A!2/«*)  (25) 


C(z)  =  2  (a/3z  +  t]2  J  [s\Ahk(z')\2  f2(z)  +  bf3(z)\Alik(z)\4]dz 


where 


Eqs.  (24)  and  (25)  can  be  readily  decoupled  to  obtain: 

^  +  (‘f  “  ^  -  vlV’WlAi.I.Uv  =  0  (26) 

The  solution  to  Eq.  (26)  is  then  used  along  with  Eqs.  (13),  (18),  and  (19)  to  determine 
the  written  grating.  We  note  that  the  general  solution  to  Eq.  (26)  allows  for  some  film 
generated  second  harmonic  during  the  writing  process,  and  the  equations  used  to  determine  the 
X^2)  grating  take  this  into  account.  However,  if  b  is  real  and  b\AXyk\A  is  sufficiently  small,  the  seed 
second  harmonic  wave  undergoes  a  phase  shift  but  no  growth  along  the  direction  of  propagation. 
Neglecting  self-  and  cross-phase  modulation,  and  linear  losses,  one  can  readily  show  in  the  plane 
wave  limit  that  this  holds  true  for  (u>x^&|A|4)/c  <  Afc,  where  A k  is  the  wavenumber  mismatch. 

The  grating  is  then  read  with  only  a  fundamental  wave  propagating  in  a  particular  mode  m. 
The  spatial  parts  of  the  field  components  for  the  generated  green  are  expressed  as: 


4V>W  =  E^.A MK i(*,*)e««' 


(27) 
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^2z)(x)  =  Y1  A2,J'tp2,J(y)e2iJ(x,  zyh’lZ 


(28) 


where  we  sum  over  all  possible  green  waveguiding  modes.  Applying  normal  mode  analysis  again,  one 
can  obtain  the  equations  describing  the  complex  amplitudes  of  the  reading  fundamental  wave  and 
the  generated  second  harmonic.  We  neglect  the  pump  wave  depletion  due  to  SHG.  This  is  justified 
since  the  maximum  conversion  rate  is  on  the  order  of  a  fraction  of  one  percent.  The  amplitude  for 
the  reading  fundamental  wave  is  described  by  an  equation  similar  to  Eq.  (13),  with  the  overlap 
integral  involving  the  waveguiding  mode  m.  For  the  generated  second  harmonic  propagating  ,  for 
example,  in  mode  n,  the  evolution  of  the  complex  amplitude  along  the  direction  of  propagation  is 
given  by: 


dz 


^3 


( 


^  +  |A,,m(.)|  !*/£(*)  W  +  (1/3 ilfflV*'-*'! 


(29) 


where  773  =  (3a;2)/(c2K2,n/32,n)>  Ki2)  ^as  same  form  as  /2(z)  except  that  it  involves  mode  m  for 
the  fundamental  wave  and  mode  n  for  the  second  harmonic,  and  where  the  overlap  integral  f${z) 
is  given  by: 

o  00 

/»(*)  =  /  /  x£ +  (2/3)$!,mf 

—a  —00 

Using  the  initial  condition  A2jTI  =  0  at  z  =  0,  one  can  readily  obtain  the  solution  to  Eq.  (29): 


A2,n  =  J  fs(z)Alttn(z)  exp  [$(*')  -  i{/32,n  -  2Pi,m)z 

0 

where  ^(z)  is  given  by: 

z 

t(z)  =YZ~  iV3  /  ^(2,)l^l,m(2')|2<iz' 


dz 


(30) 


Results  of  the  numerical  simulation  are  presented  in  the  next  section. 


3  Results 

The  problem  is  simulated  numerically  by  launching  an  infra-red  fundamental  beam  along  with  a 
green  second  harmonic  beam  in  a  3  cm  long,  2  fim  deep,  germanium-doped  waveguide  sandwiched 
between  silica  on  one  side  and  air  on  the  other  side.  The  indices  of  refraction  for  the  fundamental 
wave  for  the  waveguide,  silica  and  air  are  1.532,  1.4496,  and  1.0,  respectively.  The  corresponding 
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indices  of  refraction  for  the  second-harmonic  are  1.546,  1.4607  and  1.0.  Linear  losses  for  the  infra¬ 
red  wave  and  its  second  harmonic  are  taken  to  be  0.5  dB/cm  and  1.0  dB/cm,  respectively.  For  all 
beams,  reading,  writing  and  generated,  the  spot  size  is  taken  to  be  30/xm.  Except  where  noted, 
both  the  power  of  the  infra-red  writing  beam  and  the  power  of  the  green  writing  beam  are  kept 
constant,  with  the  power  of  the  reading  infra-red  beam  set  equal  to  that  of  the  infra-red  writing 
beam  (20  KW  peak  power).  The  power  of  the  green  writing  beam  is  taken  to  be  three  orders 
of  magnitude  less  than  that  of  the  infra-red  beam.  For  a  given  mode  oscillating  at  a  particular 
frequency,  the  power  is  given  by: 

oo  oo 

PiJ  =  (1/4)  J  J  Em-  x  H Ijdxdy  +  c.c.  (31) 

— oo  — oo 

The  constant  b  in  Eq.  (19)  is  assumed  to  be  real  and  its  value  is  adjusted  to  provide  a  maximum 
Edc  of  about  107V/m  inside  the  waveguide. 

We  launch  both  infra-red  and  green  seeds  into  the  waveguide  in  several  different  waveguiding 
TM  modes  and  compare  the  results  of  the  different  cases.  The  coupled  equations  (24)  and  (25)  are 
then  integrated  using  a  fourth-order  Runga-Kutta  method.  This  is  more  numerically  accurate  than 
solving  the  second-order  differential  equation  (26).  For  the  parameters  used  in  our  simulations,  the 
SHG  polarization  leads  to  a  phase  shift  in  the  second  harmonic  seed  (see  the  discussion  below  Eq. 
26).  The  resulting  x ^  grating  is  then  used  during  the  reading  process.  Different  TM  modes  are 
also  used  during  the  reading  process.  The  amplitude  of  the  generated  green  is  given  by  Eq.  (30) 
and  is  computed  for  all  possible  waveguiding  modes. 

In  all  numerical  simulations,  a  relatively  significant  amount  of  SHG  is  obtained  when  both 
the  reading  and  writing  infra-red  beams  are  launched  in  the  same  TM  mode.  Furthermore,  the 
generated  green  is  found  to  propagate  mainly  in  the  same  mode  as  the  seed  green.  All  other 
combinations  generate  a  few  orders  of  magnitude  less  power  with  a  rapidly  varying  sin2  length 
dependence.  As  an  example,  Figure  (la)  shows  the  growth  of  the  second-harmonic  power  generated 
in  the  TMo  mode  as  a  function  of  distance  along  the  direction  of  propagation,  z,  for  the  case  of  a 
waveguide  seeded  by  infra-red  and  green  beams,  both  launched  in  the  TMo  mode,  and  read  with  an 
infra-red  beam  also  launched  in  the  TMo  mode.  Note  that  due  to  linear  losses  and  diffraction,  the 
power  generated  is  not  quadratic  in  length  but  rather  reaches  a  maximum  value  and  then  falls  off  as 
the  SHG  power  obtained  from  the  decreasing  x^  grating  becomes  less  than  the  linear  loss.  If  there 
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are  no  linear  losses,  the  power  then  saturates  due  to  diffraction.  This  case  is  shown  in  Figure  (lb). 
It  should  be  pointed  out  that  the  SHG  power  actually  fluctuates  very  slightly  with  a  periodicity  of 
about  10.5 fim.  This  is  the  quasi-phase  matching  period  for  this  combination  of  infra-red  and  green 
light  modes.  These  observations  can  be  explained  by  the  phase-matching  condition  required  for 
efficient  SHG.  In  the  cases  in  which  relatively  poor  SHG  is  obtained,  the  phase-mismatch  between 
the  reading  infra-red  and  the  generated  green  is  not  cancelled  out  by  the  phase  variation  of  the 
effective  second-order  susceptibiliy  x ^  induced  during  the  reading  process.  It  should  be  pointed 
out,  however,  that  it  possible  for  a  waveguide  with  a  different  thickness  and/or  different  indices  of 
refraction  for  the  waveguide  and  dielectrics  to  satisfy  the  phase-matching  condition  even  though 
the  infra-red  reading  beam  and  the  dominant  mode  in  which  the  generated  green  propagates  do  not 
correspond  to  the  modes  of  the  seed  infra-red  and  green  beams.  Of  course,  in  this  case,  the  overlap 
between  the  x ^  grating,  the  fundamental  reading  wave,  and  the  generated  second-harmonic  wave 
is  not  as  favorable. 

A  result  that  is  different  than  conventional  SHG  is  that  the  combination  of  different  modes  of 
infra-red  and  green  seed  gave  the  same  order  of  magnitude  of  SHG  power  as  long  as  the  infra-red 
used  during  the  reading  process  is  in  the  same  TM  mode  as  the  writing  mode  even  though  the 
higher  order  modes  are  not  as  confined  to  the  waveguide,  where  the  SHG  occurs,  as  the  lowest 
order  mode.  For  example,  when  we  combined  a  TM\  for  the  seed  infra-red  with  a  TM\  mode 
for  the  seed  green  and  read  that  with  a  TAfi  infra-red,  the  same  order  of  magnitude  of  SHG  was 
obtained  as  that  of  the  case  when  all  beams  are  in  TMo>  This  is  in  agreement  with  experimental 
results  [17].  Figure  (2)  shows  the  growth  of  the  SHG  power  in  the  TM\  mode,  the  dominant  mode 
for  this  case,  as  a  function  of  distance  along  the  direction  of  propagation.  The  SHG  power  for  this 
case  is  about  0.4  of  that  when  all  beams  are  in  the  lowest  order  mode. 

Table  (1)  compares  the  SHG  for  various  waveguiding  modes  normalized  to  the  case  in  which 
both  the  infra-red  and  green  are  launched  in  the  TMo  mode.  Note  that,  for  an  infra-red  beam 
launched  in  a  given  mode,  the  maximum  SHG  obtained  is  when  the  green  seed  is  launched  in  the 
same  mode. 

The  results  obtained  in  the  table  above  can  be  explained  by  noting  that  for  this  particular 
process  of  SHG,  and  with  the  reading  and  writing  infra-red  beams  in  the  same  mode,  the  overlap 
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integral  has  the  same  form  as  /3 (z)  given  in  Eq.  (21).  The  first  term  in  the  integrand  is  the  dominant 
term  and  it  involves  a  multiple  of  even  powers  of  the  infra-red  and  green  modes.  Therefore,  unlike 
other  SHG  processes,  there  is  no  cancellation  in  the  overlap  integrals  when  different  modes  are 
used,  and  thus  these  integrals  are  comparable  in  value.  Still  the  best  overlap  occurs  when  both  the 
infra-red  and  green  are  in  the  same  mode  as  the  modes  will  almost  match  each  other.  Figures  (3a) 
and  (3b)  show  plots,  properly  normalized,  of  the  integrand  vs.  depth  in  waveguide  for  two  different 
combinations  of  modes,  with  y  =  -1/zm  representing  the  waveguide-silica  substrate  interface,  and 
y  =  1  jim  representing  the  waveguide-air  interface.  The  area  underneath  each  of  the  plots  represents 
the  overlap  in  y. 

We  also  note  that  as  long  as  the  infra-red  power  used  during  the  reading  process  is  the  same 
as  that  during  the  writing  process,  the  effects  of  SPM  and  XPM  are  insignificant.  If  the  intensities 
are  considerably  different,  however,  the  efficiency  is  decreased.  This  can  be  attributed  to  that 
different  intensities  produce  different  SPM  and  XPM  and  the  net  effect  in  the  phase  mismatch  is 
nonvanishing  and  builds  up  along  the  direction  of  propagation.  To  illustrate  this,  SPM  and  XPM 
were  turned  off  during  the  reading  process  and  the  results  were  compared  to  the  case  when  both 
SPM  and  XPM  are  included.  Thus  the  case  when  the  grating  is  written  under  high  intensity  but 
read  with  a  much  weaker  beam  is  modeled.  Figure  (4)  shows  a  plot  of  SHG  power  vs.  distance 
along  the  direction  of  propagation  for  both  cases. 

4  Conclusion 

We  have  formulated  the  problem  of  SHG  in  germanium-doped  planar  waveguides  that  possess 
an  effective  second-order  susceptibility  induced  by  seeding.  This  formalism  is  based  on  a  normal 
mode  analysis  that  includes  the  effects  of  self-  and  cross-phase  modulation  of  the  fundamental  and 
its  second-harmonic.  The  resultant  equations  describing  the  amplitude  evolution  of  the  second- 
harmonic  were  solved  for  various  combinations  of  waveguiding  modes.  Unlike  the  traditional  means 
of  producing  SHG,  we  note  that  considerable  power  is  generated  in  the  higher  order  modes.  The 
conversion  efficiency  is  limited  by  diffraction. 
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6  List  of  Captions  for  Figures 

Figure  (1)  SHG  Power  in  TM0  vs.  distance  along  z  for  a  grating  written  with  both  fundamental  and 
second-harmonic  in  TMq  mode  and  read  with  fundamental  in  TMq  mode  a)  linear  losses  included 
b)  zero  linear  losses. 

Figure  (2)  SHG  Power  in  TM\  vs.  distance  along  z  for  a  grating  written  with  both  fundamental 
and  second-harmonic  in  TM\  mode  and  read  with  fundamental  in  TM\  mode. 

Figure  (3)  Overlap  vs.  waveguide  depth  for  a)  infra-red  and  green  in  TMq  b)  infra-red  in  TMq  and 
green  in  TM2  ■ 

Figure  (4):  SHG  power  vs.  distance  along  z  showing  effect  of  SPM  and  XPM. 


Fundamental  Mode 

Second  Harmonic 

Normalized  Power 

TMO 

TMO 

1.00 

TMO 

TM1 

0.30 

TMO 

TM2 

0.36 

TM1 

TMO 

0.13 

TM1 

TM1 

0.40 

TM1 

TM2 

0.21 

Table  1:  SHG  power  for  some  of  the  different  infra-red  and  green  mode  combinations  normalized 
the  case  in  which  all  beams  are  in  the  lowest  order  mode. 
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The  level  of  second  harmonic  generation  (SHG)  in  poled  fused  silica  was  altered  by  pre-annealing  in 
aqueous  and  dry  nitrogen  atmospheres.  The  effect  of  water  and  heat  was  found  to  be  dependent  upon  the 
lot  of  material.  The  effect  on  the  SHG  of  either  a  wet  or  dry  pre-anneal  was  sometimes  found  to  be 
reversible  by  performing  the  opposite  pre-annealing  process.  For  both  wet  and  dry  pre-annealed  samples, 
the  SHG  was  found  to  have  both  near-surface  (<25pm)  and  bulk  components  («1.5mm).  Electron 
paramagnetic  resonance  (EPR)  data  reveals  the  presence  of  Si  and  Ge  E'  defect  sites  and  shows  that  the 
relative  concentrations  of  these  sites  are  altered  by  the  pre-annealing  processes.  EPR,  along  with  SHG  data, 
provide  evidence  for  a  second-order  "bond  effect"  along  with  a  third-order  "hole  filling  effect"  as  sources  of 
the  nonlinearity. 
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1.  INTRODUCTION 


Second  harmonic  generation  (SHG)  in  poled  fused  silica,  which  was  first  reported  by  Myers,  Mukherjee 
and  Brueck(1)  in  1991,  has  increased  in  importance  because  of  a  multitude  of  potential  device  applications. 
Such  applications  include  electro-optic  switching,  frequency  doubling  crystals,  parametric  amplification, 
parametric  oscillation,  and  linear  electro-optic  modulation/frequency  conversion  which  can  be 
monolithically  integrated  into  optical  fibers  and  planar  integrated  circuit  geometry’s.  Because  of  these 
potential  applications,  much  work  has  been  done  to  understand  the  phenomenon  of  SHG  in  poled  fused 
silica. 

In  Myers  original  paper  on  the  SHG  in  poled  fused  silica(I),  an  optically  active  layer  of  thickness  5-10  pm 
on  the  positively  biased  face  of  the  sample  was  reported.  In  addition,  the  nonlinear  optical  coefficient,  d33, 
was  found  to  be  on  the  order  of  1  pm/V  in  poled  natural  fused  quartz.  But  for  higher  purity  UV  grade 
synthetic  fused  silica  having  only  10  percent  as  many  impurities  as  natural  fused  quartz,  the  magnitude  of 
d33  was  found  to  be  an  order  of  magnitude  lower.  Based  on  this  data,  Myers  felt  that  the  SHG  originated 
from  a  space  charge  region  created  by  migrating  impurity  ions  and  electrons.  Mukherjee,  et  al.(2),  further 
proposed  that  the  electrons  were  trapped  in  non-bridging  oxygen  hole  centers  associated  with  silicon  which 
aligned  upon  poling. 

In  a  later  paper,  Nasu  et  al.(3),  reported  a  direct  dependence  of  the  SHG  on  the  hydroxyl  content  in  Si02 
glasses  prepared  by  different  methods.  Nasu  also  reported  second-order  nonlinear  susceptibilities  on  the 
order  of  .05  pm/V  in  addition  to  well-developed  fringes  in  Maker-SHG  scans  of  poled  sol-gel  derived  fused 
silica.  Based  on  the  number  and  angular  positions  of  the  fringes  along  with  the  indices  of  refraction  for 
1.06  jitm  and  532  nm  radiation  in  fused  silica,  the  interaction  layer  was  found  to  be  the  same  thickness  as  the 
sample  thickness.  Finally,  Kazansky,  et  al.(4)  reported  the  erasure  of  SHG  in  fused  silica  by  electron-beam 
irradiation  and  speculated  that  this  effect  may  be  the  result  of  electrons  neutralizing  positive  charges 
resulting  in  the  elimination  of  the  space  charge  region.  Kazansky  also  found  the  optically  active  region  to 
reside  from  12  to  19  pms  below  the  positively  biased  surface,  in  contradiction  to  that  reported  by  Myers0 1 
Because  the  poling  process  involves  a  charging  phenomenon,  work  by  Lai,  et  al.(5),  Nicollian,  et  al.(6,7),  and 
Olthuis,  et  al.(8)  on  the  role  of  water  in  the  charging  of  fused  silica  is  also  relevant.  By  dry  annealing  Si02  in 
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N2  for  1  hour  at  1000°C,  Lai  found  that  the  number  of  electron  traps  were  reduced  nearly  fourfold.  Lai 
speculated  that  these  traps  were  associated  with  water-related  centers.  Along  the  same  lines,  Nicollian 
reported  that  the  ability  to  charge  Si02  by  avalanche  injection  was  dependent  upon  the  presence  of  water  in 
the  material.  Nicollian  also  found  a  one-to-one  relationship  between  the  level  of  electron  charging  and  the 
amount  of  hydrogen  evolving  from  the  material.  Based  on  this  data,  Nicollian  speculated  the  electron 
trapping  center  to  be  the  non-bridging  oxygen  hole  center  (NBOHC).  Finally,  Olthuis,  also  having 
observed  hydrogen  evolution  upon  charging  of  Si02,  proposed  the  following  mechanism. 


■  Si-O-Si  i  +  H2Oo2i  Si-OH 

(1) 

=  Si-OH  +  H20  <=>  ■  Si-O'  +  H30+ 

(2) 

=  Si-O'  +  H30+  +  e~  <=>  =  Si-O'  +  H20  +  H. 

(3) 

Currently,  the  mechanism  for  second  harmonic  generation  in  poled  fused  silica  and  its  relationship  to 
defect  sites  is  not  well  understood.  We  will  show  that  the  level  of  second  harmonic  generation  along  with 
the  relative  concentrations  of  Si  and  Ge  E'  defect  sites  can  be  altered  by  wet  and  dry  pre-annealing  fused 
silica  over  a  range  of  temperatures.  In  addition,  we  will  present  a  model  for  the  SHG  in  poled  fused  silica 
and  show  experimental  evidence  in  support  of  this  model. 

2.  EXPERIMENTAL  METHODS 

Manipulation  of  the  level  of  SHG  and  defect  sites  in  fused  silica  was  accomplished  by  pre-annealing  the 
material  in  an  aqueous  or  dry  nitrogen  environment  at  temperatures  up  to  900°C.  Infrared  grade 
commercial  fused  silica  having  low  water  content  was  used  for  all  of  the  experiments  described  in  this 
paper.  To  determine  the  effect  on  the  SHG  of  pre-annealing  under  different  conditions,  the  nonlinear 
optical  coefficient,  d33,  of  the  pre-annealed  material  after  poling  was  compared  with  d33  of  the  poled 
untouched,  i.e.,  as  received ,  material.  To  correlate  the  level  of  SHG  with  specific  paramagnetic  defect  sites 
in  fused  silica,  EPR  was  carried  out  on  the  pre-annealed  materials. 
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A  total  of  thirty  six  as  received  samples,  each  1.59  mm  thick  and  25  mm  in  diameter,  was  used  for  the 
pre-anneals.  Three  as  received  samples,  one  each  from  lots  1,  2  and  3  for  the  wet  pre-anneals  and  from  lots 
2,  3  and  4  for  the  dry  pre-anneals,  were  used  at  each  individual  pre-annealing  temperature.  Each  sample 
underwent  only  one  pre-annealing  process,  either  wet  or  dry,  and  was  never  reused  for  a  second  pre-anneal. 
A  “lot”  of  material  was  considered  to  be  a  group  of  samples  received  from  the  manufacturer  at  one  time. 
Such  groups  of  samples  were  generally  purchased  months  apart. 

The  as  received  samples  were  initially  parallel-plate  poled  to  generate  the  nonlinear  optical  effect.  The 
poling  process  involved  heating  a  sample  in  an  oven  from  room  temperature  to  250°C  with  7000  V  applied 
across  the  plates.  This  typically  took  1/2  hour.  Following  this,  the  sample  was  maintained  at  250°C  for  1 
hour  after  which  it  was  cooled  to  room  temperature  with  the  voltage  still  applied. 

The  baseline  SHG  susceptibility  of  the  poled  sample  was  obtained  using  the  Maker-Fringe  technique(9). 
The  experimental  setup  for  this  is  shown  in  Figure  1 .  1.06  pm  radiation  from  a  Q-switched  Nd/YAG  laser, 
which  was  passed  through  a  RG695  filter  to  remove  visible  radiation  as  well  as  attenuation  filters  to  reduce 
the  beam  power  to  approximately  10  mW,  was  focused  on  each  sample  individually  to  maximize  the  second 
harmonic  signal.  Upon  irradiation  of  the  sample  with  1.06  pm  radiation,  the  resultant  frequency  doubled 
radiation  and  residual  1.06  pm  radiation  passed  through  a  BG18  filter  and  a  530  nm  bandpass  filter  for 
removal  of  the  1.06  pm  component  before  detection  by  the  photomultiplier  tube.  Typically,  3000  Q- 
switched  pulses  were  averaged  by  the  boxcar  averager  and  recorded  in  the  form  of  a  dc  voltage  by  a 
computer. 

To  generate  the  Maker-Fringe  pattern,  the  sample  was  rotated  about  an  axis  perpendicular  to  the  p- 
polarized  laser  beam  through  the  angular  range  -80°  <  0  <  +80°  in  one  degree  steps.  The  resulting  Maker- 
Fringe  pattern,  Figures  2  and  3,  is  the  result  of  a  phase  mismatch,  Ak,  between  the  fundamental  and  second 
harmonic  waves,  equation  4: 

Ak  ss  (47t/X)(n£Ocos0O3'  -  n2oicos02(o')  (4) 

where  nw  and  n2(0  are  the  indices  of  refraction  and  0</  and  02o)'  are  the  angles  of  refraction  for  the 
fundamental  and  second  harmonic  waves,  respectively,  and  X  is  the  wavelength  of  the  fundamental  beam. 
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Theoretically,  the  power  in  the  second  harmonic  is  of  the  form:  P2(0  =  f(0)d2sin2T(L,0)  where  0  is  the 
rotation  angle  as  indicated  in  Figure  1,  d  is  the  nonlinear  optical  coefficient  and  L  is  the  interaction  length. <9) 
Generally,  f(0)d2  forms  the  envelope  describing  the  low-frequency  structure  associated  with  the  Maker- 
Fringe  curve  while  sin^fL,©)  creates  the  oscillations.  In  addition,  the  envelope,  f(0)d2,  is  independent  of 
the  phase  mismatch,  equation  4.  Because  a  sample  of  poled  fused  silica  exhibits  uniaxial  or  C„,v  symmetry, 
the  nonlinear  polarization,  Pnl,  is  of  the  following  form: 


Pnl  = 


/0 


(d3, 


0 

0 

d3i 


0  0 
0  d3j 
d33  0 


d3i 

0 

0 


O' 

0 

0 


>1 


^Ei2 

e22 

e32 

2E2E3 

2E]E3 

^e,e2 


(5) 


Since  p-polarized  radiation  was  used  (i.e.  Ei^Ecos©^',  E2=  0,  E3=Esin0a)/),  and  d31  =  1/3  d33  for  fused 
silica(10),  a  nonlinear  polarization  proportional  to  d33  was  obtained: 


PNL  =  EZd 


33 


r 


-\ 

2/3sin0(O'cos0tl/ 

0 

1/3cos20co'  +  sin2©^'  ^ 


(6) 


Because  P2cD  is  proportional  to  IPnlI2,  the  nonlinear  optical  coefficient,  d33  was  evaluated.  To  get  an 
absolute  value  for  d33,  a  y-cut  quartz  standard  was  used  as  a  reference.  d33  was  then  computed  from  the 
following  equation 


where  the  subscripts  s  and  y  represent  the  sample  and  y-cut  quartz  standard,  respectively,  and  dn>y 
represents  the  nonlinear  optical  coefficient  associated  with  y-cut  quartz.  Since  an  interaction  length  could 
not  be  obtained  from  the  Maker-Fringe  scans  in  the  absence  of  fringes,  calculations  of  d33  were  based  on  an 
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interaction  length  equivalent  to  the  coherence  length.  As  a  result,  the  values  of  d33  reported  represent  the 
lower  bound.  Finally,  standard  deviations  associated  with  variations  in  the  sample  surfaces  for  the  values  of 
d33  reported  were  estimated  from  measurements  to  be  on  the  order  of  ±10  percent.  Since  a  lower  bound  of 
d33  is  reported,  the  actual  errors  associated  with  the  values  of  d33  reported  are  almost  certainly  greater  than  ± 
10  percent 

After  the  SHG  measurement,  the  sample  was  placed  between  grounded  parallel  plates  for  depoling  to 
remove  the  nonlinear  effect.  As  the  plates  were  brought  into  close  proximity  to  the  sample,  an  image  charge 
developed  on  the  plates  in  proportion  to  the  level  of  polarization  in  the  sample.  As  the  sample  was  heated 
to  250°C  and  the  depoling  process  proceeded,  the  polarization  of  the  sample  was  eliminated  and  the  amount 
of  image  charge  on  the  plates  decreased  in  proportion.  This  loss  of  image  charge  was  recorded  as  a 
function  of  time. 

After  the  sample  was  depoled,  it  was  either  wet  or  dry  pre-annealed  in  a  one-zone  tube  furnace.  The 
initial  stage  of  the  pre-annealing  process  involved  centering  the  sample  in  the  tube  furnace  and  allowing 
nitrogen  to  purge  the  system  of  atmospheric  gases  for  20  minutes.  The  temperature  was  then  increased  to 
the  specified  pre-annealing  temperature  between  400°C  to  900°C  for  1.5  hours.  Upon  completion  of  the 
pre-annealing  process,  the  sample  was  allowed  to  cool,  within  the  tube  furnace,  to  room  temperature.  The 
only  difference  between  the  wet  and  dry  pre-annealing  processes  was  that  dry  pre-annealing  involved  direct 
entry  of  nitrogen  into  the  tube  furnace  whereas  wet  pre-annealing  involved  bubbling  nitrogen  through  water 
at  approximately  75°C  before  entry  into  the  tube  furnace. 

Once  the  sample  had  been  pre-annealed,  it  was  repoled  and  the  SHG  was  recorded  again  using  a  Maker- 
SHG  scan.  To  determine  the  extent  of  the  polarization  within  the  sample,  the  sample  was  again  depoled. 

FTIR  spectra  of  both  the  as  received  and  pre-annealed  samples  were  obtained  with  a  Nicolet  51  OP  FT-IR 
spectrometer  in  reflectance  mode  between  400  cm"1  to  4800  cm'1.  1064  scans  were  averaged  in  order  to 
reduce  the  noise  level  and  allow  quantitative  comparison  of  spectra  from  different  pre-annealing 
temperatures.  EPR  spectra  were  obtained  at  room  temperature  on  a  Bruker  ESP-300  X-band  spectrometer 
using  the  out-of-phase  high  power  second  harmonic  mode  technique.00  The  EPR  parameters  utilized  are 
listed  in  Table  1  below. 
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3.  RESULTS 


A.  EFFECT  OF  PRE-ANNEALING  ON  THE  LEVEL  OF  SHG 

The  ability  to  alter  the  level  of  d33  in  a  sample  of  fused  silica  by  wet  or  dry  pre-annealing  was  found  to  be 
highly  dependent  on  the  particular  lot  of  material.  Pre-annealing  data  for  four  lots  of  fused  silica  are  shown 
in  Figures  4  and  5.  Samples  from  lots  2  and  3  have  been  both  wet  and  dry  pre-annealed  whereas  samples 
from  lot  1  were  only  wet  pre-annealed  and  samples  from  lot  4  were  only  dry  pre-annealed.  From  variation 
of  d33  across  the  sample  surface  and  variation  in  d33  among  samples  from  the  same  lot  of  material  (typically 
10  percent),  errors  associated  the  percentage  changes  in  d33  shown  in  Figures  4  and  5  are  estimated  to  be  on 
the  order  of  ±12  to  ±16  percent.  These  errors  are  shown  by  a  few  representative  error  bars  in  the  figures.  It 
must  be  pointed  out  that  although  the  variation  of  d33  among  samples  from  the  same  lot  is  approximately  10 
percent,  the  variation  between  samples  from  different  lots  can  be  much  more. 

For  wet  pre-annealing,  shown  in  Figure  4,  the  nonlinearity  was  decreased  for  all  lots  of  material  studied 
over  the  temperature  range  of  pre-anneal  investigated  with  the  possible  exception  of  lot  3  in  the  temperature 
range  around  600°C.  In  addition,  the  decrease  in  d33  was  greater  the  further  the  pre-annealing  occurred 
from  the  500-600°C  temperature  range.  Unlike  lots  1  and  2  which  retained  a  sizable  d33  when  pre-annealed 
at  temperatures  greater  than  600°C,  the  SHG  associated  with  lot  3  was  virtually  quenched  at  higher  pre¬ 
annealing  temperatures. 

While  the  trends  for  wet  pre-annealing  were  somewhat  similar  for  all  lots  of  fused  silica  studied,  the 
trends  for  dry  pre-annealing  were  not.  Samples  from  lots  2  and  4,  when  dry  pre-annealed,  shown  in  Figure 
5,  exhibited  linear  behavior  in  the  percentage  change  in  d33  as  a  function  of  temperature.  When  dry  pre- 
annealed,  the  behavior  of  samples  from  lot  2  was  quite  distinct  from  the  behavior  when  wet  pre-annealed. 

On  the  other  hand,  the  behavior  of  samples  from  lot  3  when  dry  pre-annealed  was  very  similar  to  the 
behavior  when  wet  pre-annealed.  Apparently  the  presence  of  water  was  able  to  significantly  affect  the 
defect  sites  generating  second  harmonic  in  lot  2  whereas  in  lot  3,  it  was  not.  As  is  apparent  from  the 
figures,  depending  on  the  temperature  range  studied  and  the  lot  of  fused  silica,  an  increase  or  decrease  in  d33 
may  occur  by  both  wet  and  dry  pre-annealing.  Finally,  the  more  extreme  changes  in  d33  brought  about  by 
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pre-annealing  are  typically  much  greater  than  the  individual  deviations  in  d33  observed  between  as  received 
samples  from  different  lots. 

The  increase  or  decrease  of  d33  brought  about  by  either  wet  or  dry  pre-annealing  has  been  found  to  be 
reversible  to  a  certain  extent  upon  performance  of  the  opposite  pre-annealing  process.  As  an  example,  a 
sample  of  fused  silica  having  a  d33  of  .34  pm /V  after  the  initial  poling,  obtained  a  d33  of  .16  pm /V  upon  dry 
pre-annealing  in  N2  for  1  hour  at  550°C.  Upon  subsequent  wet  pre-annealing  under  the  same  conditions,  a 
d33  of  .35  pm/V,  which  was  nearly  equal  to  the  initial  as  received  level,  was  obtained.  In  summary,  the  d33 
after  being  decreased  over  50  percent  by  a  dry  pre-anneal,  was  reestablished  to  its  initial  level  upon  wet  pre¬ 
annealing.  The  fact  that  the  presence  of  water  was  the  only  difference  between  the  two  pre-annealing 
processes  supports  the  idea  that  water  is  a  component  in  an  equilibrium  process  between  a  defect  site  in  its 
optically  active  and  inactive  configurations.  The  same  reversibility  was  also  seen  when  d33  was  initially 
increased  by  wet  pre-annealing  and  then  decreased  by  a  subsequent  dry  pre-anneal.  As  an  example,  a 
sample  of  fused  silica  having  an  initial  d33  of  .32  pm/V,  obtained  a  d33  of  .39  pm/V  when  wet  pre-annealed 
for  1  hour  at  550°C  ,  a  21  percent  increase.  Upon  subsequent  dry  pre-annealing  under  the  same  conditions, 
the  d33  was  decreased  to  .12  pm/V,  a  decrease  of  69  percent. 

Although  d33  was  fully  reversible  for  some  lots  of  fused  silica ,  for  other  lots  of  material,  d33  was  only 
partially  reversible.  As  an  example,  a  sample  of  fused  silica  having  an  initial  d33  of  .32  pm/V,  obtained  a 
d33  of  .19  pm/V  when  wet  pre-annealed,  a  40  percent  decrease.  When  the  sample  was  subsequently  dry  pre¬ 
annealed,  the  d33  was  only  increased  to  .20  pm/V,  a  6  percent  increase.  Apparently  an  irreversible  change 
of  some  sort  occurred.  According  to  Lee(12),  fused  silica  contains  both  permanent  hydroxide  which  is 
formed  in  the  material  during  manufacture  and  metastable  hydroxide  which  is  formed  upon  reaction  of 
fused  silica  with  water  vapor.  According  to  both  Lee(12)  and  Shelby(13),  metastable  hydroxide  which  has 
been  added  to  a  material  may  not  be  totally  removable.  This  argument  coupled  with  evidence  to  be 
presented  later  on  the  role  of  hydroxide  in  the  generation  of  second  harmonic  may  explain  the 
nonreversibility  of  the  pre-annealing  processes. 

Upon  poling  a  sample  that  had  been  wet  or  dry  pre-annealed,  multiple  fringes  typically  emerge  on  the 
Maker-SHG  scan  as  shown  in  Figures  2  and  3  for  wet  pre-annealed  samples  from  lots  2  and  3.  These 
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fringes  are  indicative  of  an  optically  active  region  throughout  the  entire  sample  bulk.  Such  fringes  do  not 
usually  appear  in  poled  as  received  samples,  as  shown  in  Figure  6.  The  two-lobed  pattern  shown  in  Figure 
6  is  characteristic  of  an  interaction  layer  thickness  less  than  the  coherence  length,  i.e.,  a  near-surface  effect. 
The  largest  bulk  effect  relative  to  the  near-surface  effect  has  been  found  to  occur  when  the  SHG  had  been 
significantly  altered.  Upon  dehydration,  the  emergence  of  the  bulk  effect  relative  to  the  near-surface  effect 
is  thought  to  occur  from  removal  of  the  strong  near-surface  layer  by  the  pre-annealing  process.  The 
strength  of  the  SHG  associated  with  the  bulk  effect  generated  in  this  way  is  very  low.  This  is  seen  in 
Figures  3(d-f)  for  samples  from  lot  3  that  had  been  wet  pre-annealed  at  temperatures  between  700°C  and 
900°C.  When  the  SHG  remained  strong  and  a  bulk  effect  emerged,  hydration  throughout  the  sample  bulk  is 
thought  to  be  the  leading  factor.  It  has  been  shown  by  Holmberg,  et  al.(14)  that,  even  though  the  near-surface 
component  of  the  hydration  profile  is  an  order  of  magnitude  greater  than  that  of  the  bulk,  wet  pre-annealing 
for  longer  times  and  at  higher  temperatures  will  result  in  a  growth  of  the  bulk  component  relative  to  the 
near-surface  component.  This  is  evident  in  Figures  2(c-f)  for  samples  from  lot  2  that  had  been  wet  pre¬ 
annealed  between  600°C  and  900°C. 

To  determine  the  thickness  of  the  optically  active  region,  the  experimental  Maker-SHG  pattern,  shown  in 
Figure  3f,  was  fit  theoretically,  Figure  7.  Shown  in  Figure  7(a)  is  a  theoretical  Maker-SHG  pattern, 
generated  using  an  index  of  refraction  value  of  1.46071  for  532  nm  radiation  and  an  index  of  refraction 
value  of  1.44963  for  1.06  (im  radiation.  An  interaction  layer  thickness  of  1.501  mm,  equivalent  to  the 
sample  thickness  was  used.  Upon  comparing  Figure  7(a)  with  Figure  7(b),  it  is  apparent  that  the  angular 
position  of  the  fringes  in  the  experimental  curve  coincide  with  those  in  the  theoretical  curve.  This 
coincidence  of  angles  implies  that  the  SHG  effect  extends  throughout  the  entire  sample  bulk.  By  examining 
the  fringing  patterns  in  Figures  2  and  3,  it  is  apparent  that  bulk  and  near-surface  components  of  the 
nonlinearity  are  simultaneously  present.  The  relative  component  of  the  SHG  from  the  near-surface  and  bulk 
may  vary  from  an  almost  pure  near-surface  effect,  Figure  2a,  to  an  almost  pure  bulk  effect,  Figure  3f. 
Finally,  one  conclusion  that  can  be  drawn  is  that  at  least  one  microscopic  cause  of  the  SHG  is  not  inherent 
to  the  near-surface. 
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B.  EFFECT  OF  PRE-ANNEALING  ON  THE  LEVEL  OF  DEFECT  SITES 

High  power  second  harmonic  mode  electron  paramagnetic  resonance(11)  was  performed  on  wet  and  dry 
pre-annealed  samples  from  lots  2  and  3.  All  samples  were  depoled  except  the  samples  from  lot  2  that  had 
been  wet  pre-annealed.  For  all  lots  of  fused  silica  studied,  both  Si  and  Ge  E'  centers  were  observed.  Figure 
8.  Note  the  presence  of  two  peaks  in  very  close  proximity  representing  two  distinct  paramagnetic  electron 
centers,  each  without  hyperfine  interaction.  The  peak  located  at  a  magnetic  field  of  3492.20  G  with  an 
associated  electron  Zeeman  factor,  g,  of  1.9995  has  been  attributed  to  a  germanium  E'  center  and  that 
located  at  a  magnetic  field  of  3487.91  G  with  an  associated  g  value  of  2.001  has  been  attributed  to  the  Si  E' 
center.(15)  Because  measurements  were  made  in  second  harmonic  mode,  no  absolute  germanium  density 
could  be  computed.  From  the  very  weak  EPR  signals  obtained,  concentrations  of  Si  and  Ge  E’  defect  sites 
were  probably  less  than  1015/cm3.(16) 

Upon  comparing  the  level  of  defect  sites  in  polarized  (poled)  and  nonpolarized  (depoled)  fused  silica,  no 
significant  differences  were  noted.  This  was  to  be  expected  since  the  SHG  was  largely  confined  to  the  near¬ 
surface  while  EPR  looked  primarily  at  the  bulk. 

Shown  in  Figure  9  are  the  relative  concentrations  of  the  Si  and  Ge  E'  centers  along  with  d33  and  FTIR 
data  for  wet  pre-annealed  samples  from  lot  2.  The  behavior  of  d33  can  be  explained  in  two  ways,  one  of 
which  points  to  metastable  hydroxyl  defects  and  the  other  to  permanent  hydroxyl  defects  as  sources  of  the 
nonlinearity.  First  since  there  appears  to  be  an  inverse  relationship  between  the  level  of  Si  and  Ge  E’  defect 
sites  and  d33,  it  is  speculated  that  the  optically  active  centers  could  be  the  precursors  to  the  Si  and  Ge  E1 
defect  sites.  Since  water  appears  to  be  a  factor  causing  a  decrease  in  d33  on  either  side  of  a  600°C  pre¬ 
annealing  temperature  for  this  lot  of  fused  silica,  as  shown  in  Figures  4  and  5,  one  may  speculate  that  water 
is  involved  in  some  sort  of  equilibrium  reaction  between  an  optically  active  defect  site  involving  metastable 
hydroxyl  and  the  E’  center.  In  the  literature,  Devine(17)  speculates  that  the  Si  and  Ge  E'  centers  arise  from 
their  precursors  via  the  following  reactions: 

s  Si+  .Si  =  +  H20  <=>  =  Si+  HO-Si  s  +  H.  (8) 

=  Si+  .Ge  s  +  H20  <=>  =  Si+  HO-Ge  s  +  H.  (9) 
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where  =  Si+  .Si  s  represents  a  E'  center.  [Note:  the  symbol  s  represents  the  bonding  of  Si  or  Ge  to  three 

other  atoms  in  the  lattice.]  Since  the  s  Si-OH  and  =  Ge-OH  defects  are  the  precursors  to  the  non-bridging 
oxygen  hole  centers,  these  reactions  would  be  consistent  with  the  speculation  by  Mukherjee(2)  concerning 
the  role  of  non-bridging  oxygen  hole  centers  in  the  generation  of  SHG  and  the  reactions  (1-3)  presented  by 
01thuis(8)  for  the  charging  of  fused  silica.  Furthermore,  FTIR  spectra  of  wet  pre-annealed  samples  from  lot 
2,  as  shown  in  Figure  9,  support  the  above  idea  since  the  3649  cm'1  peak,  which  corresponds  to  vibrations 
associated  with  internal  Si-OH  or  pairs  of  Si-OH  mutually  linked  by  a  hydrogen  bond^8'20*,  is  a  maximum  at 
a  pre-annealing  temperature  of  700°C  when  the  Si  E’  centers  are  at  a  minimum.  The  fact  that  the  percentage 
change  in  the  FTIR  peak  was  less  than  the  corresponding  percentage  change  in  d33  may  be  explained  by  the 
fact  that  the  FT  IR  peak  was  composed  of  both  metastable  and  permanent  hydroxyl. 

Secondly,  the  behavior  of  d33  may  be  explained  by  the  presence  of  permanent  hydroxyl  and  water  in 
fused  silica.  According  to  Lee(12),  at  pre-annealing  temperatures  less  than  600°C,  water  diffuses  into  the 
fused  silica  lattice  and  remains  in  its  molecular  form.  It  is  speculated  that  the  observed  decrease  of  d33  in 
the  400°C  to  600°C  pre-annealing  temperature  range  results  from  the  alteration  of  the  permanent  hydroxyl 
defect  sites,  ^  Si-OH  and  =  Ge-OH,  through  hydrogen  bonding  with  water: 


(10) 


It  is  thought  that  water  stabilizes  the  hydroxyl  groups,  making  them  immune  to  alteration  through  poling. 
This  would  explain  the  observed  decrease  in  optical  activity  in  this  temperature  range.  As  the  temperature 
is  increased  from  400°C  to  600°C,  there  is  less  molecular  diffusion  and  hydrogen  bonding  but  more 
reaction  of  water  with  the  lattice  to  form  metastable  hydroxyl,  which  would  account  for  the  observed 
increase  in  d33.  Finally  according  to  Lee(12),  at  pre-annealing  temperatures  greater  than  700°C,  desorption 
of  metastable  hydroxide  occurs  which  would  explain  the  small  observed  decrease  in  d33  at  pre-annealing 
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temperatures  in  excess  of  600°C.  Since  the  maximum  d33  observed  for  wet  pre-annealed  fused  silica  was, 
within  experimental  error,  equal  to  that  for  the  dry  pre-annealed  samples,  it  is  speculated  that  metastable 
hydroxyl  played  an  insignificant  role  in  the  generation  of  second  harmonic  in  this  material. 

For  samples  from  lot  2  that  had  been  dry  pre-annealed,  behavior  similar  to  that  observed  for  the 
corresponding  wet  pre-annealed  samples  from  lot  2  was  found  for  the  relative  concentrations  of  Si  and  Ge 
E'  defect  sites,  as  shown  in  Figure  10.  But  while  d33  decreased  on  either  side  of  a  600°C  pre-annealing 
temperature  for  wet  pre-annealed  samples  from  lot  2,  the  level  of  d33  remained  almost  unchanged 
throughout  the  pre-annealing  temperature  range  for  the  corresponding  dry  pre-annealed  samples. 

Obviously,  no  correlation  exists  between  d33  and  either  the  E'  centers  or  their  removable  metastable 
precursors  for  the  dry  pre-annealed  samples  from  lot  2  except  possibly  for  pre-annealing  temperatures 
greater  than  800°C.  Because  no  other  paramagnetic  defect  sites  were  found  in  this  material,  it  is  speculated 
that  permanent  hydroxide  or  a  bond  effect  was  the  contributor  to  the  SHG.  Upon  examination  of  the  image 
charge  data  in  Figure  10,  it  is  apparent  that  the  level  of  image  charge  on  the  parallel  plates  decreased  as  the 
pre-annealing  temperature  increased.  This  was  unusual  since  the  image  charge  usually  followed  the  trends 
associated  with  d33,  see  Figures  1 1  and  12.  Since  the  image  charge  is  related  to  the  magnitude  and 
distribution  of  the  polarization  within  the  sample,  either  the  geometric  distribution  of  the  SHG  was  changing 
and/or  the  type  of  defect  site  generating  the  SHG  was  evolving.  Upon  examination  of  Maker-Fringe  data,  it 
is  apparent  that  the  geometric  distribution  was  not  changing  as  the  pre-annealing  temperature  increased.  So 
the  most  likely  explanation  of  the  data  is  that  the  microscopic  source  of  the  SHG  was  evolving  as  the  pre¬ 
annealing  temperature  increased. 

Upon  examination  of  Figures  1 1  and  12,  it  is  apparent  that  samples  from  lot  3  reacted  in  a  similar  way  to 
both  wet  and  dry  pre-annealing.  In  addition,  for  both  wet  and  dry  pre-annealed  samples,  d33  was  quenched 
at  pre-annealing  temperatures  above  600°C.  From  Figures  1 1  and  12,  it  is  apparent  that  the  level  of  Ge  and 
Si  E'  defect  sites  remained  relatively  constant  for  pre-annealing  temperatures  up  to  700°C.  d33  was 
obviously  not  correlated  to  the  level  of  Si  and  Ge  E'  defect  sites  or  their  precursors.  In  addition,  the  relative 
concentration  of  the  defect  sites,  within  experimental  error,  were  identical  for  both  the  wet  and  dry  pre- 
annealed  samples  in  the  pre-annealing  temperature  range.  This  was  unlike  the  behavior  seen  in  lot  2  where 
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the  relative  concentration  of  E’  defect  sites  for  the  dry  pre-annealed  material  was,  at  a  minimum,  double 
those  in  the  wet  pre-annealed  material.  The  relatively  constant  concentration  of  E'  defect  sites  over  the  pre¬ 
annealing  temperature  range  in  addition  to  the  lack  of  sensitivity  of  the  d33  to  water  provides  evidence  that 
water  did  not  react  with  the  E  defect  sites  as  in  lot  2.  For  pre-annealing  temperatures  above  800°C,  thermal 
growth  of  the  Si  and  Ge  E  defect  sites  was  seen.  Since  d33  was  virtually  zero  at  these  pre-annealing 
temperatures,  the  thermally  generated  defect  sites  were  obviously  not  active  in  the  generation  of  the  second 
harmonic.  Based  on  the  fact  that  there  was  no  correlation  between  d33  and  the  paramagnetic  defect  sites  in 
addition  to  the  fact  that  d33  was  quenched  at  pre-annealing  temperatures  above  600°C,  a  “bond  effect”  is 
proposed  to  be  a  source  of  SHG  in  this  lot  of  fused  silica.  This  is  further  supported  by  the  fact  that 
permanent  hydroxide  is  virtually  nonremovable  at  pre-annealing  temperatures  less  than  1000°C(12),  in 
addition,  to  the  relatively  extensive  bulk  contribution  to  the  SHG  compared  to  lot  2  material,  Figures  2  and 
3.  It  is  apparent  though  by  examining  Figure  3  that  when  the  SHG  reached  its  maximum  in  the  500-600°C 
pre-annealing  temperature  range,  a  near-surface  effect  was  predominant.  This  could  be  the  result  of  water 
uptake  to  form  metastable  hydroxide.  Based  upon  the  literature07, 19'21) ,  it  is  speculated  that  the  peroxy 
bridge  structure  in  the  trans  configuration,  which  is  converted  to  the  dy-like  configuration  lacking  an 
inversion  center(22)  upon  poling,  equation  6,  could  be  the  optically  active  bond  structure  in  lot  3  material. 


trans 


* 


O -  O 


(11) 


cis 


According  to  Devine,  et  al.°9, 20),  peroxy-radical  defect  sites,  which  form  from  the  reaction  of  diffusing  02 


with  oxygen  deficient  sites  within  the  material,  pre-anneal  to  form  peroxy  bridges  in  the  temperature  range 
300-600°C.  This  is  consistent  and  could  be  another  explanation  for  the  growth  of  SHG  in  this  temperature 
range  besides  metastable  hydroxyl.  In  addition,  observations  by  Freund,  et  al.(21)  indicate  a  thermal 
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dissociation  of  the  peroxy  bridge  at  pre-annealing  temperatures  above  500°C  into  two  =  Si-O.  groups  which 
are  suggested  by  Devine  et  al.(17)  to  pre-anneal  back  into  the  lattice  in  the  form  =  Si-O-Si  =.  This  would 
support  the  sudden  decline  of  the  SHG  to  nearly  zero  between  500-600°C.  Finally,  Freund^  ^  also  reports 
that  the  peroxy  bond  is  expected  to  orientate  itself  when  an  external  stress  is  applied.  This  would  support 
speculation  that  the  trans  configuration  is  converted  to  its  optically  active  form,  the  cw-like  configuration, 
upon  poling. 

4.  MODEL 

Based  on  data  associated  with  d33,  the  image  charge,  EPR  and  FTIR,  it  is  felt  that  the  SHG  in  poled  fused 
silica  is  likely  coming  from  two  sources.  It  is  proposed  that  a  “bond”  effect  in  addition  to  a  “hole-filling” 
effect  plays  a  role  in  the  generation  of  second  harmonic.  For  samples  of  fused  silica  showing  no  correlation 
of  d33  with  paramagnetic  defect  sites,  in  addition  to  quenching  at  pre-annealing  temperatures  greater  than 
600°C,  a  2nd  order  “bond”  effect  is  thought  to  be  generating  the  second  harmonic.  The  optically  active 
defect  structure  is  felt  to  the  peroxy  bridge  defect  in  the  cw-like  configuration.  It  is  proposed  that  upon 
poling,  the  trans-configuration  having  zero  net  dipole  moment  undergoes  rotation  about  the  0-0  bond  to 
form  the  c/s-like  configuration  which  has  a  net  dipole  moment,  equation  (6).  Since  the  ds-like 
configuration  is  energetically  unstable,  it  is  felt  that  the  during  the  poling  process  the  electric  field  aligns  the 
dipoles  which  are  then  frozen  in  place  when  the  sample  is  cooled.  Upon  depoling,  the  ds-like  configuration 
most  likely  reverts  back  to  the  lower  energy  trans-configuration  since  no  electric  field  is  present  to  hold  the 
structure  in  place.  The  net  result  of  such  a  structure  would  be  a  second  order  SHG  effect. 

Finally,  for  lots  of  fused  silica  that  maintain  an  appreciable  d33  at  higher  pre-annealing  temperatures,  it  is 
speculated  that  hydroxyl  defect  sites  associated  with  Si  and  Ge  are  the  sources  of  the  SHG.  One 
explanation  for  the  behavior  of  d33  upon  wet  and  dry  pre-annealing  is  represented  by  an  equilibrium 
between  E’  and  metastable  hydroxyl  defect  sites,  equations  (4)  and  (5).  The  SHG  is  then  the  greatest  when 
the  number  of  metastable  hydroxyl  defect  sites  are  maximized.  The  trends  in  d33  can  also  be  explained  by 
the  stabilization  of  permanent  hydroxyl  defect  sites  at  lower  pre-annealing  temperatures  through  hydrogen 
bonding  with  molecular  water  and  the  desorption  of  metastable  hydroxyl  at  higher  temperature  pre-anneals. 
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As  a  result,  at  both  lower  and  higher  pre-annealing  temperatures,  the  number  of  hydroxyl  defect  sites 
available  for  poling  is  decreased.  Finally,  when  a  sample  of  fused  silica  containing  =  Si-OH  and  =  Ge-OH 
defects  is  then  poled,  it  is  speculated  that  molecular  hydrogen  is  released  and  the  resultant  non-bridging 
oxygen  hole  centers  (=  Si-O.  or  s  Ge-O.)  are  charged. 

5.  CONCLUSIONS 

The  level  of  SHG  and  concentration  of  defect  sites  in  fused  silica  can  be  manipulated  by  wet  or  dry  pre¬ 
annealing  in  nitrogen.  In  some  lots  of  material,  the  effects  of  wet  and  dry  pre-annealing  do  not  differ 
significantly  and  the  SHG  goes  virtually  to  zero  when  wet  or  dry  pre-annealed  above  500-600°C.  In  other 
lots  of  fused  silica,  the  effects  of  wet  and  dry  pre-annealing  are  quite  distinct  in  that  wet  pre-annealing 
decreases  the  SHG  significantly  on  either  side  of  600°C  while  dry  pre-annealing  has  virtually  no  effect  over 
the  range  of  pre-annealing  temperatures  from  400-900°C. 

It  has  also  been  found  that  the  effect  on  the  SHG  of  either  a  wet  or  dry  pre-anneal  can  be,  in  some  cases, 
fully  reversed  by  the  opposite  pre-annealing  process. 

After  either  wet  or  dry  pre-annealing,  a  weak  bulk  SHG  effect  has  been  found  to  emerge.  This 
demonstrates  that  the  source  of  the  SHG  is  not  limited  to  defects  inherent  only  to  the  sample  near-surface. 

EPR  analysis  has  shown  that  fused  silica  contains  Si  and  Ge  E*  centers.  d33,  image  charge,  EPR  and  FTIR 
data  reveals  the  presence  of  one  possible  second  order  “bond  effect”  in  addition  to  a  third  order  “hole¬ 
filling  effect”.  The  molecular  source  of  the  “bond  effect”  is  proposed  to  be  the  peroxy  bridge  structure 
which  is  converted  from  the  trans  configuration  to  the  a  cw-like  configuration  upon  poling.  The  third  order 
“hole-filling  effect”  is  speculated  to  originate  from  the  charging  of  the  =  Si-O.  and  s  Ge-O.  non-bridging 
oxygen  hole  centers,  which  are  derived  from  metastable  and  permanent  hydroxyl  defect  sites,  to  set  up  a 
depletion  region. 

Finally,  there  is  still  little  known  about  the  molecular  cause  of  SHG  in  fused  silica.  Future  work  will 
focus  on  designing  experiments  to  confirm  or  deny  the  role  of  hydroxide  defect  sites  and  peroxy  bridge 
structures  in  the  generation  of  second  harmonic.  In  addition,  the  energetics  associated  potentially  optically 
active  defect  sites  will  be  studied  with  semi-empirical  and  ab-initio  quantum  chemistry  calculations. 
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Figure  Captions 


Figure  1  -  Experimental  setup  of  Maker-SHG  scan.  Inset  shows  the  geometrical  arrangement  of  the 
p-polarized  input  1.06  pm  beam  with  respect  to  the  poled  sample. 

A  =  RG695  filter,  B  =  attenuation  filters,  C  =  BG18  filter,  D  =  530  nm  band  pass  filter. 

Figure  2  -  Relative  intensity  of  SHG  [a.u.]  vs  rotation  angle  [degrees]  for  poled  samples  from  lot  2  that 
had  been  wet  pre-annealed  at  a  series  of  temperatures:  a)  400°C,  b)500°C,  c)600°C,  d)700°C, 
e)800°C,  f)900°C.  [Note:  the  amplitudes  of  the  Maker-Fringe  scans  are  not  drawn  to  scale.  The 
amplitudes  of  (b)-(f)  relative  to  (a)  are  indicated  at  the  far  right  of  each  Maker-Fringe  scan.] 

Figure  3  -  Relative  intensity  of  SHG  [a.u.]  vs  rotation  angle  [degrees]  for  poled  samples  from  lot  3  that 
had  been  wet  pre-annealed  at  a  series  of  temperatures:  a)  400°C,  b)500°C,  c)600°C,  d)700°C, 
e)800°C,  f)900°C.  [Note:  the  amplitudes  of  the  Maker-Fringe  scans  are  not  drawn  to  scale.  The 
amplitudes  of  (b)-(f)  relative  to  (a)  are  indicated  at  the  far  right  of  each  Maker-Fringe  scan.] 

Figure  4  -  Percentage  change  in  d33  for  a  series  of  samples  from  lots  1-3  that  have  been  wet  pre-annealed 
over  the  temperature  range  400°C  to  900°C.  [Note:  samples  from  lot  1  have  only  been  wet 
pre-annealed  whereas  samples  from  lots  2  and  3  have  been  both  wet  and  dry  pre-annealed,  see 
Figure  2.] 

Figure  5  -  Percentage  change  in  d33  for  a  series  of  samples  from  lots  2-4  that  have  been  dry  pre-annealed 
over  the  temperature  range  400°C  to  900°C.  Note  that  samples  from  lot  4  have  only  been  dry 
pre-annealed  whereas  samples  from  lots  2  and  3  have  been  both  wet  and  dry  pre-annealed,  see 
Figure  1. 

Figure  6  -  Relative  intensity  of  SHG  [a.u.]  vs  rotation  angle  [degrees]  for  a  poled  as  received  sample. 

Figure  7  -  a.  Theoretical  Maker-SHG  fit  using  an  index  of  refraction  value  of  1.46071  for  532  nm  radiation 
and  an  index  of  refraction  value  of  1.44963  for  1.06  pm  radiation.  An  interaction  layer  thickness 
of  1.65  mm,  which  was  equivalent  to  the  sample  thickness,  was  also  used,  b.  Experimental  maker 
SHG  scan  of  a  sample  after  dry  pre-annealing 
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Figure  8  -  High  power  electron  paramagnetic  resonance  spectrum  of  a  sample  of  fused  silica.  The  peak 
located  at  a  magnetic  field  value  of  3485.23  G  has  been  attributed  to  the  Ge  E'  center  and  that 
located  at  3480.96  G  has  been  attributed  to  the  Si  E'  center. 

Figure  9  -  Samples  from  lot  2  that  have  been  wet  pre-annealed  over  a  series  of  temperatures  from  400°C  to 
900°C,  a.  d33  [pm/V]  vs  pre-annealing  temperature  [°C].,  b.  relative  intensity  of  3649  cm'1  peak 
[a.u.],  c.  Relative  concentration  of  Si  E’  centers  [a.u.]  vs  pre-annealing  temperature  [°C],  d. 
Relative  concentration  of  Ge  E'  centers  [a.u.]  vs  pre-annealing  temperature  [°C]. 

Figure  10  -  Samples  from  lot  2  that  have  been  dry  pre-annealed  over  a  series  of  temperatures  from  400°C  to 
900°C,  a.  d33  [pm/V]  vs  pre-annealing  temperature  [°C].,  b.  Image  charge  *[105]  [°C]  vs 
pre-annealing  temperature  [°C],  c.  Relative  concentration  of  Si  E'  centers  [a.u.]  vs  pre-annealing 
temperature  [°C],  d.  Relative  concentration  of  Ge  E'  centers  [a.u.]  vs  pre-annealing 
temperature  [°C]. 

Figure  1 1  -  Samples  from  lot  3  that  had  been  wet  pre-annealed  over  a  series  of  temperatures  from  400°C  to 
900°C,  a,  d33  [pm/V]  vs  pre-annealing  temperature  [°C],  b.  Image  charge  *105  [C]  vs  pre¬ 
annealing  temperature  [°C],  c.  Relative  concentration  of  Ge  E’  centers  [a.u.]  vs  pre-annealing 
temperature  [°C],.d.  Relative  concentration  of  Si  E’  centers  [a.u.]  vs  pre-annealing 
temperature  [°C 

Figure  12  -  Samples  from  lot  3  that  had  been  dry  pre-annealed  over  a  series  of  temperatures  from  400°C  to 
900°C,  a.  d33  [pm/V]  vs  pre-annealing  temperature  [°C],  b.  Image  charge  *105  [C]  vs  pre¬ 
annealing  temperature  [°C],  c.  Relative  concentration  of  Si  E1  centers  [a.u.]  vs  pre-annealing 
temperature  [°C],  d.  Relative  concentration  of  Ge  E'  centers  [a.u.]  vs  pre-annealing 
temperature  [°C]. 
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Tables 


Table  1  -EPR  parameters  used  for  analysis  of  defect  sites  in  fused  silica  by  the  high  power  second 
harmonic  method 


EPR  parameter 

Value 

Receiver  gain 

1.0  X  105 

Phase 

229° 

Harmonic 

2 

Modulation  frequency 

50  kHz 

Modulation  amplitude 

6.232  G 

Conversion 

81.92  msec 

Time  constant 

163.84  msec 

Sweep  time 

167.772  sec 

Center  field 

3490  G 

Sweep  width 

300  G 

Resolution 

2048 

Frequency 

9.773  GHz 

Power 

20  mW 

lot  4 


x  1.52 

N 
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ABSTRACT 

Laser  produced  plumes  from  a  Ge02  target  were  investigated  using  both  optical  emission 
spectroscopy  and  laser  induced  fluorescence.  Elemental  neutral  and  ionic  Ge  and  diatomic  GeO 
were  monitored  as  a  function  of  02  pressure  at  4  cm  from  the  target  surface.  Data  obtained  from 
laser  fluorescence  experiments  showed  no  indication  of  chemical  production  of  GeO  in  these 
plumes  under  any  conditions  studied;  however,  the  atomic  emission  intensities  of  several  ionic 
and  neutral  Ge  lines  appeared  to  increase  with  02  pressure.  This  signified  plume  interactions 
with  the  background  gas.  The  trends  in  the  species  concentrations  were  determined  and  the 
results  are  interpreted  with  respect  to  mechanisms  for  thin  film  growth. 


INTRODUCTION 

Our  recent  work  in  growing  thin  films  of  germanium  dioxide  (Ge02)  using  laser  ablation 
deposition  showed  that  stoichiometric  thin  films  can  be  fabricated  under  two  different  deposition 
conditions,  and  the  results  of  these  studies  suggest  the  possibility  of  two  separate  kinetic 
mechanisms  for  growing  these  films.  First,  stoichiometric  films  could  be  produced  by  ablating  a 
Ge02  target  in  150  mTorr  of  ambient  02  [1].  At  these  relatively  high  pressures,  elemental  Ge  in 
the  ablation  plume  can  react  gas  kinetically  [2]  with  02  to  form  GeO,  which  can  subsequently 
recombine  with  excess  02  at  the  substrate  surface  in  a  disproportionation  reaction.  Generally,  the 
oxygen  content  in  the  films  increased  with  02  pressure  from  10*5  to  0.15  Torn 

Performing  depositions  at  100's  mTorr  of  an  ambient  gas  causes  the  plume  to  become  very 
highly  forward  directed  [3]  and  leads  to  films  of  nonuniform  thickness.  For  that  reason,  we 
explored  film  depositions  at  low  pressures  to  produce  plumes  that  are  spatially  more  extended 
and  allow  a  better  opportunity  to  produce  flat  films  for  optical  waveguides.  Stoichiometric  films 
were  prepared  at  lower  background  pressures  when  the  ablation  plume  was  passed  through  a  dc 
discharge  in  a  30  mTorr  oxygen  background  [4].  In  these  studies,  we  found  that  the  combination 
of  the  substrate  temperature  and  the  dc  discharge  enabled  us  to  incorporate  the  correct  oxygen 
content  in  the  film.  At  these  lower  pressures,  chemical  reactions  would  appear  unlikely  because 
of  thejow  collision  rates  between  Ge  and  any  ambient  gas.  However,  the  discharge  could  activate 
the  material  in  the  plume  or  the  ambient  gas  providing  additional  energetic  species  to  promote 
surface  recombination  reactions. 

Reactions  between  elemental  species  in  laser  produced  plumes  and  a  reactive  background 
gas  are  known  to  be  important  for  producing  good  quality  high  Tc  superconducting  thin 
films  [5].  An  understanding  of  the  potential  chemical  processes  occurring  between  the  plume 
material  and  any  reactive  gas  in  between  the  target  and  the  substrate  is  needed  to  elucidate  the 
important  thin  film  growth  mechanisms.  Thus,  the  aim  of  this  work  is  to  explore  the  microscopic 
mechanisms  responsible  for  forming  high  quality  dielectric  thin  films  under  both  growth 
conditions  discussed  above. 


347 

Mat.  Res.  Soc.  Symp.  Proc.  Vol.  334.  ©1994  Materials  Research  Society 


EXPERIMENT 

Figure  1  shows  a  schematic  diagram  of  the  experimental  apparatus.  A  puked  ArF jxcimer 

laser  tt  =  193  nm,  pulse  width  of  28  ns  (FWHM)]  was  focused  onto  the  surface  of  a 

laser  la.  iy  ’  F  .  n«itv  of  1  J  cm*2  The  target  was  housed  in  a  chamber  that 

at  normal  incidence  with  an  e^rgy  de^stty  of Ua m  ^  1  8  the  chamber  with  a 

routinely  achieved  a  base  pressure  of  5  x  0  T«r.  < Oxygen ^  ^  ^  path  tQ  determine 

WmiimMifM 

features  originating  from  electronically  excted  output  of  a 

“w  was  vied  between  1  -  »  -  »  “  f 'f"oMA  ^ 

OMA  in  a  single  trace. 


results 

^atefso?ely°f fmm  excited gTItZ  emission  intensi Js ^eSe!"'  m^Ton^s 

« ^,S“*S3,  excited  states!,  the  plame 
decrease  widTutcreln^O^preswrw'as'illtistraied  in  Fig.  J*rhe  '^ray.v'l|e^i  convened  to^  rate 


348 


RING  ELECTRODE 


Figure  1 .  Diagram  of  the  experimental  system.  The  diagram  shows  the  arrangement  for  the  LIF 
studies.  The  set-up  for  the  optical  emission  studies  changes  by  removing  the  dye  laser  beam  and 
operating  the  OMA  in  an  ungated  mode. 
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Figure  2.  Time- integrated  plasma  emission  spectrum  recorded  at  an  oxygen  pressure  of 
30  mTorr  and  at  4.0  cm  from  the  target  surface.  Trace  (a)  shows  Ge  emission  with  (w)  and 
without  (wo)  the  discharge  while  trace  (b)  shows  the  behavior  of  atomic  oxygen. 
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than  that  for  02.  This  may  indicate  that  the  decays  are  not  caused  by  a  chemical  depletion  of 
ground  state  atoms.  Electronic  quenching  could  account  for  these  decays,  but  the  radiative 
lifetimes  of  the  states  under  study  (-250  ns)  prohibit  such  encounters  [7].  For  example,  the  time 
between  collisions  in  an  equilibrium  mixture  is  about  12  jxs  so  the  electronic  states  will  never  see 
a  collision.  These  pressure  experiments  were  repeated  with  the  dc  discharge  with  no  difference  in 
the  results. 

Although  the  emission  studies  did  not  reveal  the  presence  of  GeO*,  the  removal  of  ground 
state  Ge  I  by  reaction  is  not  totally  ruled  out.  Reactions  between  Ge  and  02  produce  GeO  in  the 
electronic  ground  state  with  AE  =  +1.2  eV  [7],  The  exothermicity  of  the  reaction  is  sufficient  to 
populate  GeO(X)  up  to  vu  =  5  so  no  visible  emission  would  be  produced.  We  searched  for 
vibrational  excitation  in  GeO(X),  which  would  indicate  chemical  formation  of  GeO,  by  pumping 
from  various  X  state  vibrational  levels  (0  <  v"  <  5)  to  the  A  state  using  the  strongest  Frank  - 
Condon  allowed  transitions.  We  discovered  strong  emission  lines  from  the  molecule  upon 
excitation  from  only  v"  =  0.  No  LIF  emission  was  observed  after  exciting  from  v"  >  1  even  in  the 
presence  of  various  pressures  of  02.  In  experiments  similar  to  those  performed  with  Ge  I,  the 
GeO  LIF  emission  intensities  decreased  with  increasing  02  pressures.  The  GeO(A  — >  X)  LIF 
signals  became  undetectable  at  02  pressures  near  40  mTorr.  The  radiative  lifetime  for  GeO(A)  is 
not  known,  but,  assuming  typical  values  of  3-5  (is,  a  significant  number  of  collisions  would  not 
occur  eliminating  electronic  quenching  as  the  source  of  the  intensity  decrease.  Although  the 
mechanism  for  these  intensity  decreases  at  various  02  pressures  is  currently  not  understood,  the 
results  do  indicate  that  chemical  production  of  GeO  via  interactions  of  the  plume  material  with 
the  background  02  does  not  occur. 

Since  the  chemical  production  of  GeO  via  plume-ambient  gas  interactions  does  not  appear  to 
be  important,  we  traced  the  origin  of  diatomic  GeO  by  monitoring  LIF  signals  as  a  function  of 
distance  from  the  target  surface.  These  LEF  experiments  showed  that  GeO  molecules  were 
present  in  the  plume  at  distances  as  near  as  2.5  mm  from  the  target  surface.  Molecular  LIF 
emission  was  also  observed  in  vacuum.  These  results  indicate  that  GeO  is  either  produced  as  a 
result  of  recombination  of  material  in  the  plume  near  the  target  surface  where  the  ion  and  neutral 
number  densities  are  high  or  the  molecule  dissociated  directly  from  the  target. 


DISCUSSION 

These  experiments  showed  that  interactions  between  the  plume  material  and  the  ambient  02 
gas  directly  affected  the  concentrations  of  the  excited  neutral  and  ionic  states  of  Ge  and  O  atoms 
as  well  as  molecular  GeO.  Two  possible  mechanisms  can  be  postulated  to  explain  the  observed 
emission  and  LIF  behavior  described  above.  The  first  mechanism  involves  non-reactive  plume  - 
ambient  gas  interactions.  Plume  interactions  with  the  background  gas,  especially  at  the  higher  gas 
pressures  (  >  50  mTorr),  spatially  confine  the  plume  and  should  cause  both  the  atomic  emission 
features  and  the  LEF  intensities  to  increase.  However,  we  observe  increases  in  only  the  atomic 
emission  intensities  and  not  in  the  LEF  signals.  Therefore,  this  mechanism  may  not  be  plausible. 
The  second  mechanism  involves  reactive  encounters  between  the  plasma  and  the  ambient  gas. 
Scattering  events  between  electrons  and  atoms  or  molecules  and  high  kinetic  energy  collisions 
between  atoms  could  produce  neutral  electronic  excitation  as  well  as  additional  ionization. 
Therefore,  atomic  emission  experiments  would  show  an  increase  in  emission  intensities  as  was 
observed  here.  These  interactions  would  also  result  in  a  concurrent  decrease  in  ground  state 
species  which  is  consistent  with  the  LIF  results. 
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Figure  3.  LIF  emission  spectra  from  Ge  as  a  function  of  O2  pressure:  (a)  P(02)  -  25  mTorr,  (b) 
P(02)  =  0  mTorr,  and  (c)  P(02)  =  10  mTorr.  These  spectra  were  recorded  4  cm  from  the  target 
surface  with  a  10.5  \is  delay  time  between  the  excimer  and  the  probe  lasers. 


Figure  4.  Plot  of  the  natural  log  of  the  relative  fluorescence  intensities  vs.  02  pressure.  Decays 
from  two  LIF  emission  lines  are  plotted:  (a)  X  =  326.9  nm  and  (b)  X  =  275.5  nm.  The  cause  of 
this  decay  may  be  due  to  a  combination  of  reactive  and  non-reactive  collisions  with  the  ambient 
oxygen. 
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These  results  may  have  significance  in  determining  the  mechanisms  for  Ge02  thin  film 
growth.  Since  chemical  reactions  do  not  appear  to  dominate  the  plume-ambient  gas  processes, 
surface  reactions  must  control  the  thin  film  growth.  At  high  ambient  gas  pressures  (i.e.,  P(02)  = 
150  mTorr),  excess  02  may  be  required  to  complete  recombination  reactions  on  the  substrate 
surface,  whereas  at  lower  pressures,  the  discharge  produces  excited  O  atoms  which  can  then 
interact  more  energetically  with  the  substrate.  These  processes  produce  the  same  result: 
stoichiometric  thin  films.  The  importance  of  GeO  in  the  film  formation  process  is  not  completely 
understood.  Diatomic  GeO  is  abundant  in  vacuum  but  disappears  as  the  02  pressure  is  increased. 
Since  we  cannot  detect  population  in  other  vibrational  levels  of  the  ground  state,  GeO  may  either 
be  scattered  from  the  excitation  region  or  the  molecule  could  be  dissociated  as  a  result  of  plasma- 
02  interactions. 


SUMMARY 

Constituents  of  a  laser-produced  plume  from  a  Ge02  target  have  been  spectroscopically 
examined.  The  results  indicate  complicated  processes  that  occur  when  the  plume  interacts  with 
the  ambient  02.  High  kinetic  energy  species  are  produced  and  the  behavior  of  individual  species 
in  the  plume  may  be  governed  by  several  reactive  and  non-reactive  encounters.  Chemical 
production  of  diatomic  species,  which  are  critical  for  growing  good  quality  high  Tc  thin  films, 
*does  not  occur  for  laser  ablated  Ge02  plumes  in  02»  and  thus  does  not  appear  to  be  important  for 
growing  stoichiometric  Ge02  thin  films. 
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